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Chapter 1 Introduction 


1.1 About This Book 


This guide describes the 3d MPC software product which is intended for 
multivariable feedback control and optimization of an industrial process that have 
many input and output signals. Inputs in the form of manipulated variables affecting 
actuators, and outputs in the form of process variables from sensors. MPC stands for 
model predictive control, a widely used technique for advanced control. For 
processes with strong interaction between different signals this technique can offer 
substantial performance improvement compared with traditional single-input 
single-output control strategies. 


The 3dMPC product consists of both on-line and off-line components. The main 
function of the on-line components is a model predictive controller. The on-line 
components also include graphical user interfaces for operator interaction and start 
and stop. The off-line engineering tools for configuration, modeling, tuning, and 
analysis are available as a MATLAB toolbox. 


The 3dMPC product is intended to run on a PC under Windows NT. It is connected 
to an OCS/DCS through an OPC server. The graphical interface can run on the same 
PC or on another PC connected through a network. 


If you are a user responsible for installation, configuration, runtime operation or 
maintenance of the 3dMPC, this manual is intended for you. 


To configure and build applications, you are required to have knowledge of 
multivariable control. Basic Windows NT skills are also required. 


This book is not intended to be the sole source of instruction for the 3d MPC 
software product. Those people involved in configuration should attend the 
applicable course offered by ABB 


3BSE 011 572R101 1 


3dMPC User’s Guide 
Chapter 1 Introduction 


Ho - 


Use of Warning, Caution, Information, and Tip 


This publication includes Warning, Caution, and Information symbols where 
appropriate to draw the readers attention to safety related or other important 
information. It also includes a symbol Tips to provide useful hints to the reader. 
The symbols should be interpreted as follows: 


Warning indicates the presence of a hazard which could result in personal 
injury. 


Caution indicates the presence of a hazard which could result in equipment or 
property damage. 


Information alerts the reader to pertinent facts and conditions. 


Tip indicates advice on, for example, how to design your project or how to use a 
certain function. 


Although Warning hazards are related to personal injury, and Caution hazards are 
associated with equipment or property damage, it should be understood that 
operation of damaged equipment could, under certain operational conditions, result 
in degraded process performance leading to personal injury or death. Therefore, it is 
important to comply fully with all Warning and Caution notices. 


1.2 Product Overview 


The product 3dMPC controls a multivariable process by a combination of feedback 
and feed-forward. The control system under consideration is shown in Figure 1-1. 
The controller determines the manipulated outputs, OUT, based on actual 
measurements of process variables, PV, and feed forward signals, FF, and of 
operator defined parameters, and external inputs. The process variables can be 
assigned set-points that are the target for the feedback control law. The feed forward 
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signals are measurable disturbances acting on the process that can be used for feed 
forward 


operator defined 


parameters disturbances 
! f | FF 
external inputs y 
> 
OUT PV 
-—| CONTROLLER ;+———-®| PROCESS 
—-| 


Figure 1-1. Feedback control system. 


The number of variables handled by the controller is limited as follows: 
¢ Number of manipulated outputs (OUT) is limited to 16. 

¢« Number of process variables (PV) is limited to 32. 

¢ Number of feed-forward variables (FF) is limited to 16. 


The 3dMPC product is based on the model predictive controller (MPC) concept. 
This means that, at each sample the manipulated outputs are obtained from 
constrained optimizations. The loss functions in these problems penalizes a 
weighted sum of squared control-errors and moves in the manipulated outputs. The 
constraints in the optimization comes from limits on the manipulated outputs and 
from limits on the process variables. 


The 3dMPC product consist of both on-line and off-line components. The on-line 
components consist of algorithms for control, of a data-logging function, and a 
function for the generation of excitation signals for system identification. The on- 
line components also include graphical user interfaces for operator interaction and 
start and stop of the on-line algorithms. 
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The 3dMPC product is intended to run on a PC under Windows NT 4. It is 
connected to an OCS/DCS through an OPC server running in the same PC as the on- 
line controller. The graphical interface can run on the same PC or on a network PC. 


The off-line engineering tools are implemented as a MATLAB toolbox. There are 
tools for configuration, data treatment, modeling, tuning, and analysis. 


1.2.1 Functional Overview 


3dMPC consists of on-line components for control and operator interaction, and of 
off-line engineering components for configuration, modeling, tuning, and analysis. 


The on-line components of the 3d MPC consist of a controller object and graphical 
interfaces for user interaction. The on-line controller object has an OPC interface 
for data exchange with an OCS/DCS. 


The off-line engineering components are implemented in MATLAB. They are 
accessed using graphical user interfaces. 


The components of 3d MPC is shown in Figure 1-2. The figure also shows the 
interfaces to the OCS/DCS using OPC servers residing in the local NT station. 


1.2.1.1 Model Predictive Control 


Model predictive control is a fairly general class of algorithms for feedback and 
feed-forward control. They have that in common that a model of the process, subject 
for control, is used to predict the effects in process variables due to actual and future 
changes in manipulated outputs and feed forward variables. 


The dynamic model 
y(k) = G(q)u(k) 


that describes the real world process relates how manipulated outputs u(k), from the 
operator or from a controller, affects the actual and future process variables y(k). 
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On-line Components Off-line Components 


Engineering Environment 
User 
Interfaces Configuration Tool 
a Tuning Tools 
On-line Modeling Tools 
Controller 
OPC interface . 
Pre-processing Tool 


OPC server 


Figure 1-2. Overview of components of 3d MPC. 


The process variables are measured signals that should follow desired set-points 
r(k). The model predictive control objective could be defined as follows: 


Find a sequence of current and future increments for manipulated 
outputs that minimizes a weighted sum of future squared control 

errors and a weighted sum of squared increments in the sequence 
of manipulated outputs, while limits for manipulated outputs and 
limits for predicted process variables are considered. 


3BSE 011 572R101 5 


3dMPC User’s Guide 
Chapter 1 Introduction 


This objective can be expressed more precise mathematically: Find a sequence of 
increments for manipulated outputs Au(k+i) over a control horizon of m samples to 
minimize the objective function 


Dp m—1 
id ae 
ie y Pe(k + ill, + y |AAu(k + all 
i=1 i=0 
subject to constraints 
Vig SMR +I SY, Vie [1,p] 
Uj SUK +I) Sup; Vie [0,m-1] 
Au, <Au(k +i) < Au, ; Vie [0,m—1] 


In this minimization the control errors e(k)=r(k)-y(k) are considered over a 
prediction horizon of p future samples. The constraints for process variables y(k) are 
considered over a prediction horizon of p samples, and the constraints for 
manipulated outputs u(k) and their increments are considered over a control horizon 
of m future samples. In this procedure, the dynamic model is used for the prediction 
of future values for the process variables. Important design parameters are the 
weights A and I° To weights used in the 3dMPC product are defined from 
equivalent control error factors and from equivalent move factors. These tuning 
variables define the trade off between control errors and actuator efforts for the 
control system. 


1.2.1.2 Observer Based Model Predictive Control 


In an observer based model predictive controller, like this one, the dynamic model is 
described as a discrete time state-space model 
x(k+ 1) = Ax(k) + B, uk) + By d(k) 
y(k) = Cx(k) 


where x(k) is the state vector, d(k) are measurable disturbances used for feed 
forward. The matrices A, B,, and B, define the model. These parameters are 
determined using the modeling tools provided in the 3d MPC product. 
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1.2.2 On-line Controller Functionality 


The on-line components of the 3d MPC consist of a controller object and graphical 
interfaces for user interaction. 


1.2.2.1 Operational Modes 


The inputs and outputs of the controller can operate in independent modes. This 
means that some signals in the controller can work in automatic mode while other 
signals can work in manual mode. The modes for the signals can be changed on- 
line. 


Operational Modes for Manipulated Outputs 


The operational mode for a manipulated output determines if the actual value is 
determined by the MPC algorithm, or by the operator, or by an external input to the 
controller. They are described in Table 1-1. 


Table 1-1. Operational modes for manipulated outputs. 


Mode Description 
Auto The manipulated output is determined by the control 
algorithm 
Man The manipulated output is defined by an element in the 


manOUT vector, given by the operator. 


E3 The manipulated output is defined by an element in the 
extRef3 external input vector. 


Bal The manipulated output is defined by an element in the 
balRef external input vector. 


A difference between Bal and E3 is that speed limitations are considered in E3 but 
not in Bal. Further, excitation signals for process identification can not be added in 
mode Bal. The way a manipulated output is determined is summarized in 

Figure 1-3. The value of the variable modeOUT determines how the manipulated 
output OUT is determined. A manipulated output in Bal, E3, or Man is normally 
used as a feed-forward variable in the controller. 
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The limits for manipulated outputs are explicitly considered in the control 
algorithm. The achieved value of the addOUT variable is adjusted after the 
saturation so it matches the actual value of the OUT variable. The reason for this is 
to have consistent data in presentations and log-files. 


‘Speed limits: ‘Absolute limits: 

‘External or operator ‘External or operator 
‘supplied, depending + (7. 
‘on configuration i! 


*  {supplied, depending 


excitation signal (i) 
fp \shah cues setae in ose foe Bal 

‘Control | 
‘Algorithm . 
oe eae ayes Speed Amplitude | QUT(i) 
mano Uiti) Limiter Limiter 
extRef3(i) 
balRef(i) 

m= Auto 

m= Man 

m £3 

m= Bal 


Figure 1-3. Overview of how the i:th manipulated output is determined 


Operational Modes for Process Variables 


The operational mode for a process variable determines how the variable is used in 
the controller and how the corresponding set-point is determined. The set-point 
could be determined by the operator or by an external input to the controller. 
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The operational modes for the process variables are described in Table 1-2. 


Table 1-2. Operational modes for process variables. 


Mode Description 

Auto The process variable is controlled to follow a set-point 
defined by the operator. 

E1 The process variable is controlled to follow an external set- 
point. 

Predict The process variable is not controlled. It is used only for 
state estimation in the observer. 

Off The process variable is not used. 


Figure 1-4 shows how the i:th working set-point is determined. 
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autoSP(i) 
extRef1 (i) 


PV(i) 


' Speed limits: 
‘External or operator 
‘supplied, depending 
‘on configuration 


———_ 


—— |) 
1 


m= Auto 
mm F1 


= Predict, Off 


Table 1-3. Operational modes for feed-forward variables. 


Speed 
Limiter 


‘Absolute limits: 
‘External or operator 
‘supplied, depending 


Amplitude 
Limiter 


Figure 1-4. Overview of how the i:th working set-point is determined. 


Operational Modes for Feed-Forward Variables 


WSP(i) 


The operational mode for a feed-forward variable determines how the variable is 
used in the controller. The operational modes for the feed-forward variables are 
described in Table 1-3. 


Mode Description 
Active The feed-forward variable is used for feed-forward. 
Predict The feed-forward variable is used only for state estimation 
in the observer. 
Off The feed-forward variable is not used. 
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Critical Variables 


Any variable in the controller can be marked as critical. This means that the variable 
is So important that it must take active part in the MPC algorithm, otherwise it is not 
meaningful or allowed to use automatic control at all. The controller will then have 
all manipulated outputs operating in manual mode. 


Three Degree of Freedom Dynamic Optimization 

The controller can be configured to have different dynamic responses to 
* — set-point changes 

e changes in measurable disturbances 

¢ feedback control. 


To achieve this, three optimization problems, of the type described in Section 
1.2.1.1, Model Predictive Control, are solved with appropriate control error 
formulations. 


Static Optimization 


After the dynamic optimizations has been performed a static optimization algorithm 
is employed to drive the manipulated outputs towards desired values. This 
algorithm will only have effect in cases where there is ambiguity in how 
manipulated outputs are chosen for desired targets for the process variables and 
when desired values are supplied. The static optimization algorithm penalizes 
weighted differences between the actual and desired values for manipulated outputs. 


Constraint Handling 


A main feature of model predictive controllers is the ability to handle constraints, 
not only for manipulated outputs but also for process variables. 


Two sets of soft constraints with priorities can be defined for each process variable: 
One set defines constraints as absolute values, and one set defines constraints as 
values relative to the set-points. The constraints could be defined either by external 
inputs or by the operator. 


Soft constraints are defined in terms of priorities as described in Table 1-4. 
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Table 1-4. Definition of soft constraint 


Priority Description 
1-9 Nine levels for highly penalized constraints. Priority 1 is the 

highest. 

10 The penalty for constraint violation is 10 times the penalty on 
feedback control errors. 

11 The penalty for constraint violation is 5 times the penalty on 
feedback control errors. 

12 The penalty for constraint violation is twice the penalty on 
feedback control errors. 

13 The penalty for constraint violation is the same as the penalty on 
feedback control errors. 

14 The penalty for constraint violation is half the penalty on feedback 
control errors. 

15 The penalty for constraint violation is 1/5 of the penalty on 
feedback control errors. 


Constraints can be defined for each manipulated output: There is one set of hard 
absolute constraints, and one set of hard speed constraints. There is also one set of 
soft constraints with priorities that can be used to define the most desired 
operational range for manipulated outputs. The constraints can be defined either by 
external inputs or by the operator. Table 1-5 describes possible modes for 


constraints. 
Table 1-5. Operational modes for constraints. 
Mode Description 
Ext The constraint is defined by external input. 
Opr The constraint is defined by the operator. 
Off The constraint is not used. 
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The operational mode of each constraint can be individually configured on-line, 
either to be defined by an external input or by the operator, or not to be accounted 
for at all. 


Each soft constraint is given a priority that defines how important it is to satisfy a 
particular constraint in a situation where not all soft constraints can be satisfied. In 
such situations the optimization algorithm provides a control signal that minimizes a 
weighted sum of squared constraint violations. The weights in this optimization 
procedure depend on which of the fifteen priority levels that has been chosen. 


Signal Linearization 


The controller integrates functions for non-linear transformations on each signal. 
This is useful when the process is composed of a linear part and non-linear parts as 
shown in Figure 1-5. There is one non-linear function, ,;, for manipulated 
outputs, one non-linear function, Yr, for feed-forward variables, and one non-linear 


function, Yow for the process variables. . 


Linear Dynamics PV 


eg Control 
oun Algorithn ~“—_ Yr <——— 
WSP 
<< 


<q] (pv)! 


Figure 1-5. Linearization function 
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The functions ‘Pou, ‘Pre, and ‘¥,, have diagonal structures. Hence they handle non- 
linearities for each individual variable. The functions ‘Poy, Pr, and ‘¥,, are strictly 
monotonic increasing or decreasing functions in order to be invertible. The non- 
linearities are defined in the pre-processing tool, prior to filtering 


Parameter Scheduling 


The controller can be configured to use an external input signal for parameter 
scheduling. One of four pre-defined sets of controller parameters is then used in the 
algorithm of the controller. The set of controller parameters accounted for in the 
scheduling involves for example a complete process model, equivalent error and 
move intervals, and limit definitions. 


Hysteresis bands are used to avoid unnecessary switching between two adjacent 
intervals. The hysteresis band between two intervals is defined as 10% of the 
smallest of the two intervals. 


When the parameter scheduling is disabled the operator can manually select one of 
four pre-defined specification files for the operation of the controller. 


On-line Design Calculation 


The controller can be reconfigured on-line. This happens either when a new set of 
controller parameters is loaded or when the operational mode of a signal is changed. 


PRBS Generator 


The controller has functions for generation of excitation signals for process 
identification. Pseudo Random Binary Sequences can simultaneously be added to 
several manipulated outputs. This integrated function provides a user with means to 
easily do experiments for process identification. 


Data Logging 


The controller has a data logging function that stores selected signals and variables 
in a log-file. This integrated function provides a user with means to easily do 
experiments for process identification. 
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Backup, Restore, and Warm-start 


A backup function allows that data, defined by the operator, can be stored on disk 
for later retrieval by the restore function. 


A backup can also be used for warm-start of the on-line controller. In such a case the 
controller is loaded with the stored data direct at start up. Backups for warm-starts 
could be done on a request from the operator, or after any parameter subject for 
backup has changed. With this approach the on-line controller could be started with 
well defined parameters, or with the latest used parameters. 


The on-line controller will have all manipulated outputs operating in manual after a 
restore or a warm-start. 


1.2.3 On-line Graphical User Interfaces 


1.2.3.1 Interface for start-up 


The start-up interface is used to point out an application that will be controlled. The 
application has earlier been prepared by the off-line engineering tools. The start-up 
interface is also used to select the OPC-server that will be used. 


1.2.3.2 Operator Interface 


The operator interface presents actual data for the on-line controller and allows the 
operator to set parameters in the on-line controller. 


It can run on the same node as the on-line controller or it could run on any NT 
workstation in the local area network. 
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1.2.4 Off-line Engineering Tools 


The engineering tools are available in MATLAB. MATLAB is launched with the 
proper settings from the NT Start-Menu. The different engineering tools are then 
accessed using pull-down menus in the engineering tool main window that is 
available when MATLAB has been started. The following off-line engineering tools 
are available: 


¢ Configuration Tool 
e  Pre-processing Tool 
* Modeling Tools 

¢ Tuning Tool. 


These are described in the following sections. They have all advanced graphical 
user interfaces. 


1.2.4.1 Configuration Tool 


The configuration tool is used to create and modify the actual configuration of the 
control system. A configuration has to be defined before any of the other parts of the 
3dMPC is used. The configuration is a definition of what signals in the OCS/DCS 
that are used as manipulated outputs, as process variables, and as feed-forward 
variables. 


1.2.4.2 Pre-Processing Tool 


The pre-processing tool provides functions for signal processing to generate data 
sequences suitable for the process identification procedure in the modeling tools. 
Data sequences are initially created from log-files that contain data that has been 
recorded by the on-line controller logging function. 


The functionality available for signal processing is described below. Properties of 
the data sequences are shown graphically. The plots are updated after each action. 
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Linearization 


Static non-linearities can be introduced for each individual signal to compensate for 
non-linear sensors and non-linear actuators. These non-linear functions are 
preserved in controllers based on models built from the data sequences where the 
non-linear functions originally were defined. 


In the pre-processing tool the non-linear transformations of a data sequence are 
performed prior to the linear filtering. 


Linear Filtering 


A data sequence can be filtered by a low-pass, a high-pass, a band-pass, or by a 
band-stop filter. The filtering function also involve removal of mean-values and 
trends. 


Data Editing 


The data editing functionality allows removal of outliers in the data. A cut function 
removes a part of a data sequence, a split function splits a data sequence in two 
parts, and a merge function merges data sequences. There is also a function to re- 
sample a data sequence with a new sampling interval. 


1.2.4.3 Modeling Tools 


The modeling tools are used for building models of the process from data 
sequences. The models are intended for the controller tuning described in Section 
1.2.4.4, Tuning Tools. There are modeling tools for process identification and for 
merging of models obtained from different process identifications. The tools also 
contain elaborate analysis functions to assist the user. The modeling tools are 
accessed using advanced graphical user interfaces. 


Parametric Identification 


The main part of the modeling tools are the parametric identification tool that is 
used for building state-space models from data sequences. The used two step state- 
of-the art algorithm for parametric identification combines ease of use with the 
ability to use data sequences from both open and closed loop identification 
experiments. 
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Model Merging 


There is a tool to merge models created at different occasions and covering different 
parts that is containing different sets of inputs and outputs of a process. 


Given two models, a new model is constructed where each transfer function is taken 
from either of the two models. The user selects the used parts. 


The main use of this tool is to combine identified models with simple user defined 
models. 


Definition of Simple Models 


This tool allows a user to define models that are already known. For each 
component of a multivariable transfer function, the user can enter a gain, a time 
constant and a delay. 


Model Evaluation 


In addition to the model creation tools, there are also tools to examine the models in 
different respects. The following analysis functions are available: 


¢ Simulated/Predicted output 


The measured outputs of a given evaluation data set are plotted together with 
either the simulated or the predicted model outputs. 


° Prediction Error Variances 


The prediction error variances as a function of prediction length are plotted for 
given models and a given evaluation data set. 


¢ Noise Spectrum 
The noise spectrum for each output of given models are plotted. 
e Frequency Response 
Gain and phase curves for all transfer functions of given models are plotted. 
¢ — Transient Response 
Step or impulse responses for all transfer functions of given models are plotted. 
¢  Poles/Zeros 


The poles and zeros for all transfer functions of given models are plotted. 
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1.2.4.4 Tuning Tools 


The tuning tools are used for the controller tuning. The tools also contain elaborate 
analysis functions to assist the user. 


MPC Tuning 


This tool is used to create and modify data for the controller tuning. 


MPC Tuning Analysis 


The MPC tuning analysis tool is used to evaluate the behavior of the feedback 
control system. The following functions for analysis are available for all relevant 
open and closed loop transfer functions: 


Step responses 
Bode diagrams 
Nyquist diagrams 


Singular value frequency response plots. 


Special emphasize is put upon transfer functions that are relevant for robustness 
considerations. 


1.2.5 Product Release History 


Table 1-6 lists the major milestones in the development of 3d MPC. 


Table 1-6. Release History 


Version Description 


1.0 


Internal prototype. 


1.1 


First official release. 
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1.3 Prerequisites and Requirements 


The following minimum recommendations apply for the 3dMPC product: 


Personal computer (PC) Pentium IT 350 MHz 


Primary memory 128 Mb 
Operating system Windows NT version 4 with service pack 4 (SP4) or later. 
Monitor 19 inch, 256 colors, 1280x1024 pixels resolution 


An OPC server must be available for on-line control. 


The 3dMPC product runs under Windows NT4, with SP4. The engineering tools 
require MATLAB 5.2 or later. 


It is a requirement that the selected personal computer must be designed to 
withstand climatic, mechanical and electrical effects in the intended environment. 


1.4 Related Documentation 


Table 1-7 lists documentation related to the 3d MPC product. 


Table 1-7. Related Documentation 


Category Title Description 
Software Advant OPC Server for Advant Fieldbus | This book provides instructions for using 
100 User’s Guide the OPC server for Advant Fieldbus 
100. 
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The following is a list of associated abbreviations that you should be familiar with. 


Term 


AC 


AF100 


COM 


DCOM 
FF 


GUI 
MIMO 
MPC 


OCS 
OLE 


OPC 
OUT 
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Description 


Advant Controller, the common term for the series of 
controllers used in ABB automation systems, for 
example AC110, AC450. 


Advant Fieldbus 100. AF100 is a high speed 
communication link intended for communication 
between controllers in the AC 400 and AC 100 series 
and other equipment adapted for the bus. 

The transmission rate is 1.5 Mbit/s and it supports three 
transmission media: twisted pair, coaxial and optical 
media. 


Component Object Model, a specification that defines 
how individual software components can interact and 
share data under Windows. Developed by Microsoft. 


Distributed COM. Extends COM to networks. 


Feed-forward variables in the controller. These signals 
are measurable disturbances that are used for feed- 
forward. 


Graphical user interface 
Multi-Input Multi-Output 


Model predictive control. The type of control algorithm 
used in 3dMPC 


Open Control System 


Object Linking and Embedding. A technology, based on 
COM, developed by Microsoft. 


OLE for Process Control. 


Manipulated outputs from the controller. These are 
inputs to the controlled process. 
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Term Description 
PRBS Pseudo Random Binary Sequence 
PV Process variables in the controller. These are outputs 


from the controlled process. 


The following is a list of associated terms that you should be familiar with. 


Term Description 


Advant OCS The term refers to one or more of the open control 
systems from ABB, that is, ABB Advant, ABB Master 
and ABB MOD 300. 


MATLAB An interactive system for advanced technical computing, 
developed by The MathWorks. 
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Chapter 2 Installation 


2.1 Site Planning 


When you are site planning, it is important that you check that requirements for 
environmental conditions (temperature, humidity, chemical), grounding and power 
supply are met, according to standards for your application. These requirements are 
stated in the hardware documentation for your PC. 


2.2 Installation Procedures 


Before you start the installation make sure that you log in as a user with 
administration privileges. 


Run the 3d MPC Setup.exe program from the CD if it does not start automatically. 


1. 
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Welcome 
Click Next. 
Software License Agreement 


Please read the license agreement. 
If you accept all the terms of the preceding license agreement, click Yes. 
If you click No, the Setup will close. 


User Registration 

Please register your product now, click Next. 

Registration Confirmation 

Confirm that the registration information is correct, click Yes. 
If you click No, the Setup will close. 

Choose Advant Base Directory Destination Location 
Choose default directory. 

Click Next. 
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6. 3dMPC 
Select the options you want to install. 
Click Next. 
7. ABB Advant Setup Log 
Gives information about the configuration settings. 
Click Next. 
8. Restart Windows 


Choose Yes or No. If new components were added to your system it is 
necessary to restart Windows for a proper operation of the 3dMPC product. 


Click OK. 


If the MATLAB executable is installed in another location than 
C:AMATLAB\bin\matlab.exe the Start Menu shortcut Advant > 34MPC > 
Engineering Tools must be modified to point out the actual location of the MATLAB 
executable. 


2.2.1 Environment Variable 


The environment variable MPCApplicationsRoot defines the directory in which 
each application has a separate directory tree with all its data. 


In the installation, MPCApplicationsRoot is defined to point out a directory where 
the example applications are located. 


The variable MPCApplicationsRoot can be changed to point out another location 
for the applications. 


2.2.2 License Installation 
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The product 3d MPC requires valid license files to work. It needs one license file for 
the engineering tools and one license for each application. The license files are 
supplied on floppy disks. The floppy also contains an installation program. 


License files are installed after the software is installed. 
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The engineering tool license is installed by the command 
>> A:\lic_e 


It places the engineering license in the directory 3MPCApplicationsRoot§, that 
is the directory that is pointed out by the environment variable 
MPCApplicationsRoot. 


The on-line controller license is installed by the command 
>> A:\lic_a ApplicationName 


The name ApplicationName is the name of the application for this license. The 
license file is placed in the directory 
SMPCApplicationsRoot%\ApplicationName. 


2.3 Uninstallation Procedures 
Uninstallation removes the installed product from your machine, that is it removes: 
e¢ Files 
e« Registry entries 
e  Start-menu icons 
e Unregister ActiveX components 
To run uninstallation you do as follows: 
1. Stop the 3dMPC operator interface clients and the on-line controllers. 
2. Open the Control Panel 


3. Double-click the icon Add/Remove Programs to view the Add/Remove 
Programs Properties window, see Figure 2-1. 


4. Select the product you want to uninstall and press the button Add/Remove.... 


A window asking if you are sure to uninstall the specified product will appear. 
Click Yes if you are sure, otherwise click No. 


6. Wait for the uninstallation to finish. 
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Add/Remove Programs Properties AES 
Install/Uninstall Windows NT Setup | 


es) To install a new program from a floppy disk or CD-ROM 
—_ 


drive, click Install. 


(22) The following software can be automatically removed by 

Windows. To remove 4 program or to modify its installed 
components, select it from the list and click 
Add/Remove. 


AC 450 AMPL Libraries 1.2/1B2/0 
Adobe Acrobat 3.01 

Adobe Acrobat PDFWriter 3.03 
Adobe FrameMaker v5.5 


Bus Configuration Builder 2.6/0B5/0 
Compag Insight Management Agents hd 


Add/Remove... | 
OK | Cancel | Apply | 


Figure 2-1. Add/Remove Programs Properties 


2.4 Start-up and Shut-down Procedures 


The functions of 3d MPC are implemented as ordinary Windows NT applications. 
They are not automatically started when Windows is started. They are available 
from the windows Start Menu: Advant > 3dMPC. 


The functions of 3d MPC are stopped on request from the user. 
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2.5 Product Verification 


After you have performed the installation of the 3d MPC product, you are 
recommended to verify that the software has been properly installed. This means 
that all components could be started and that the on-line controller can connect to an 
OPC server. 


If the option “3d MPC examples” has been installed and if the environment variable 
MPCApplicationsRoot is unchanged since the installation, the following steps 
verifies the installation: 


1. 
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Check that the 3dMPC Engineering Tool is available from the Start menu. 
Select Advant > 3dMPC > Engineering Tool from the Start menu. This 
should start MATLAB, and when you are prompted for an application, select 
PDistiller. You should now see the window “Engineering Interface: PDistiller’. 
Exit by File > Exit. 


Check that the 3d MPC On-Line controller is available from the Start menu. 
Select Advant > 3dMPC > On-Line Controller from the Start menu. You will 
now see the “3dMPC Start-Up and Shut-Down” window. Select the PDistiller 
application, then check the simulate process box, and click Start. In the “Enter 
password for the controller” window, click OK. You should get the message 
“The PDistiller is OK!”, and a DOS window appears where license information 
is shown. Let the application run. 


Check that the 3dMPC Operator Interface is available from the Start menu. 
Select Advant > 3dMPC > Operator Interface from the Start menu. Select 
Local Machine and click Find Controllers. Select the PDistiller application, 
and click Connect. The Operator Interface should now appear. Exit by 

File > Exit. 


Check that the 3d MPC On-Line controller can find OPC servers. Go back to 
the “3d MPC Start-Up and Shut-Down” window. Click Stop and confirm by 
clicking Yes. Remove the check-mark from the simulate process box. 

Click Start. You will now see the “Select OPC Server” window. Verify that 
you can see the available OPC servers in the pop-up menu. Click Cancel. Then 
click Close in the “3d MPC Start-Up and Shut-Down” window. 
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2.6 Uninstallation Verification 
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After you have performed the uninstallation of the 3d MPC product, you are 
recommended to verify that the software has been properly uninstalled. To 
accomplish this, do the following: 


1. Verify that 3d MPC is not shown in the Add/Remove Programs Properties 
window that is launched from the Control Panel. See Figure 2-1. 


2. Verify that the directory SystemModules/3dMPC is empty. It is located under 
Advant Base Directory Destination Location that was chosen during 
installation, see Chapter 2.2, Installation. 
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3.1 Considerations 


This chapter describes how multi-variable controller applications are built for 
3dMPC. Building an application could follow the scheme outlined in Figure 3-1. 
Control system design is in many perspectives an iterative process which is 
illustrated by the revise directives in Figure 3-1. The user is assumed to have 
knowledge of necessary revisions when models or tuning specifications not are 
satisfactory. The dashed arrows in Figure 3-1 indicate data usage, and should not be 
confused with the workflow indicated by solid arrows. 


3.2 Capacity and Performance 


The 3dMPC product can be configured to handle various numbers of signals. 
Table 3-1 shows the maximum number of each signal type that can be handled. 


Table 3-1. Maximal number of signals 


Type of Signal Maximum number 
Manipulated outputs 16 
Process variables 32 
Feed-forward variables 16 


The memory consumption and the execution effort for the on-line controller depend 
mainly on the number of signals, on the complexity of the process model, and on 
design parameters as prediction horizon and control horizon. 


dd) Make sure that the PC is powerful enough to execute the on-line controller with the 
desired sampling interval for the most demanding tuning specification. 
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Experiment design 
On-line controller 

- open/closed loop 

- PRBS amplitude 

- PRBS sampling int. 


Data pre-processing 
Pre-Processing Tool 

- remove outliers 

- static non-linearities 

- remove mean and trend 
- filtering 


Subspace identification 
Modelling Tool 
- select model set 


Configuration 
Configuration Tool 
- signal selection 

- sampling interval 


Model evaluation 
Modelling Tool 

- prediction errors 

- frequency responses 


4 / OK 
/ 
Y 7 
Prediction error model 
identification / 


Modelling Tool 
- starting model 


- simulation 


Revise 


A 


Not OK 


Revise 


Model evaluation 
Modelling Tool 

- prediction errors 

- frequency responses 


WW 


MPC tuning 
Tuning Tool 

- model 

- error intervals 

- move intervals 
- observer design 
- constraints 


> Tuning specifications 


OK 


OK, ready for commissioning 


- simulation 


Not OK 


a 


MPC analysis 
Tuning Tool 

- simulation 

- frequency responses 


Revise 


a 


Not OK 


Revise 


Full simulation 
On-line Controller 
(Internal simulator) 


Not OK 


Ee 


Figure 3-1. Outline for building a 3dMPC application. 
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3.2.1 Memory Usage 


Figure 3-2 shows the estimated memory usage as function of the problem 
complexity. The solid line shows the estimated memory usage for continuous 
operation. The dashed line shows the estimated peak memory usage that occurs 
during an on-line redesign of the controller. 


The complexity parameter n, is given by 
My = Noy M+Nyy* P 


where 71,,,; is the number of manipulated outputs, m is the control horizon, n 


pv 18 the 
number of process variables, and p is the prediction horizon. 


Memory usage, MB 


o 66006100 «150 «200 280 300 350 400 450-500 
Problem complexity, n, 


Figure 3-2. Estimated memory usage. 
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3.2.2 Execution Time 


The execution time per sample, for a certain hardware configuration, depend on the 
problem complexity. Figure 3-3 shows the estimated execution time per sample as 
function of the problem complexity for a PC with a 400 Mhz Pentium II processor 


The complexity n> is given by 


Ay = Noy NTN, +Ngy: P 


and n,; is the number of manipulated outputs, m is the control horizon, n, is the 
total number of soft constraints for process variables, deviations, and manipulated 
outputs, 14, is the number of dead-bands used for process variables, and p is the 
prediction horizon. 


Execution time, seconds 
co 


Q 50 100 160 200 250 
Problem complexity, nj 


Figure 3-3. Estimated execution time. 
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3.3 Engineering Tools Overview 


The engineering tools are started from the windows start-menu. That is, click the 
Start button, then click on Advant, then click on Engineering Tools as shown in 
Figure 3-4. 


Documentation 
LS} Engineering ey Engineering Tool 
Compaq Information Center > (5) MIMO Development Project & On-Line Controller 


ro} Software Administration & Operator Interface 


Netscape SmartUpdate 
(3) WinSDE 1.21 


Figure 3-4. Engineering interface initialization from the windows start menu. 


Then a window appears, where an MPC-application is chosen from a pop-up menu 
with the available MPC-applications on this PC. The engineering tools control panel 
then appears as shown in Figure 3-5. 


Engineering Interface: PDistiller |. (OP x} 
File Tools AMPL ‘Windows Help 
Licensed to: Development Team, SEAPR 


Figure 3-5. Engineering Tool control panel. 


The available MPC-applications on a PC are located in the MPC-application root 
directory, that is defined by the Windows NT environment variable 
MPCApplicationsRoot!. The MPC-application root directory contains one 
directory for each MPC-application. An MPC-application is recognized by the 
name of this directory. See Section 5.2.2, Application Directory Structure, for 
details. 


3.4 Tutorial for Configuration and Application Building 


A tutorial is found in Appendix B. 


1. This variable is set by the installation procedure. 
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3.5 Application Procedures 


3.5.1 Configuration Tool 


A configuration has to be defined before any of the other parts of the controller is 
used. The configuration data consists of parameters describing the controller 
structure, for example number of signals, OCS tag-names for object connections, 
signal ranges and so on. The configuration data also contains variables with 
symbolic names that are used to identify the different signals in the controller. The 
use of these variables makes it possible to use the controller and the models in other 
configurations with the same symbolic names. 


The Configuration Tool is basically an editor for the configuration data. The 
configuration data is stored in the file config.mxc in the application root directory. 
The Configuration Tool consists of a main configuration window and one window 
for each type of signal, that is manipulated outputs, process variables and feed 
forward variables. These windows are used to modify properties of particular 
signals. 


3.5.1.1 Main configuration window 
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Figure 3-6 shows the main configuration window that gives an overview of all the 
signals in the controller, the basic sampling interval, the execution mode and the 
parameter scheduling. 


The user defines the following configuration data directly from the main 
configuration window: 


¢ Description of configuration (editable field) 

¢ Basic sampling interval (editable field) 

¢ Units for basic sampling interval (pop-up menu) 

e Execution mode (pop-up menu) 

e« Name of parameter-scheduling variable (editable field) 


¢ Description of parameter-scheduling variable (editable field) 


¢ Units for parameter-scheduling variable (editable field) 
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Figure 3-6. Main configuration window 


¢« Connection of parameter-scheduling variable (editable field) 
e¢ Range for parameter-scheduling variable, low and high (editable fields). 


The basic sampling interval, units for basic sampling interval and execution mode 
must be filled in. If no values are entered in these fields a warning will be displayed. 
Optional parameters are marked Opt. Basic sampling interval and parameter 
scheduling ranges are numeric values, if a character string is entered in theses fields 
a warning will be displayed 
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All changes for each separate signal is invoked from pull-down menus in the main 
configuration window. The pull-down menu functions are summarized in Table 3-2. 


Table 3-2. Pull-down menus for the main configuration window 


Menu Menu item Function 
File Save Save all the data to config.mxc. 
Exit Close the window and exit. 
Edit PV Select a process variable from a sub-menu and 


open a process variable configuration window for 
entering properties for the selected process 
variable. 


OUT Select a manipulated output from a sub-menu and 
open a manipulated output configuration window for 
entering properties for the selected manipulated 
output. 


FF Select a feed forward variable from a sub-menu and 
open a feed forward variable configuration window 
for entering properties for the selected feed forward 
variable. 


Insert PV Select a process variable from a sub-menu after 
which a new process variable is inserted. 


OUT Select a manipulated output from a sub-menu after 
which a new manipulated output is inserted. 


FF Select a feed forward variable from a sub-menu 
after which a new feed forward variable is inserted. 


Delete PV Select a process variable from a sub-menu for 
removal. 
OUT Select a manipulated output from a sub-menu for 
removal. 
FF Select a feed forward variable from a sub-menu for 
removal. 
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Table 3-2. Pull-down menus for the main configuration window (Continued) 


Menu Menu item Function 
Move PV Select a process variable from a sub-menu that 
should be moved. 
OUT Select a manipulated output from a sub-menu that 
should be moved 
FF Select a feed forward variable from a sub-menu that 
should be moved. 
Windows | Existing MATLAB | Enables easy switch to another MATLAB figure 
figures 
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3.5.1.2 Manipulated output configuration windows 
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shows a manipulated output configuration window that is used for 
modification of configuration parameters for a selected manipulated output. 


a Configuration Edit: OUT1 


Figure 3-7. Manipulated outputs configuration window. 


The user enters data in the editable fields. The general part contains symbolic name, 
description, unit, and ranges for this signal. The connections part contains OCS tag- 
names for signals connected to the controller. OCS tag-names are entered as strings. 
The format of the string depends on the naming convention used by the OPC server 
the controller will be connected to. 
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Name and range in the general part must be filled in and a connection must be 
defined for the OUT attribute. Optional parameters are marked Opt. 


Range parameters are numeric values, if a character string is entered in theses fields 
a warning will be displayed. 


3.5.1.3 Process variable configuration windows 


Figure 3-8 shows a process variable configuration window that is used for 
modification of configuration parameters for a selected process variable. 


Configuration Edit: P¥1 


Figure 3-8. Process variable configuration window. 
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The user enters data in the editable fields. The general part contains symbolic name, 
description, unit, and ranges for this signal. The connections part contains OCS tag- 
names for objects connected to the controller. OCS tag-names are entered as strings. 
The format of the string depends on the naming convention used by the OPC server 
the controller will be connected to. The supervision part contains limits that are used 
for the process variable supervision function. 


Name and range in the general part must be filled in and a connection must be 
defined for the PV attribute. If no values are entered in these fields a warning will be 
displayed. Optional parameters are marked Opt. 


Range and supervision parameters are numeric values, if a character string is 
entered in theses fields a warning will be displayed 
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3.5.1.4 Feed-forward variable configuration windows 


Figure 3-9 shows a feed-forward variable configuration window that is used for 
modification of configuration parameters for a selected feed-forward variable. 


Configuration Edit: FF1 


Figure 3-9. Feed forward variable configuration window. 


The user enters data in the editable fields. The general part contains symbolic name, 
description, unit, and ranges for this signal. The connections part contains OCS tag- 
names for objects connected to the AdvaControl MIMO Controller. OCS tag-names 
are entered as strings. The format of the string depends on the naming convention 
used by the OPC server the controller will be connected to. The supervision part 
contains limits that are used for the feed-forward variable supervision function. 


Name and range in the general part must be filled in and a connection must be 
defined for the FF attribute. If no values are entered in these fields a warning will be 
displayed. Optional parameters are marked Opt. 
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3.5.2 Pre-processing Tool 


Range and supervision parameters are numeric values, if a character string is 
entered in theses fields a warning will be displayed 


This section describes the user interface of the Pre-processing Tool of 3d MPC. In 
addition to the material here, Section B.4, Data pre-processing in Chapter B, 
Configuration Tutorial contains an example of the use of the tool. 


The purpose of the Pre-processing Tool is to provide functions for signal processing 
to generate data sequences suitable for the process identification procedures in the 
Modeling Tools (Section 3.5.3, Modeling Tools). 


3.5.2.1 Overview of the Pre-processing Tool 
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The Pre-processing Tool contains functions to perform the following tasks: 


Filtering: 


De-mean/De-trend: 


Signal Close-up: 


Cut Data: 


Split Data: 


Resample: 


Apply linear filters on the data, Section 3.5.2.2, How 
to filter data. 


Remove trends and/or the mean values from a data 
sequence, Section 3.5.2.3, How to remove mean and 
trends from data. 


Provides an enlarged view of one of the signals in a 
data sequence, Section 3.5.2.4, Overview of the Signal 
Close-up function and Section 3.5.2.5, How to use the 
Signal Close-up function. 


Remove parts of a data sequence, Section 3.5.2.6, 
Overview of the Cut Data and Split Data Tools and 
Section 3.5.2.7, How to use the Cut Data Tool. 


Split a data sequence in two parts, Section 3.5.2.6, 
Overview of the Cut Data and Split Data Tools and 
Section 3.5.2.8, How to use the Split Data Tool. 


Resample a data sequence Section 3.5.2.9, How to 
resample data. 
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Non-linear Transformation: Create or edit static non-linear transformations (used 
in the Pre-processing Tool and the On-line Controller) 
Section 3.5.2.10, Overview of the Non-linear 
Transformation Editor and Section 3.5.2.11, How to 
use the Non-linear Transformation Editor. 


Remove Outliers: Remove outliers from the data, Section 3.5.2.12, 
Overview of the Remove Outliers Tool and Section 
3.5.2.13, How to use the Remove Outliers Tool. 


Edit description and history: Edit the data descriptions and history logs saved 
together with data sequences, Section 3.5.2.14, How to 
edit data descriptions and history logs. 


Log-file Import: Convert a log-file to a data sequence file, Section 
3.5.2.15, Overview of the Log-file Import Tool and 
Section 3.5.2.16, How to use the Log-file Import Tool. 


Merge Data: Merge up to ten data sequences into one, Section 
3.5.2.17, Overview of the Merge Data Tool and 
Section 3.5.2.18, How to use the Merge Data Tool. 


The tool is started from the Tools > Pre-processing menu of the Engineering Tools 
control panel (see Section 3.3, Engineering Tools Overview). The main window of 
the tool is shown in Section 3-10, Overview of the Pre-processing Tool with 
descriptions of the items in Table 3-3 and of the menus in Section, The menu bar of 
the Pre-processing Tool main window on page 48. 


Each data file contains both the original data sequence as contained in the log file 
and the processed data. The processing is always performed on the original data, so 
it is not possible to apply a series of filters by opening a processed file and continue 
the processing. 


The following steps are performed each time a specification is changed in the tool 
(only the signals currently displayed are influenced) or when a file is saved (all 
signals are influenced). 


1. Apply any non-linear transformations. 


2. Remove mean and trend (if specified). Any specified trend break-points (see 
Section 3.5.2.3, How to remove mean and trends from data) and data borders 
between merged data sets (see Section 3.5.2.17, Overview of the Merge Data 
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Tool and Section 3.5.2.18, How to use the Merge Data Tool) divides the data 
sequence into smaller parts, in which the mean and trend values are removed 
separately. 


Filtering (if specified). Data set borders divides the data sequence into smaller 
parts, which are filtered separately. 


In every data-file there is a history log, in which events and important manipulations 
of the data are logged. The following functions are possible to trace in the log: 


Save and save as 

Cut Data 

Split Data 

Changed non-linear transformations 
Remove Outliers 
De-mean/De-trend 


Filtering. 
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¢ MIMO Pre-processing 


GCE ese 
File Edit Tools Windows 


60 80 
Spectrum (rad/s) 


Sampling interval: [CER nin Fiter order: [ER 


Figure 3-10. Pre-processing Tool main window. 
Descriptions of the items are given in 


3BSE 011 572R101 45 


3dMPC User’s Guide 
Chapter 3 Configuration/Application Building 


46 


Table 3-3. Descriptions of the items of the Pre-processing Tool 


main window, Figure 3-10 


Item 


Description 


The menu bar. Described in Section, The menu bar of the Pre- 
processing Tool main window on page 48. 


Small plot areas showing the power spectra of signals. The left 
column of plot areas display Process Variables and the right column 
display Manipulated Outputs and Feed-Forward variables. Five 
Process Variables and a total of five Manipulated Outputs and Feed- 
Forward variables can be displayed at the same time, which ones 
selected by the pop-up menus above the plot areas (item 3). Both the 
raw\') and filtered data are displayed. One pair of Process Variables 
and Manipulated Outputs/Feed-Forward variables can also be 
displayed in the larger plot areas in the left part of the window (items 
6 and 7). To select a signal to be displayed in the larger plot area, the 
small plot area displaying the signal is clicked with the mouse. The 
background in the small plot will change from black to grey to indicate 
the selection. 


Pop-up menus where the signals to display in the small plot areas 
(item 2) are selected. The left menu contains all Process Variables 
contained in the data sequence and the right menu contains all the 
Manipulated Outputs and Feed-Forward variables. 


Large plot areas showing the time series of a pair of Process 
Variables and Manipulated Outputs/Feed-Forward variables selected 
in the small plot areas in the right part of the window (item 2). Both the 
raw'") and filtered data are displayed. Any break-points defined for 
the data (see Section 3.5.2.3 on page 50) and the borders between 
merged data sets (See Section 3.5.2.17 on page 73 and Section 
3.5.2.18 on page 75) are indicated by vertical dashed lines in the 
plots. 


Pop-up menus used to remove mean and trends from individual 
signals (see Section 3.5.2.3 on page 50). The options are: Keep 
mean and trend, Remove mean and trend, Remove mean. 
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Table 3-3. Descriptions of the items of the Pre-processing Tool 


main window, Figure 3-10 (Continued) 


Item 


Description 


Editable fields where break-points can be entered. Break-points split 
the data in parts where each part is treated separately when trends 
and means are removed (see Section 3.5.2.3 on page 50). 


Large plot areas showing the power spectra of a pair of Process 
Variables and Manipulated Outputs/Feed-Forward variables selected 
in the small plot areas in the right part of the window (item 2). Both the 
raw\') and filtered data are displayed. Filter cut-off frequencies are 
indicated by vertical lines in the plot. In addition to indicate the cut-off 
frequencies, the lines can also be used to change the frequencies by 
dragging them with the mouse (see Section 3.5.2.2 on page 49). 


Pop-up menu to select filter type (see Section 3.5.2.2 on page 49). 
The options are: Off (no filter applied), low (low-pass filter), high 
(high-pass filter), pass (band-pass filter), stop (band-stop filter). 


Editable fields where filter cut-off frequencies are entered (see 
Section 3.5.2.2 on page 49). For low-pass and high-pass filters, only 
one of the fields is active. 


10. 


Pop-up menus used to remove mean and trends from all signals (see 
Section 3.5.2.3 on page 50). The options are: Keep mean and trend, 
Remove mean and trend, Remove mean. 


11. 


Button that brings up a dialog where a sampling interval to use for 
resampling is entered (see Section 3.5.2.9 on page 57). 


12. 


Button that brings up a dialog where a new filter order can be entered 
(see Section 3.5.2.3 on page 50). 


(1) Any non-linear transformations and any requested removal of mean and trends are applied also 
to the “raw” data. 
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The menu bar of the Pre-processing Tool main window 


Each heading corresponds to a menu label. 


File 


Edit 


Open 

Open a data file. 

Save 

Save all data to current file. 
Save as 

Save all data to a new file. 
Exit 


Close the window. The user is prompted to save if unchanged data exist. 


History 


Edit the history log for the current file (see Section 3.5.2.14, How to edit 
data descriptions and history logs). 


Description 


Edit the description of the current file (see Section 3.5.2.14, How to edit 
data descriptions and history logs). 


Tools 


Cut 


Remove selection of samples from a data sequences (see Section 3.5.2.6, 
Overview of the Cut Data and Split Data Tools and Section 3.5.2.7, How 
to use the Cut Data Tool). 


Split 


Split a data sequence into two (see Section 3.5.2.6 on page 53 and Section 
3.5.2.8, How to use the Split Data Tool). 
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Resample 

Change the sampling interval (see Section 3.5.2.9, How to resample data). 
Set filter order 

Change the order of the filter (see Section 3.5.2.2, How to filter data). 
Log-file import 


Import a log file and save it as a data-file (see Section 3.5.2.15, Overview 
of the Log-file Import Tool and Section 3.5.2.16, How to use the Log-file 
Import Tool). 


Edit nonlinear functions 


Invoke the Non-linear Transformation Editor (see Section 3.5.2.10, 
Overview of the Non-linear Transformation Editor and Section 3.5.2.11, 
How to use the Non-linear Transformation Editor). 


Remove outliers 


Remove outliers in a data sequence (see Section 3.5.2.12, Overview of the 
Remove Outliers Tool and Section 3.5.2.13, How to use the Remove 
Outliers Tool). 


Merge data-files 


Merge two or more data-files into one (see Section 3.5.2.17, Overview of 
the Merge Data Tool and Section 3.5.2.18, How to use the Merge Data 
Tool). 


Windows 


Switch to another MATLAB figure or the command window. 


3.5.2.2 How to filter data 


The user selects the type of filter from the Filter popup menu (Figure 3-10, item 8). 
Depending on which type of filter that is selected, one or two vertical lines, 
representing the cut-off frequencies, appear in each large frequency plots 

(Figure 3-10, item 7). The cut-off frequencies for the filter can be edited in the Low 
and High edit fields (Figure 3-10, item 9) or changed by dragging the cut-off 
frequency lines with the mouse. 
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The filter order is changed with the Filter order push button (Figure 3-10, item 11) 
or with the Tools > Set filter order menu, (see Section, The menu bar of the Pre- 
processing Tool main window on page 48). In both cases, the dialog of Figure 3-11 
is brought up where the new filter order is entered. The filter order is the order of the 
total filter. For a band-pass filter, that is the sum of the orders of the low-pass and 
high-pass parts of the filter (and equivalently for a band-stop filter). 


MIMO Pre-processing Tool: Filter Order. log971117... FEI ES 


Fill in the new order, and press OK! 
The order will be multiplied by two for 
band pass and band stop filters. 


Filter order: 


Old value: 5 


New value: | 


OK Cancel | 


Figure 3-11. Dialog used to define a new filter order. 


The filter that is used is a butterworth filter. The data is filtered twice, once in both 
directions, which cancels out phase shifts produced by the filter. 


3.5.2.3 How to remove mean and trends from data 
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Mean values or trend values can be removed from individual signals or from all 
signals at once. The Global de-trend popup menu (Figure 3-10, item 10) operates 
on all signals, and the popup menus above the large time series plots (Figure 3-10, 
item 5) operates on the signal currently in the plot. The user has three choices in 
these menus: Keep mean and trend, Remove mean and trend and Remove mean. If 
Remove mean is chosen for a signal in a time series plot and break-points are 
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defined for the signal, the user will be asked if the break-points should be used or 
not. 


When the operation is performed on single signals, the user can define so called 
break-points that will force de-mean and de-trend to treat data between the break- 
points as separate data sequences. 


Break-points are entered for one signal at a time into the Break-points edit field 
above the time series plots Figure 3-10, item 6). The break-point values are entered, 
separated by spaces or commas, between the brackets in the break-points edit fields. 
The unit of the values is the same that is used for the x-axis in the time series plot. 
There can not be more than 10 break-points entered for one signal. The break-points 
are presented as vertical dashed lines in the time series plot. 


3.5.2.4 Overview of the Signal Close-up function 


A close-up view of one of the signals, either as a time-series or a power spectrum, is 
obtained when the user clicks in one of the large time series or frequency plots 
(Figure 3-10 items 4 and 7). The close-up is presented in a window like that in 
Figure 3-12. The menu bar of the window is described in Section, The menu bar of 
the Signal Close-up window on page 53. 
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MIMO Pre-processing Tool - 1og971117_1510_2250.mxd - Close Up Display 


File Tools Windows 


65101262.GIV. 


Figure 3-12. Signal Close-up window. 
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The menu bar of the Signal Close-up window 
Each heading corresponds to a menu label. 
File 
Exit 
Close the window. 
Tools 
Remove outliers 


Open the remove outliers window with current signal (see Section 
3.5.2.12, Overview of the Remove Outliers Tool and Section 3.5.2.13, 
How to use the Remove Outliers Tool). 


Windows 


Switch to another MATLAB figure or the command window. 


3.5.2.5 How to use the Signal Close-up function 


In a Signal Close-up window, the ordinary zoom functionality of MATLAB is 
activated which means that the user can zoom in on certain parts of the plots by 
dragging with the mouse over it (see MATLAB documentation for details). It is also 
possible to start the Remove Outliers Tool (see Section 3.5.2.12, Overview of the 
Remove Outliers Tool and Section 3.5.2.13, How to use the Remove Outliers Tool) 
from the menu of the window. The tool is then initialized with the data in the 
window. 


3.5.2.6 Overview of the Cut Data and Split Data Tools 


The Cut Data Tool is used to remove data in selected time intervals simultaneously 
from all signals in a data sequence. The Split Data Tool is used to split a data 
sequence into two new sequences. 


The tools are started from the Tools > Cut and Tools > Split menus in the Pre- 
processing Tool main window, see Section, The menu bar of the Pre-processing 
Tool main window on page 48. 
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When any of the tools are started, one of the windows shown in Figures 3-13 and 
3-14 is brought up. The items of the windows are described in Table 3-4. 


MIMO Pre-processing Tool: Cut Data from File log971117_1510_2250.mxd. [Joi] 
File Plot Windows 


Change mouse mode to ‘Cut’, select samples and press OK Current mouse mode: Cut Change 5 


‘min 


Figure 3-13. Cut Data window. 
The items of the window are described in Table 3-4. 
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MIMO Pre-processing Tool: Split File log971117_1510_2250.mxd into 2 New Data Sets. [Jol] 
File Plot Windows 


Filln the Sample time, and press OK. 


Current mouse mode: Zoom 


min 
Spit et Sample time: 


—— 


OK Cancel 


7 6 8 


Figure 3-14. Split Data window. 
The items of the window are described in Table 3-4. 
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Table 3-4. The items of the Cut Data and Split Data windows, 
Figures 3-13 and 3-14. 


Item Description 
1. Plot of selected Manipulated Output or Feed-Forward variable. 
2. Pop-up menu to select a Manipulated Output or Feed-Forward 
variable to view. 
3. Plot of selected Process Variable. 
4. Pop-up menu to select a Process Variable to view. 
5. Only in Cut Data windows. Button to change mouse mode. The 


mouse can either be used to zoom in on specific parts of the viewed 
signals, or to select a time-interval in which data should be removed. 
The button toggles the mode. In zoom mode, the zooming is 
synchronized in the two plot areas. In cut mode, a selected time- 
interval will make the data in that interval being removed from all the 
signals in the file, not only the two ones viewed. 


6. Only in Split Data windows. Editable field where the time (counted 
from the beginning of the data) where to split the data is entered. If 
the Return key is pressed, or the mouse is clicked outside the field, a 
solid vertical line is drawn in the two plot areas at the time instant 
where the split will take place. 


7. OK button. When pressed the requested action (selected data 
removed, or the file being splitted) is performed and a file dialog is 
brought up where the user selects the file to put the result in (if the 
action is Split, two dialogs are presented in sequence). The Cut/Split 
Data window is then closed and control is returned to the main 
window. If the action was Split, the window will contain the data of the 
original file. If the action was Cut, the window will contain the data of 
the resulted file, with vertical dashed lines indicating positions where 
data has been removed"), 


8. Cancel button. The Cut/Split Data window is closed without any 
action taking place. 


(1) The indication is the same as for files created by the merging of several data sequences (see 
Section 3.5.2.17 and Section 3.5.2.18 ) and this is how the file actually will be treated. 
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3.5.2.7 How to use the Cut Data Tool 


The default signals in the window are the same as in the time series plots in the Pre- 
processing Tool main window (Figure 3-10, item 4). The user can change the 
signals to be viewed with the popup menus above the plots, Manipulated 
Outputs/Feed-Forward variables in the upper plot and Process Variables in the lower 
plot. The default mouse mode is zoom, which allows the user to zoom in on a 
specific area of the time series plot. To select data to cut, the mouse mode is 
changed to Cut, with the Change push button (Figure 3-13, item 5). When cut mode 
is selected, it is possible to select bad parts, in the data sequences, for removal. The 
selected parts will be marked with Xs. When OK is pushed (Figure 3-13, item 7), 
the selected samples will be removed from all signals and a file dialog is presented 
where the user selects a file to put the result in. The Cut Data window is closed and 
the Pre-processing Tool main window is updated with the data of the resulted file. 
The area where the cut has been made will be marked with a vertical dashed line in 
the time series plots. The Cancel button (Figure 3-13, item 8) closes the Cut Data 
window without performing any action. 


3.5.2.8 How to use the Split Data Tool 


The user enters the split time in the editable field (Figure 3-14, item 6). When the 
split time is entered (after the Return key is hit, or the mouse is clicked outside the 
field), a solid vertical line in each plot shows where the time series will be split. 
When OK is pushed (Figure 3-14, item 7), the time series to the left of the split line 
will be saved into a new file and the time series to the right of the split line will be 
saved to another new file. Two file dialogs are presented in sequence where the user 
selects the files to save the results in. The Cancel button (Figure 3-14, item 8) 
closes the Split Data window without performing any action. 


3.5.2.9 How to resample data 


By pushing the Sampling interval push-button in the Pre-processing Tool main 
window (Figure 3-10, item 11) or selecting the Tools > Resample menu of the main 
window (see Section The menu bar of the Pre-processing Tool main window on 
page 48) the dialog in Figure 3-15 is presented. 


3BSE 011 572R101 57 


3dMPC User’s Guide 
Chapter 3 Configuration/Application Building 


MIMO Pre-processing Tool: Resample. log971117_... i=] E3| 


Figure 3-15. Dialog where a sampling interval for resampling is entered. 


The user enters the new sampling interval in the New value edit field. When the OK 
button is pushed, the dialog is closed, the resampling is performed and a file dialog 
is brought up where the user selects a file to save the result in. The Pre-processing 
Tool main window is updated with the resampled data. The Cancel button closes 
the dialog with no resampling performed. 
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3.5.2.10 Overview of the Non-linear Transformation Editor 


The Non-linear Transformation Editor is started from the Tools > Edit nonlinear 
functions menu in the Pre-processing Tool main window, see Section , Overview of 
the Pre-processing Tool. 


The editor is used to define and edit transformations that linearize signals contained 
in data sequences. The transformations are used both in the Pre-processing Tool to 
produce linearized data sequences and in the On-line Controller. 


The class of transformations allowed is that of invertible piece-wise linear 
functions. The transformations are defined by two vectors specifying the break- 
points of the piece-wise linear functions, given as x-y coordinates. For the functions 
to be invertible, they have to be monotonous (either increasing or decreasing). 


When the editor starts, the main window of Figure 3-16 is shown. The different 
items in the window are described in Table 3-5 and the menus are described in 
Section, The menu bar of the Non-linear Transformation Editor main window on 
page 63. 


During a session with the editor, all transformations are held in the internal memory 
of the editor and are not transferred to the Pre-processing Tool and applied to the 
data until the editor is closed. 
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Figure 3-16. Non-linear Transformation Editor main window. 
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Table 3-5. Items of the Non-linear Transformation Editor 
main window, Figure 3-16. 


Item 


Description 


The menu bar. Described in Section, The menu bar of the Non-linear 
Transformation Editor main window on page 63. 


Edit area where the transformations are displayed and also edited 
graphically with the mouse. The break-points defining the 
transformation are displayed by bullets. 


Push button that reverts the transformation currently edited in the edit 
area (item 2) to that saved internally in the editor. 


Push button that clears the edit area (item 2). 


Toggle button. When pressed, new break-points are added to a 
transformation by clicking with the left mouse button in the edit area 
(item 2). If the new break-point would make the function non- 
monotonous it is not accepted. By clicking with the right mouse, the 
break-point closest to where clicked is removed. 


Push button that when pushed brings up a dialog where the break- 
points can be entered as two MATLAB row vectors"). One for the x- 
values and one for the y-values. 


Push button that when pushed brings up a dialog where any valid 

MATLAB expression returning a real numeric matrix with two rows can 
be entered. The matrix will define the break-points, with the x-values 
taken from the first row and the y-values taken from the second row. 


Push button that when pushed brings up the dialog of Figure 3-17 
where a standard non-linear function can be selected. 


(1) 
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Square root ha 


Figure 3-17. Standard non-linear function dialog. 
The dialog is brought up from the Non-linear Transformation Editor main window, 
Figure 3-16. 


Table 3-6. Items of the standard non-linear function dialog, 


Figure 3-17 
Item Description 
1. Pop-up menu from where the non-linear function is chosen. 
2. String displaying the parameterization of the selected non-linear 
function. 
3. Fields where values for the parameters of the chosen non-linearity 
are entered. 
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Table 3-6. Items of the standard non-linear function dialog, 
Figure 3-17 (Continued) 


Item Description 
4. Field where the number of break-points to use in the piece-wise linear 
function approximation is entered. 
5. Ok button that when pushed computes and displays the piece-wise 
linear function approximation and closes the dialog. 
6. Cancel button closing the dialog. 


The menu bar of the Non-linear Transformation Editor main window 


Each heading corresponds to a menu label. 
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File 


Edit 


Save 


Saves the transformation currently edited in the edit area (Figure 3-16, 
item 2) in the internal memory of the editor. 


Close 


Returns the changed non-linear transformations to the Pre-processing Tool 
main window and closes the editor. 


The sub-menus contain together all the signals available. The menus are used 
to select which transformation to display and edit in the edit area (Figure 3-16, 
item 2). 


PV 

Select a Process Variable. 
OUT 

Select a Manipulated Output. 
FF 


Select a Feed-Forward variable. 
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Change y-range 


Change the range for the y-axis. 


3.5.2.11 How to use the Non-linear Transformation Editor 
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When the editor is started, all defined transformations are loaded into the internal 
memory of the editor and the transformation for the first Process Variable (PV 1) 
will be displayed in the edit area of the main window of the editor (Figure 3-16, 
item 2). If no transformation is defined for the signal, the edit area will be empty. 


Another transformation to edit can be selected with the Edit > PV, Edit > OUT and 
Edit > FF menus (see Section, The menu bar of the Non-linear Transformation 
Editor main window on page 63). The default x and y-ranges are the ranges 
specified in the configuration file (see Section 3.5.1, Configuration Tool). If there 
are no ranges specified in the configuration file, 0 to 10 are used. If the signal to edit 
contains a non-linear definition, it is expanded to cover the whole range before 
being displayed. 


All the edited or created non-linear transformations are stored locally in the editor 
until it is closed. 


There are four different ways to edit the non-linear transformations: 


To enter or delete points in the non-linear transformation, push the Add/delete 
points button (Figure 3-16, item 5) and enter points with the left mouse button 
and delete the point that is closest to the mouse pointer with the right mouse 
button. To stop the Add/delete points mode push the button a second time. It is 
not possible to add points that make the transformation non-monotonous. 


The Edit curve button (Figure 3-16, item 6) shows the values for the points in 
the non-linear transformation. The values are presented as two vectors with the 
same length, one for the y-values and one for the x-values, which can be edited. 


The MATLAB expression button (Figure 3-16, item 7) makes it possible to 
specify the transformation in the form of a MATLAB statement. The statement 
should return a real numeric matrix with two rows. The first row will be used 
as the x-coordinates of the break-points defining the transformation, and the 
second row will be used as the y-coordinates. The statement could for instance 
be a function that returns the required matrix, making it possible to create 
libraries of non-linear characteristics of actuators and measurement devices. 
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Extra break-points are added at the ends of the line to make the definition cover 
the whole range. 


e¢ The Standard function button (Figure 3-16, item 8) gives the possibility to 
choose the non-linearity as a standard function. The button brings up the dialog 
of Figure 3-17 from where one of the following non-linear functions can be 
chosen: 


=atbj(x+c) 


= atb(x+c/y 


y 

y 

y = at+blIn(x+0c) 
y = a+blog(x+c) 
y 


= atbe"*® 
with a, b, c and p user supplied parameters. In is the natural logarithm and log 
the base 10 logarithm. 


If any break-points contain complex values or “Not a Numbers”, they are 
removed with a warning message displayed. If the chosen function is non- 
monotonous in the current interval, a dialog is brought up where it is possible 
to choose if the definition should be restricted to be monotonous by keeping 
only the first part of the curve, the last part of the curve, or keeping all the 
points irrespective of the fact that it is non-monotonous. If the third alternative 
is chosen, the Edit curve or Add/delete points functions must then be used to 
make the function monotonous before it is possible to save it. 


An edited transformation can be saved into the internal memory of the editor 
through the File > Save menu. When a new signal is selected through the Edit 
menu, or when the editor is exited, it is checked if the transformation of the current 
signal has been modified, and if so, a dialog is brought up where the saving of the 
transformation can be acknowledged or cancelled. Before a transformation is saved, 
it is checked that it is monotonous, and if not, the saving is prohibited. 


The Revert definition button (Figure 3-16, item 3), replaces the non-linear 
transformation currently edited in the edit area of the main window of the editor 
(Figure 3-16, item 2) with the one saved internally in the editor. The Delete curve 
button (Figure 3-16, item 4) clears the edit area of the editor main window. 
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When the File > Close menu is selected, the signals with changed non-linear 
transformations are updated in the Pre-processing Tool main window and the editor 
is closed. 


3.5.2.12 Overview of the Remove Outliers Tool 


The Remove Outliers Tool is started from the Tools > Remove outliers menu in the 
Pre-processing Tool main window, see Section, The menu bar of the Pre-processing 
Tool main window on page 48, or from the same menu in a Signal Close-up 
window, see Section, The menu bar of the Signal Close-up window on page 53. 


Figure 3-18 shows the Remove Outliers Tool main window which is used to remove 
outliers in one or more time series. 
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MIMO Pre-processing Tool - 1og971117_1510_2250.mxd - Remove Outliers BEE 
File Windows 
Current mouse mode: — Remove outlier 
851Q1262.G1V 
a 20 40 60 80 100 120 140 min 
OK | Cancel 


Figure 3-18. Remove Outliers Tool main window. 
Descriptions of the different items are given in Table 3-7. 


The different items of the window are described in Table 3-7. 
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Table 3-7. Descriptions of the Remove Outliers Tool main window, Figure 3-18. 


Item 


Description 


1. 


Area where the currently treated signal is displayed. 


2. 


Pop-up menu where a Process Variable to treat is selected. When a 
new signal is selected, any Manipulated Output or Feed-Forward 
variable that may have been selected through the pop-up menu, item 
3, is deselected. 


Pop-up menu where a Manipulated Output or Feed-Forward variable 
to treat is selected. When a new signal is selected, any Process 
Variable that may have been selected through the pop-up menu, item 
2, is deselected. 


Button to change mouse mode. The mouse can either be used to 
zoom in on specific parts of the viewed signals, or to select samples 
to be removed. The button toggles the mode. 


OK button. When the button is pushed, all signals treated during the 
current Remove Outlier session are transferred back to the Pre- 
processing Tool main window and the Remove Outlier main window 
is closed. 


Cancel button. When the button is pushed, the Remove Outlier main 
window is closed without influencing the data in the Pre-processing 


Tool main window. 


3.5.2.13 How to use the Remove Outliers Tool 
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If the Remove Outliers Tool is started from a Signal Close-up window (see Section 
3.5.2.4, Overview of the Signal Close-up function and Section 3.5.2.5, How to use 
the Signal Close-up function) the signal displayed in the main window of the tool 
(Figure 3-18, item 1) is the same as that in the Close-up window. If the tool is 
started from Pre-processing Tool main window, no signal is displayed. The user can 
change the signals to be viewed with the popup menus above the plot (Figure 3-18, 
items 2 and 3), Process Variables to the left and Manipulated Outputs/Feed-Forward 
variables to the right. The default mouse mode is zoom, which allows the user to 
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zoom in on a specific area of the time series plot. To select outliers, the mouse mode 
is changed to Remove outliers, with the Change push button (Figure 3-18, item 4). 
When Remove outliers mode is selected, it is possible to select one or more 
samples, in the time sequence. The selected samples will be replaced with a straight 
line between the point to the left and to the right of them (the original appearance of 
the modified part of the curve is shown as a dotted line). The user can remove one or 
more outliers for one signal and then change signal and continue to remove outliers. 
The changes are stored locally in the Remove Outliers main window until either the 
OK button or the Cancel button is pushed (Figure 3-18, items 5 and 6). When OK 
is pushed, the data in the Pre-processing Tool main window is updated. Cancel 
discards the changes for all signals. In both cases, the Remove Outliers main 
window is closed. 


3.5.2.14 How to edit data descriptions and history logs 


Saved together with each data sequence are a description of the data and a history 
log. The description is always entered by a user, either when the logging of an 
identification experiment is started in the 3dMPC On-line controller (see Section 
4.4, Operator Interface) or here in the Pre-processing Tool. The history log is 
updated automatically by the Pre-processing Tool to document the different 
operations performed on the data (see Section 3.5.2.1, Overview of the Pre- 
processing Tool), but it is still possible for the user to manually edit it. 


The editing is performed in edit windows opened from any of the Edit > History or 
Edit > Description menus in the Pre-processing Tool main window (see Section, 
The menu bar of the Pre-processing Tool main window on page 48). The window is 
shown in Figure 3-19. 
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4 MIMO Pre-processing Tool: Edit Description in log971117_1510_2250.mxd._ 


Description of the files merged into log971117_1510_2250.mxd 


log971117_1510.mxd: 


log971117_2250.mxd: 


Figure 3-19. Window used to edit data descriptions. 
An equivalent window is used to edit the history log of the data. 


The existing description or history log is shown in the window. To scroll the 
window use the Page Up, Page Down, the arrow keys, Ctrl-End and Ctrl-Home 
on the keyboard. The user can edit the existing history inside the edit field. The OK 
button updates the history in the Pre-processing Tool main window and the Cancel 
button discards the changes. 


3.5.2.15 Overview of the Log-file Import Tool 


The Log-file Import Tool is used to convert log files created by the 3d MPC On-line 
controller (see Section 4.4, Operator Interface) to data sequence files of the format 
used by the Engineering Tools. 


The tool is started from the Tools > Log-file import menu of the Pre-processing 
Tool main window (see Section, The menu bar of the Pre-processing Tool main 
window on page 48). The main window of the tool is shown in Figure 3-20 with 
descriptions in Table 3-8. The menu of the window is described in Section, The 
menu bar of the Log-file Import Tool main window on page 72. 
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Log-file Import 


Figure 3-20. Log-file Import Tool main window. 
Descriptions in Table 3-8. 
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Table 3-8. 


Descriptions of the Log-file Import Tool main window, Figure 3-20. 


Item 


Description 


1. 


The menu bar. Described in Section, The menu bar of the Log-file 
Import Tool main window on page 72. 


Button that when pushed brings up a file dialog where the file to 
convert is selected. 


Editable field where the sampling interval is entered. Each saved data 
sequence is assumed to contain equidistantly sampled data and the 
used sampling interval is saved together with the data. The field is 
initialized with the median value as measured from the log-file (see 
item 4). 


Statistics regarding the sampling interval used when the log-file was 
created. Each value in the log-file is time stamped and from these 
time stamps the median value and the bounds for the low and high 
quartiles of the obtained sampling intervals are calculated and 
displayed. 


The number of samples in the log-file. 


Window displaying some information about the log-file. 


The menu bar of the Log-file Import Tool main window 


Each heading corresponds to a menu label. 


File 


Save as... 


Save the converted log in a selected file. 


Exit 


Close the window. 


Windows 


Switch to another MATLAB figure or the command window. 
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3.5.2.16 How to use the Log-file Import Tool 


The user selects the log-file to convert in a file dialog that is presented when the 
Log-file push button (Figure 3-20, item 2) is pushed. The sampling interval is 
entered in the Select sampling interval edit field (Figure 3-20, item 3). The entered 
sampling interval is the defined time between the samples in the data-file. The 
measured sampling interval in the log-file is presented as the upper bound of the low 
quartile, the median and the lower bound of the high quartile. The median value is 
the default value for the Select sampling interval edit field. The number of samples 
is also presented (Figure 3-20, item 5) and the Info text field (Figure 3-20, item 6) 
contains the header of the log-file. 


The File > Save as... menu (see Section, The menu bar of the Log-file Import Tool 
main window on page 72) presents a file dialog in which the user can enter the data- 
file name, where the data will be saved. 


3.5.2.17 Overview of the Merge Data Tool 


The Merge Data Tool is used to merge data sequences collected at different 
occasions and possibly covering different parts of the process. Up to ten files can be 
merged together at a time into one new file. The files that are merged can 
themselves be the result of a merge, making it possible to merge an unlimited 
number of files together. 


The tool is started from the Tools > Merge data files menu of the Pre-processing 
Tool main window (see Section, The menu bar of the Pre-processing Tool main 
window on page 48). 


The main window of the tool is shown in Figure 3-21 with the menus described in 
Section, The menu bar of the Merge Data Tool main window on page 75. 
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MIMO Pre-processing Tool - Merge Data-files 


Figure 3-21. Merge Data Tool main window. 
Descriptions in Table 3-9. 


Table 3-9. Descriptions of The Merge Data Tool main window, Figure 3-21. 


Item Description 
1 The menu bar. Described in Section, The menu bar of the Merge 
Data Tool main window on page 75. 
2. Buttons that when pushed brings up file dialogs where the files to 


merge are selected. 
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The menu bar of the Merge Data Tool main window 
Each heading corresponds to a menu label. 
File 
Save as 
Merge the selected files into one and save it in a new file. 
Exit 
Close the window (no merge). 
Windows 


Switch to another MATLAB figure or the command window. 


3.5.2.18 How to use the Merge Data Tool 


When any of the File # push buttons (Figure 3-21, item 2) is pushed, a file dialog is 
displayed in which the user can enter one of the files to be merged. Ten files is the 
maximum number of files that can be merged, and File 1 will be the first in the new 
data file. The files should have the same sampling interval, otherwise it is not 
possible to merge them. The File > Save as menu (see Section, The menu bar of 
the Merge Data Tool main window on page 75) will display a file dialog in which 
the user can specify a file where the merged data will be saved. When a file that are 
the result of a merge is displayed in the Pre-processing Tool main window, the 
borders between the original files are marked with vertical dashed lines in the time 
series plots (Figure 3-10 on page 45, item 4). 
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3.5.3 Modeling Tools 


This section describes the user interface of the Modeling Tools of 3dMPC. 

In addition to the material here, Section B.5, Identification in Chapter B, 
Configuration Tutorial contains a worked identification example and treats some of 
the most common considerations that has to be taken when creating a model from 
logged data. 


The purpose of the modeling tools is to create the models needed to design the MPC 
controller (and also used on-line by the 3dMPC On-line controller). The models are 
either identified from process data or given as simple user defined models. 


In addition to the model creation tools there are also tools to examine the models in 
different respects and there is a tool to merge models created at different occasions 
and covering different parts (that is, containing different sets of in-, and outputs) of 
a process. 


The contents of the section are summarized below 
Identification page 79 


Describes the Identification Tool that is used to create state-space models from 
collected process data. 


Overview of the Identification Tool page 79 
Gives an overview of the functionality of the Identification Tool 
How to use the Identification Tool page 87 
Gives detailed descriptions on the use of the Identification Tool 
Basic file handling in the Identification Tool page 87 
Describes the basic file handling functionality of the Identification Tool 
Model scan page 89 


Describes how to perform a scan over the parameters n, beta and gamma 
that are inputs to the Subspace identification algorithm. 
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Batch identification page 95 


Describes how to create a number of models in one run with the Subspace 
identification algorithm, with the models having different inputs and 
outputs and where different values for n, beta and gamma (input 
parameters to the algorithm) are used. 


Prediction error identification page 96 


Describes how to use the Prediction Error identification algorithm to 
create a final model. 


Model Evaluation page 98 


Describes the Model Evaluation Tool used to examine the quality of models created 
either by the Identification Tool or as User Defined Models. 
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Overview of the Model Evaluation Tool page 98 
Gives an overview of the functionality of the Model Evaluation Tool 
How to use the Model Evaluation Tool page 106 
Gives detailed descriptions on the use of the Model Evaluation Tool 

Basic file handling in the Model Evaluation Tool page 106 


Describes the basic file handling functionality of the Model Evaluation 
Tool 


Plots page 107 
Describes how the different plots used to evaluate the models, are created. 
Simulated/Predicted Outputs page 116 


Describes how the measured outputs of a given evaluation data set are 
plotted together with either the simulated or the predicted model outputs. 


Prediction Error Variances page 117 


Describes how the prediction error variances as a function of prediction 
length are plotted for given models and a given evaluation data set. 
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Step responses page 117 
Describes how the step responses of all transfer functions of given models 
are plotted. 

Frequency Responses page 117 


Describes how gain and phase curves for all transfer functions of given 
models are plotted. 


Impulse Responses page 118 


Describes how the impulse responses of all transfer functions of given 
models are plotted. 


Noise Spectra page 118 
Describes how the noise spectrum for each output of given models are 
plotted. 

Poles and Zeros page 118 


Describes how the poles and zeros for all transfer functions of given 
models are plotted. 


Model Merging page 119 


Describes the Model Merging Tool used to merge together two models covering 
different parts of a process. 


Overview of the Model Merging Tool page 119 
Gives an overview of the functionality of the Model Merging Tool 
How to use the Model Merging Tool page 123 
Gives detailed descriptions on the use of the Model Merging Tool 
Basic file handling in the Model Merging Tool page 123 
Describes the basic file handling functionality of the Model Merging Tool 
Merging the models page 124 


Describes the actual commands used to merge two models. 
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User Defined Models page 125 


Describes the User Defined Models Tool used to define and edit simple models 
directly (without identifying from data). 


Overview of the User Defined Models Tool page 125 
Gives an overview of the functionality of the User Defined Models Tool. 
How to use the User Defined Models Tool page 129 


Gives detailed descriptions on the use of the User Defined Models Tool. 


3.5.3.1 Identification 


Overview of the Identification Tool 


The Identification Tool is started from the Tools > Modeling > Identification menu 
of the Engineering Tools control panel (see Section 3.3, Engineering Tools 
Overview). The main window of the tool is shown in Figure 3-22 on page 80 with 
descriptions of the different items given in Table 3-10 on page 81 
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(fd Add beta 


MIMO Identification 


Figure 3-22. Identification Tool main window. 
The upper view is the default one, with the batch identification functionality disabled, and the lower 
view is the appearance with the batch identification enabled. The items of the window are described 
in Table 3-10 on page 8&1. 
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Table 3-10. Descriptions of the items of the Identification Tool 


main window, Figure 3-22 


Item 


Description 


The menu bar. Described in Section, The menu bar of the 
Identification Tool main window on page 84. 


Space where the name of the loaded Identification data is displayed. 


Space where the name of the loaded Evaluation data is displayed. 


Column of check-boxes where the Process Variables to include in the 
identification are selected. The names of the variables (as given by 
the configuration, see Section 3.5.1 on page 34) are displayed at 
each box. There are equivalent columns for the Manipulated Outputs 
and Feed-Forward variables. 


Column of fields where known delays in the Process Variables are 
entered as number of samples. There are equivalent columns for the 
Manipulated Outputs and Feed-Forward variables. 


Column of fields where scaling factors to apply to the Process 
Variables are entered. Scaling factors can improve the numerical 
properties of the algorithms when the magnitudes of the signals differ 
very much. There are equivalent columns for the Manipulated 
Outputs and Feed-Forward variables. 


Column with statistical estimates of the standard deviations of the 
Process Variables. The estimates are based on the loaded 
Identification data. There are equivalent columns for the Manipulated 
Outputs and Feed-Forward variables. 


Column with the range widths (rangeHi-rangeLo) for the Process 
Variables. There are equivalent columns for the Manipulated Outputs 
and Feed-Forward variables. 


Fields where intervals of n, beta and gamma, input parameters to 
the Subspace identification algorithm, are entered. The tool can 
perform scans over all models corresponding to combinations of the 
parameters in the intervals (see Section, Model scan on page 89). 
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Table 3-10. Descriptions of the items of the Identification Tool 


main window, Figure 3-22 (Continued) 


Item 


Description 


Button that when pushed inserts standard intervals for the n, beta 
and gamma ranges of item 9. The ranges are based on the number of 
signals to be included in the models, so those must have been 
selected (see item 4) before pushing the button. 


Check-box controlling how the input gamma to the subspace algorithm 
is chosen during a model scan (see Section, Model scan on page 
89). When the box is checked, each gamma _ value is obtained as the 
sum of one value in the beta interval and one in the gamma interval 
defined by the fields of item 9. This is indicated by the string “+ beta” 
being displayed to the right of the gamma interval fields. When not 
checked, the gamma and beta intervals are independent. 


12. 


Toggle button. When the button is in its pressed state, the main 
window of the Identification Tool is expanded to include a part where 
a batch identification can be defined (see Section, Batch 
identification on page 95). The File and Subspace identification 
menus will also get a couple of new items concerning the batch 
identification (see Section, The menu bar of the Identification Tool 
main window on page 84) 


3BSE 011 572R101 


3dMPC User's Guide 
Section 3.5.3 Modeling Tools 


Table 3-10. Descriptions of the items of the Identification Tool 


main window, Figure 3-22 (Continued) 


ltem 


Description 


13. 


Radio-button used to select one out of 15 rows of editable fields and 
push-buttons. Only five rows are visible at a time, which ones 
controlled by the pop-up menu, item 21. Each row logically 
corresponds to a memory slot where a connection structure can be 
saved. By a connection structure we mean all the inputs to the 
subspace identification algorithm: the parameters n, beta and 
gamma, the signals to include in the identified model and the delays 
and scalings associated with the signals. A connection structure is 
loaded into the slot by pushing the button, item 20. The loaded values 
for the parameters n, beta and gamma are those entered in the 
editable fields, item 14, of the selected row and the loaded signals, 
delays and scalings are those defined by the items 4, 5 and 6. When 
a slot contains a loaded connection structure, the string Structure 
Number <i> loaded is displayed directly to the right of the radio- 
button and when a radio-button corresponding to such a slot is 
pushed, the items 4, 5 and 6 are updated to display the values of the 
loaded structure. The tool can in one batch create all the models 
corresponding to the loaded connection structures (see Section, 
Batch identification on page 95). 


Editable field where the parameters n, beta and gamma fora 
connection structure (see item 13) are entered. 


Editable field where to enter the name of a file to save a model 
corresponding to a connection structure (see item 13) in. 


Button that when pushed brings up an edit window where to enter a 
description to save together with a model corresponding to a 
connection structure (see item 13). 


Check-boxes used to select rows from which loaded connection 
structures (see item 13) are cleared. 


18. 


Push-button that when pushed clears checked (through item 17) rows 
from loaded connection structures (see item 13). 
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Table 3-10. Descriptions of the items of the Identification Tool 
main window, Figure 3-22 (Continued) 


Item Description 

19. Push-button that when pushed clears the row currently selected by a 
radio-button (item 13) from any loaded connection structure. 

20. Push-button that when pushed loads a connection structure into the 
row currently selected by a radio-button (item 13). 

21. Pop-up menu used to select which five connection structure rows out 
of 15 to show (see item 13). 


The menu bar of the Identification Tool main window 


Each heading corresponds to a menu label. 


File 


Load identification data set... 
Load a data file to use in the identification. 
Load evaluation data set... 


Load a data file to use when ranking models calculated in the scanning 
over n, beta and gamma, see Section, Model scan on page 89. 


Load set of connection structures... 


This menu is shown only when the batch identification is enabled (the 
Batch identification toggle button in the main window is pressed, see 
Figure 3-22 on page 80). It loads a file containing a set of defined 
connection structures, see the Section, Batch identification on page 95. 


Save a set of connection structures... 


This menu is shown only when the batch identification is enabled (the 
Batch identification toggle button in the main window is pressed, see 
Figure 3-22 on page 80). It saves a set of defined connection structures to 
a file, see Section, Batch identification on page 95. 
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Delete 
delete data... 
Delete selected data files from disk. 
delete model... 
Delete selected model files from disk. 
Reload configuration file 


Reload the configuration file (see Section 3.5.1 on page 34) from disk. 
This is needed only if the file has changed and the changes are to be used 
during the current session. 


Evaluation of the models 


Bring up the user interface for the Model Evaluation Tool, see Section 
3.5.3.2 on page 98. 


Exit 
Exit the Identification Tool. 
Prediction error identification 
Start with models 
<List of available model files> 


The sub-menu has one item for each available model file. Selecting 
one of the items starts a prediction error identification initialized with 
the selected model, see Section, Prediction error identification on 
page 96. 


Subspace identification 
Scan n, beta and gamma 


Start a scan over given intervals of the parameters n, beta and gamma (see 
Section, Model scan on page 89). 
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Present last scan result 


Brings up the display showing the ranking of the models obtained in the 
last scan over the parameters n, beta and gamma, see Section, Model 
scan on page 89. 


Batch identification 


This menu is shown only when the batch identification is enabled (the 
Batch identification toggle button in the main window is pressed, see 
Figure 3-22 on page 80). 


all structures 


Make a batch subspace identification for all connection structures 
currently defined, see the Section, Batch identification on page 95. 


current structure 


Make a subspace identification for the currently selected connection 
structure, see Section, Batch identification on page 95. 


Information 
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Information about the data sets 


<List of available data files> 


The sub-menu has one item for each available data file that can be 
used for identification or model evaluation. Selecting one of the items 
brings up a window showing some information about the data. The 
shown information is: the sampling interval for the data, the names of 
all signals in the file, any non-linear transformations applied to the 
signals, the history of the data (a list of all the processing steps the 
data has gone through) and any textual description saved together 
with the data (descriptions can be entered when a log is started in the 
3dMPC On-line controller, see Section 4.4 on page 157, and in the 
Pre-processing tool, see Section 3.5.2 on page 42). 


If an identification data set is loaded, the first two items of the sub- 
menu are labeled current id-data set and current ev-data set and 
selecting any of them brings up the above information for the loaded 
identification data set or evaluation data set (when an identification 
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data set is loaded, but no evaluation data set, the identification data is 
used also as evaluation data). 


Information about the models 
<List of available model files> 


The sub-menu has one item for each available model file. Selecting 
one of the items brings up a window showing some information 
about the model. The shown information is: the sampling interval for 
which the model is valid, the names of all signals in the file, any non- 
linear transformations applied to the signals, any textual description 
saved together with the model and the history of the model (a list 
describing how the model was created). 


Windows 


The menu has items for each existing MATLAB figure window, and one item for 
the MATLAB command window. Selecting one of the items just brings the 
corresponding window to the front. 


How to use the Identification Tool 


Basic file handling in the Identification Tool 
Identification/evaluation data sets 


An identification data set is loaded from the File > Load identification data set... 
menu of the main window. A file list containing all available data files is brought up 
from where the file to load is selected. The name of the loaded file is displayed in 
the upper left corner of the window (see Figure 3-22 on page 80, item 2). 


Until an identification data set is loaded, all items below the menu bar as well as the 
menu Subspace identification in the main window are disabled (see Section, The 
menu bar of the Identification Tool main window on page 84). This applies also to 
the items Load evaluation data set..., Load set of connection structures... and 
Save set of connection structures... in the File menu (the last two menus are 
invisible if the Batch identification toggle-button, item 12 in Figure 3-22 on page 
80, is not pressed). 
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The identification data set will be used for all subsequent subspace identifications. 
When it has been loaded, all menus (except the Options menu that will be disabled 
as long as no scan over the parameters n, beta and gamma has been done, see 
Section, Model scan on page 89) will be enabled. All check-boxes, edit fields and 
pop-up menus in the display will be enabled, except check-boxes and edit fields 
corresponding to signals missing (as compared to the configuration, see Section 
3.5.1 on page 34) in the data file. 


When loaded for the first time, the identification data set will also be set as the 
evaluation data set (with the name displayed at the top of the main window, item 3 
of Figure 3-22 on page 80). The evaluation data set is used when ranking models 
calculated in the scanning over n, beta and gamma (see Section, Model scan on 
page 89). A separate data set can loaded from the File > Load evaluation data set... 
menu. As for the identification data sets, any of the data files contained in the 
current application can be loaded. 


Before a data set is loaded the Information > Information about the data sets 
menu can be used to bring up windows containing information on selected data sets. 
There is also an Information > Information about the models menu for bringing 
up the corresponding information on created models. 


Any of the data or model files contained in the current application can be deleted 
from disk using the File > delete data... or File > delete model... menus. A list of 
the available files is presented in a pop-up menu where one file at a time can be 
selected and deleted. 


Connection structure sets 


It is possible to define up to 15 different sets of connection structures (all inputs 
needed by the subspace identification algorithm) and to identify all the 
corresponding models in one batch (see the Section, Batch identification on page 
95). 


All connection structures currently loaded in memory can be saved to a file with a 
user-supplied name. This is done through the File > Save set of connection 
structures... menu. 


A saved set of connection structures can be reloaded at a later occasion using the 
File > Load set of connection structures... menu. A file dialog containing all 
connection structure sets available in the current application is brought up from 
which the set to load is selected. 
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Model scan 


The subspace identification algorithm needs three user-supplied parameters, n, 
beta and gamma. The parameter n is the number of states in the identified state- 
space model and beta and gamma are used internally in the algorithm to create a 
preliminary predictor model. The predictor uses beta old process input and output 
vectors and gamma future process input vectors to predict gamma future process 
output vectors. The parameter Gamma also puts a limit on the number of states n. 
The limit is that n must be less or equal to gamma - 1 (see also Section A.1.1 on page 
193). 


The values to use for theses parameters can be difficult to decide in advance. To 
help in the decision, there is a function for scanning over intervals of the (integer) 
parameters. The intervals are defined by entering values in the edit fields just below 
the signal selection check boxes of the main display (item 9 in Figure 3-22 on page 
80). 


By pushing the Standard intervals button (item 10 in Figure 3-22 on page 80), 
intervals based on the number of signals to include in the identified models are filled 
in. With many signals included, the upper limits of the intervals can be very large, 
and to decrease the time of the scan, it might be necessary to manually fill in smaller 
values. See also the discussion in the Configuration Tutorial, Section B.5.2 on page 
213. 


All models created during the scan uses the same input and output signals. Which 
signals to include is determined by checking the corresponding check-boxes in the 
main window (item 4 in Figure 3-22 on page 80). 


For each signal there is an edit field in which known delays in the signal is entered 
as number of samples (items 5 in Figure 3-22 on page 80). If a delay is known, the 
model quality will be improved by leaving the algorithm from the task of estimating 
it, but if a too large value is entered, it is impossible for the algorithm to correct for 
the failure. This means that in case of uncertainty, it is better to underestimate than 
overestimate the delays. 


If the magnitudes of the signals differ very much, the numerical reliability of the 
identification is improved if the signals are scaled. A multiplicative scaling factor 
can be entered for each signal (item 6 in Figure 3-22 on page 80). To support the 
selection of scaling factors, the standard deviation of the signals (calculated based 
on the identification data set) and the range of the signals (as defined by the 
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configuration file, see Section 3.5.1 on page 34) are displayed in the window (items 
7 and 8 in Figure 3-22 on page 80). 


The scanning is started from the Subspace identification > Scan n, beta and 
gamma menu. Models are then calculated for all combinations of n, beta and 
gamma. 


If either no PV or neither any OUT nor any FF signal is selected, a warning message 
is displayed informing that at least one PV and one OUT or FF signal must be 
selected, and the scanning is terminated. 


When the calculations started, the window shown in Figure 3-23 on page 90 is 
brought up. It shows how many of the models that has been calculated so far and the 
maximum number of models to calculate. The window also contains a button 
permitting the calculations to be terminated in advance. 


Progress Info BEE 
Max. number of models to calculate Done so far 
153 7 

Terminate the calculations | 


Figure 3-23. Window showing the progress of a model scan. 


The number of models can be very large so it is not feasible to save all the models to 
disk. Instead, for each calculated model a measure of the predicting quality of the 
model is computed. The measure consists of a weighted sum of the prediction error 
variances produced by the model using the evaluation data. The sum is taken both 
over the different Process Variables and over different prediction lengths. 


The result of a scan is presented in a window as that shown in Figure 3-24 on page 
91 with descriptions in Table 3-11 on page 92. 
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Figure 3-24. Window presenting the result of a model scan. 
The items are described in Table 3-11. 
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Table 3-11. Descriptions of the items of a model scan result window, 


Figure 3-24 


Item 


Description 


The File menu. It contains only a Close command that closes the 
window. A window presenting the latest scan can always be 
recreated through the Subspace identification > Present last scan 
result menu of the main window of the Identification tool. 


The Show menu. The menu controls what is plotted in the plot area 
(item 7). The menu items are: all models, 20 best, Set... . If the first 
item is selected (the default), the plot shows the prediction error 
variance measure as a function of all models included in the scan. If 
the second item is selected, only the 20 models having the lowest 
value for the measure are included in the plot and the last item brings 
up a dialog where the user can enter how many models to include in 
the plot. 


The Identification menu. It contains only the item: Start subspace 
identification with saving. Selecting it, makes the models checked 
in the list in the right part of the window (item 8) being recreated and 
saved to disk. 
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Table 3-11. Descriptions of the items of a model scan result window, 
Figure 3-24 (Continued) 


Item Description 


4. The Evaluate selected models menu. The menu contains the items 
Model Output, Prediction Error Variances, Step Response, 
Frequency Response, Impulse Response, Noise Spectrum and 
Poles/Zeros. The items are divided into two groups, the three first 
and the rest. Only the first group is visible when the window is first 
brought up. In the Options menu (see below) there is an item More 
evaluation commands that can be checked to make all the items 
available. Selecting any of the items makes all the models checked in 
the check-boxes (item 8) being plotted in a new window using the 
corresponding function of the Model Evaluation Tools (see Section 
3.5.3.2, Model Evaluation). If all the models checked already exist on 
disk, the corresponding evaluation tool is called directly. If none of the 
models exist, they are first created by running a subspace 
identification. If some of the models exists, a dialog is brought up 
where the user can choose if only the existing models should be 
plotted, or if all the models first should be identified (destroying the 
existing ones). 


5. The Options menu. Toggles the visibility of the second group of plot 
function items in the Evaluate selected models menu (see item 4). 
When all the items in the menu are visible, a check mark is shown at 
the More evaluation commands item. 


6. The Windows menu. The menu has items for each existing MATLAB 
figure window, and one item for the MATLAB command window. 
Selecting one of the items just brings the corresponding window to 
the front. 


7. The plot area. The prediction error variance measure is plotted as a 
function of the models (which models to include in the plot is 
controlled from the Show menu, item 2). 


8. A column of check-boxes used to select models to be recreated and 
saved to disk. The 20 models having the lowest value of the 
prediction error variance measure are included in the list and can be 
selected. 
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Table 3-11. Descriptions of the items of a model scan result window, 


Figure 3-24 (Continued) 


Item 


Description 


Columns displaying the values used for the parameters n, beta and 
gamma to create the 20 models having the lowest value of the 
prediction error variance measure. 


10. 


Column displaying the obtained prediction error variance measure for 
the 20 models having the lowest values. 


Editable fields where to enter the names of the files to save the 
recreated models in. The fields are initially filled in with names 
including date and time when the scan was finished, and the values 
used for the parameters n, beta and gamma. The name prefix 
containing the date and time can be changed for all the fields at once 
by entering a new string in the editable field at the bottom of the 
window (item 13). 


12. 


Push-buttons that when pushed bring up edit windows where 
descriptions of the models to recreate are entered. The descriptions 
are saved together with the models. 


13. 


An editable field where a uniform name prefix to be used for the 
models to recreate can be entered. The default prefix contains the 
date and time the scan was finished. 


The window contains a plot of the prediction error variance measure as a function of 
all the models included in the scan, and a list of the 20 models having the lowest 
value of the measure. The list is used to select models to be recreated and saved to 


disk. 


When selecting models to recreate, it is wise to take models with different value of 
n, as this is the system dimension. The parameters beta and gamma are just values 
influencing how the models are calculated but do not have any relevance for a 
created model (see also Section B.5.2 on page 213 in Chapter B, Configuration 


Tutorial). 
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Batch identification 


The tool has the possibility to create in one run a whole batch of models with 
different connection structures (with connection structure we here basically mean 
all user-supplied input to the subspace identification, that is, which signals to 
include together with assumed delays and scalings and the parameters n, beta and 
gamma). Section B.5.4 on page 221 in Chapter B, Configuration Tutorial describes a 
situation where this feature is useful. 


The connection structures to consider are loaded into internal memory using the 
different controls at the bottom of the main window of the Identification Tool (see 
Figure 3-22 on page 80). To have this part of the display shown, the Batch 
identification toggle button in the window (item 12 Figure 3-22 on page 80) must 
be pressed. 


There can be up to15 connection structures loaded. We can view the situation as 
having 15 slots available where we can insert connection structures. A slot is 
selected by clicking one of the radio buttons in the lower left corner of the window 
(item 13 Figure 3-22 on page 80). 


Only 5 radio buttons are shown at a time, but there is a pop-up menu (item 21 Figure 
3-22 on page 80) where any of three different sets (1..5, 6..10 and 11..15) of radio 
buttons can be selected for display. For each radio button there are corresponding 
edit fields to enter values for n, beta and gamma, an edit field to enter a name for 
the file where the identified model is saved and a push button that brings up an 
editable window where a description of the model can be entered. The first line of 
the description is displayed as the label of the button (except when no description is 
entered in which case the string “No description entered” is shown). 


The connection structure is loaded into memory by pushing the button labeled load 
a structure (item 20 Figure 3-22 on page 80). The data loaded is the entered values 
of n, beta and gamma, the file name and description of the model and the data 
entered in the signal selection part of the window (see Section, Model scan on page 
89 on how signals are selected and known delays and signal scalings are defined). If 
a structure already was loaded into the selected slot, a dialog prompting for the 
confirmation of the destruction of the old structure is brought up. 


When a radio button corresponding to a slot containing a connection structure is 
selected, the controls on the same row as the button and the signal selection part of 
the window changes to display the loaded structure. 
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Structures can be removed from the internal memory one at a time by selecting a 
radio button and then pushing the button labeled delete a structure (item 19 Figure 
3-22 on page 80). A dialog prompting for the confirmation of the deletion is brought 
up before the deletion takes place. 


Several structures can also be deleted at once by checking the check boxes that 
exist, one for each radio button (item 17 Figure 3-22 on page 80). When the 
structures to delete has been checked, the deletion is executed by pushing the Do 
button (item 18 Figure 3-22 on page 80). Before the deletion takes place, a dialog 
displaying a list of the slot numbers for the selected structures and prompting for the 
confirmation of the deletion is brought up. 


The Subspace identification > Batch identification menu! of the main window 
contains two commands, all structures and current structure. If the command all 
structures is selected, identifications corresponding to all connection structures in 
memory are performed, with the results saved in files with the given names. If any 
of the given file names equals already existing files, a dialog is brought up giving 
the possibility to cancel the destruction of the old file. 


If the command current structure is selected, only a model with the connection 
structure corresponding to the selected slot is created. 


Prediction error identification 


To start a prediction error identification, one of the models listed in the Prediction 
error identification > Start with model menu of the main window is selected. 


The list consists of all model files in the current application. Only models created by 
identification (subspace or prediction error) are possible to use. 


By selecting one of the files, the dialog shown in Figure 3-25 on page 97 is brought 
up. 


1. The menu is visible only when the Batch identification is pressed. 
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MIMO Identification Tool: Prediction error identification from : idmod1_P¥6. 


Name for the model : idmod1_PV6_p1 Max iterations : | 10 


Description 


This PEV model was based on: idmod1_PV6 . 


Press OK to start the prediction error identification 


OK | Cancel | 


Figure 3-25. Window from where a prediction error identification is started. 


The title of the dialog indicates which model that the search starts from (the one 
selected in the menu). There are edit fields to enter a name for the file where the 
result is saved and a description of the model. There is also a field to enter the 
maximum number of iterations allowed in the parameter search (see Section A.1.2 
on page 197). 


The file name field is initially filled in with a string constructed as 
<initModel>_p<i> where <initModel> is replaced by the name of the model that the 
search starts from and <i> is a number making the file name unique. 


The model description field is initialized by a string indicating that the model is 
created using the prediction error method and which model that was used as the 
initial one. 
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A default value, reasonable in most cases, is also given for the maximum number of 
allowed iterations. 


The calculations start when the OK button is pushed. If the given file name equals 
that of an already existing file, a dialog is brought up giving the possibility to cancel 
the destruction of the old file. 


When the new model is created, it is saved in a file with the given name. 
3.5.3.2 Model Evaluation 


Overview of the Model Evaluation Tool 


The Model Evaluation Tool is started from the Tools > Modeling > Evaluation 
menu of the Engineering Tools control panel (see Section 3.3 on page 33). 


When the tool is started the Model Evaluation control panel (Figure 3-26 on page 98 
with descriptions in Table 3-12 on page 99) and the General Plot Setup Interface 
window (Figure 3-27 on page 100 with descriptions in Table 3-13 on page 101) are 
brought up. The menus of the windows are described in Section, The menu bar of 
the Model Evaluation control panel on page 103 and Section, The menu bar of the 
General Plot Setup window on page 105. 


Evaluation of models BEI EI 
ile Plot Information Options Windows 


Evaluation data set : Update the popupmenu 


Figure 3-26. Model Evaluation control panel. 
The items are described in Table 3-12 on page 99. 


98 3BSE 011 572R101 


3dMPC User's Guide 
Section 3.5.3 Modeling Tools 


Table 3-12. Descriptions of the items of the Model Evaluation 
control panel, Figure 3-26 


Item Description 


1. The menu bar. It is described in Section, The menu bar of the Model 
Evaluation control panel on page 103. 


2. Pop-up menu selecting the data set used in the plots that need data 
in addition to the evaluated models. The menu is at the start of the 
tool initiated with all data files available in the current application. 


3. A push-button that when pushed updates the pop-up menu (item 3) 
with all data files available in the current application. This is 
necessary if new data sets are added during an evaluation session. 


The Model Evaluation control panel is from where all the different plots of the tool 
are created. 


The General Plot Setup Interface window is used to set up the default behavior of 
the plots and to administrate the matching between plotted models and colors and 
line-styles used in the plots. For the model/color/line-style matching to be consistent 
during a session, the window must not be closed. 
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Figure 3-27. General Plot ark ea ad Wanda 
The items are described in ‘1 
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Table 3-13. Descriptions of the items of the General Plot Setup 


Interface window, Figure 3-27 


Item 


Description 


The menu bar. It is described in Section, The menu bar of the 
General Plot Setup window on page 105. 


Column of check-boxes used to select which Process Variables to 
include in created plots. There are equivalent columns for the 
Manipulated Outputs and the Feed-Forward variables. 


Column of editable fields where step/impulse sizes to apply in the 
Manipulated Outputs are entered. The values are used when step-, 
and impulse-response plots are created. There is an equivalent 
column for the Feed-Forward variables. 


Pop-up menu that controls if confidence levels should be included by 
default or not in plots where that is possible (for models created by 
the prediction error method, confidence levels can be included in 
plots of simulated outputs, step-, impulse-, and frequency-response 
plots and in noise spectra plots). 


Editable field where to enter the default confidence level to use for 
confidence intervals. 


Pop-up menu that controls if the plots of the simulated or predicted 
outputs (see Section, Simulated/Predicted Outputs on page 116) 
should contain simulated or predicted outputs. 


Editable field visible only if the pop-up menu, item 6, selects predicted 
outputs. The entered (integer) value is used as the prediction length 
(as number of samples) to use in predicted output plots. 


Pop-up menu that controls if the gain or the phase should be plotted 
by default in the frequency plots (see Section, Frequency Responses 
on page 117). 
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Table 3-13. Descriptions of the items of the General Plot Setup 


Interface window, Figure 3-27 (Continued) 


Item 


Description 


Editable field where the default frequency range to use in frequency 
response and noise spectra plots is entered. The range is entered as 
a two-element MATLAB row vector (two values separated by space or 
comma surrounded by brackets), and the unit is fraction of the 
Nyquist frequency. 


10. 


Editable field where to enter default lengths (in number of samples) of 
step-, and impulse-responses. 


Editable field where to enter the maximal prediction length (in number 
of samples) to include by default in prediction error variance plots 
(see Section, Prediction Error Variances on page 117). 


12. 


Column containing the names of plotted models. The adjacent 
columns (item 13) gives the colors and line-styles that identifies the 
models in the plots. When a new model is plotted, its name is 
automatically inserted in a free row of the column. Models can be 
removed and inserted by purpose using the Models managing menu 
(see Section, The menu bar of the General Plot Setup window on 
page 105). If the window is closed and then recreated, the column is 
cleared, so that any new plot will assign the plotted model to the first 
row in the column. A closed window can be recreated through the 
File > Open the General Plot Setup Interface menu of the Model 
Evaluation control panel. 


13. 


Columns containing strings identifying colors and line-styles. Each 
row in the columns corresponds to a model name in the column, item 
12. The tool handles 40 different color/line-style combinations but 
only 20 are visible at a time in the columns. The pop-up menu, item 
14, controls which 20 to display. 


14. 


Pop-up menu controlling which 20 model/color/line-style 
combinations to display in the columns, items 12 and 13. 
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The menu bar of the Model Evaluation control panel 


Each heading corresponds to a menu label. 
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Delete model 

Delete model files from disk. 

Delete data 

Delete data files from disk. 

Parametric identification 

Brings up the Identification Tool (see Section 3.5.3.1 on page 79). 
Open the General Plot Setup Interface 


Opens the General Plot Setup Interface window, or if it is already open, 
brings it to the front. 


Exit 


Closes the Model Evaluation control panel and all open plot windows. A 
dialog with a warning message is brought up permitting to cancel the 
operation. 


<Plot function> 


<Plot function> can be any of the menu items Model Output, Prediction 
Error Variances, Step Response, Frequency Response, Impulse 
Response, Noise Spectrum and Poles/Zeros. Only the first three items 
are visible by default. In the Options menu (see below) there is an item 
Full plot that can be checked to make all the items available. 


New plot 
<List of available model files> 


The sub-menu has one item for each available model file. 
Selecting one of the items produces in a new window a plot of 
the selected type (depending on which plot function that is 
selected) using the selected model. 
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<List of open plot windows of this type> 


Below the New plot item there is one item for each open plot window 
of the selected type. 


<List of available model data files> 


The sub-menu has one item for each available model data file. 
Selecting one of the items adds to the selected window a plot of 
the selected type using the selected model. 


Information 
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Information about the models 


<List of available model data files> 


The sub-menu has one item for each available model file. Selecting 
one of the items brings up a window showing some information 
about the model. The shown information is: the sampling interval for 
which the model is valid, the names of all signals in the file, any non- 
linear transformations applied to the signals, any textual description 
saved together with the model and the history of the model (a list 
describing how the model was created). 


Information about the data sets 


<List of available data files> 


The sub-menu has one item for each available data file that can be 
used for identification or model evaluation. Selecting one of the items 
brings up a window showing some information about the data. The 
shown information is: the sampling interval for the data, the names of 
all signals in the file, any non-linear transformations applied to the 
signals, the history of the data (a list of all the processing steps the 
data has gone through) and any textual description saved together 
with the data (descriptions can be entered when a log is started in the 
3dMPC On-line controller, see Section Identification Experiment 
window4.4.7 on page 177, and in the Pre-processing tool, see Section 
3.5.2 on page 42). 
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Options 
Full plot menu 
Toggles if all the plot functions should be available in the Plot menu or 


not (see the description of the Plot menu above). When all the items in the 
Plot menu are visible, a check mark is shown at this item. 


Windows 


Contains one menu item for each open MATLAB figure window and one item 
for the MATLAB command window. By selecting an item, the corresponding 
window is made active. 


The menu bar of the General Plot Setup window 
File 
Reload the configuration file 


Reloads the configuration file (see Section 3.5.1 on page 34) from disk. 
This is needed only if the file has changed and the changes are to be used 
during the current session. 


Close 


Closes the General Plot Setup Interface window. A dialog with a warning 
message is brought up when this command is issued. The dialog has the 
possibility to cancel the command. The reason is that the setup window must 
be open for the administration of the plots to work. 


Models managing 
Add models 
<List of available model files> 


The General Plot Setup Interface window administrates the binding 
between the different models and the colors used in the plots. This is 
done through a list with each entry in the list initially containing only 
color/line-style pairs (see Figure 3-27 on page 100, items 12 and 13). 
When a plot is created for a model never plotted before, the model is 
added to the first free entry in the list making all subsequent plots 
involving that model use the same color. The 
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Models managing > Add models menu contains one item for each 
model available in the current application. When one of the items is 
selected, that model is added to the color/line-style list, making it 
possible to assign a color to a model prior to any plotting of the 
model. 


Delete models 
<List of models in the color/line style list> 


The Models managing > Delete models menu contains one item for 
each model in the color/line-style list (see the Add models sub- 
menu). When an item of the sub-menu is selected, that model is 
removed from the color/line-style list. 


Windows 


Contains one menu item for each open MATLAB figure window and one 
item for the MATLAB command window. By selecting an item, the 
corresponding window is made active. 


How to use the Model Evaluation Tool 


Basic file handling in the Model Evaluation Tool 


The evaluation data set used in all the calculations that need process data is loaded 
via the pop-up menu of the Model Evaluation control panel (see Figure 3-26 on 
page 98, item 2). 


Any of the data or model files contained in the current application can be deleted 
from disk using the File > Delete data or File > Delete model menus of the Model 
Evaluation control panel. A list of the available files is presented in a pop-up menu 
where one file at a time can be selected and deleted. 


Information about the data and model files included in the current application can be 
shown by selecting a file name in one of the sub-menus Information about the 
models or Information about the data sets in the Information menu of the Model 
Evaluation control panel. 
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Plots 


All the plots that can be created in the Model Evaluation Tool are per model. A new 
plot is created either in a new window or in an existing one. To plot in a new 
window, one of the 

Plot > [plot function] > New plot > [model name] menus of the Model Evaluation 
control panel is used. Figure 3-28 on page 107 shows such a menu selection. To 


Evaluation of models 
File Information Options ‘Windows 
Eval Model Output > 


Prediction Error Variances > 


_ Step Response » . © 


Figure 3-28. Menu selection in the Model Evaluation control panel 
Creates a Prediction Error Variance plot in a new window. 


direct a plot to an existing window, one of the 
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Plot > [plot function] > [window name] > [model name] menus, as shown in 
Figure 3-29 on page 108, is used. 


Evaluation of models 


File Information Options ‘Windows 


Eval Model Output » 
Prediction Error Variances > New plot 
——— Step Response Prediction Error Variances 1 > 


Figure 3-29. Menu selection in the Model Evaluation control panel 
Creates a Prediction Error Variance plot in an existing window. 
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The general appearance of the windows where the different plots are drawn is 
exemplified with the step response plot window showed in on page 109. 


Yrange 
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Figure 3-30. Step response plot window created in the Model Evaluation Tool. 
Descriptions of the items in the window are given in 
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Explanations of the different items in the window are given in Table 3-14 on page 
110 and the different menus are described below. 


Table 3-14. Descriptions of the items in a step response 
plot window, Figure 3-30 


Item 


Description 


The menu bar, described in Section, The menu bar of plot windows 
on page 112 


Title giving the type of plot. 


The plot area. There is usually one plot area per transfer function, 
with the areas arranged so that each column corresponds to a 
process input (Manipulated Output or Feed-Forward variable, which 
one indicated at the top of the columns as OUT<i>, FF<i> 
respectively) and each row to a Process Variable (which one 
indicated to the left of the rows as PV<i>). Plots of 
simulated/predicted outputs and noise spectra only have one plot 
area per output and the plots of the prediction error variances only 
have one plot area. If the mouse is clicked inside a plot area, the 
content of that area is repeated in a new separate window where it 
also is possible to zoom in on a specific part of the plot by dragging 
the mouse over it. 


A button used to change the plotted x-interval. By pushing the button 
a dialog is brought up where the new interval is entered. The interval 
is entered as a two-element MATLAB row vector (two values separated 
by space or comma surrounded by brackets). The unit is the 
displayed time unit for step-, and impulse-response plots, plots of 
simulated/predicted outputs and plots of prediction error variances. 
For frequency-response and noise spectra plots the unit is 
radians/second and for pole/zero plots the unit is dimensionless. 


A button used to change the plotted y-interval. By pushing the button 
a dialog is brought up where the new interval is entered. The interval 
is entered as a two-element MATLAB row vector (two values separated 
by space or comma surrounded by brackets). 
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Table 3-14. Descriptions of the items in a step response 
plot window, Figure 3-30 (Continued) 


Item 


Description 


A button used to change the step and impulse sizes, with the current 
size displayed on the button. The button is only available for step-, 
and impulse-response plots. There are separate buttons for each 
process input. When pushed, a dialog is brought up where the new 
size is entered and the plot is redrawn with the new size applied. 


String displaying the steady-state gain of the transfer function (for 
step-response curves only). When plots of several models exist in a 
window, it is usually the gains of the last plotted model that is 
displayed. The gains for another model can however be displayed by 
selecting any of the items Color/line style, Show or Highlight of the 
Curves > [Model] menu of the plot window (see Section, Curves on 
page 112). 


Area for page selector toggle buttons. Only one button can be in its 
pushed state at a time. 


For models with more than 16 process inputs and/or 16 process 
outputs, the plots are divided into several pages. The first 16 rows of 
the transfer function matrix falls on different pages than the rest of the 
rows, and correspondingly for the 16 first columns of the matrix. 


There will be one toggle button for each available page, arranged ina 
matrix fashion indicating the position in the transfer function matrix of 
the transfer functions plotted on the corresponding page. By pushing 
one of the buttons, the corresponding page is displayed. 


When there is only one page available, no buttons are visible as in 
Figure 3-30. 
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The menu bar of plot windows 
Each heading corresponds to a menu label. 
File 
Close 
Closes the plot window 
Curves 
<List of models plotted in the window> 


The Curves menu contains one sub-menu for each model plotted in the 
window. 


Color/line style: <color> <line style> 


The item has no effect when selected (except for step-response 
windows where the selection makes the steady-state gains for the 
selected model being displayed). Its purpose is just to display the 
color and line style used for the selected model. 


info string 


Brings up a window showing some information about the model. The 
shown information is: the sampling interval for which the model is 
valid, the names of all signals in the file, any non-linear 
transformations applied to the signals, any textual description saved 
together with the model and the history of the model (a list describing 
how the model was created). 


Hide/Show 


Hides or shows the curves for the selected model. When a model with 
visible curves is selected, the item shows the string Hide and 
selecting the item makes the curves for the selected model invisible. 
When a model with invisible curves are selected, the item shows the 
string Show and selecting the item makes the curves visible. 


For step-response windows, when Hide is applied for curves of a 
model that have its steady-state gains displayed, the gains of the last 
non-hidden model in the Curves menu are displayed instead. If all 
models are hidden, no gains are displayed. When the Show command 
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is given for a model, the steady-state gains for that model will be 
displayed. 


Highlight 


Highlights (makes thicker) the curves for the selected model. If a 
model with highlighted curves is selected, a check-mark is shown to 
the left of the string Highlight and selecting the item de-highlights 
the curves. 


For step-response windows, when the curves of a model are 
highlighted, the steady-state gains of that model are displayed. 


Delete 


Deletes the curves for a model from the window. At the same time, 
the corresponding sub-menu of the Curves menu is removed. 


For step-response windows, if the deleted model had its steady-state 
gains displayed, the gains of the last non-hidden model in the Curves 
menu is displayed instead. If all models are hidden, no gains are 
displayed. 


Add confidence intervals 


This item is disabled for plots where confidence intervals cannot be 
created. If the selected model already has confidence intervals 
plotted, the Remove confidence intervals and Change confidence 
level items replaced this item. 


<ij> % 

There are four items with <ij> replaced by 95, 97, 99 and 99.9. 
By selecting any of the items, the confidence level to be used 
when calculating confidence intervals is set to the selected value 


and the confidence interval for the selected model is calculated 
and plotted. 


others 


By selecting this item a dialog is brought up where an arbitrary 
confidence level to be used when calculating confidence 
intervals, can be entered. After the entered value has been 
acknowledged by pushing the OK button of the dialog, the 
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confidence interval for the selected model is calculated and 
plotted. 


Remove confidence intervals 


This item is shown only for models that have confidence intervals 
plotted. When selected, the confidence intervals are removed and the 
item in the menu is replaced by the Add confidence intervals item. 


Change confidence level 


This item is shown only for models that have confidence intervals 
plotted. 


<ij> % 
There are four items with <ij> replaced by 95, 97, 99 and 99.9. 
By selecting any of the items, the confidence level for the 


selected confidence interval is changed to the selected value and 
the interval is replotted. 


others 


By selecting this item a dialog is brought up where an arbitrary 
confidence level to be used when calculating confidence 
intervals, can be entered. After the entered value has been 
acknowledged by pushing the OK button in the dialog, the 
confidence level for the selected confidence interval is changed 
to the entered value and the interval is replotted. 


Names of variables 


PV 


<Full PV name> 


The sub-menu PV of the Names of variables menu has one item for 
each Process Variable used in the plot, with the name of the signal 
used as the item label. Selecting an item has no effect. The only 
purpose is to display the full names of the signals as space limitations 
makes it unfeasible to display them directly in the window. 
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OUT 
<Full OUT name> 


The sub-menu OUT of the Names of variables menu has one item 
for each Manipulated Output used in the plot, with the name of the 
signal used as the item label. 


FF 
<Full FF name> 


The sub-menu FF of the Names of variables menu has one item for 
each Feed-Forward variable used in the plot, with the name of the 
signal used as the item label. 


Options 


Automatic Y scaling 


When selected, the y-axis is rescaled so that all the curves fit into the plot. 
This may be needed when curves for a new model are added to an existing 
plot. 


For the plots of poles and zeros (see Section, Poles and Zeros on page 
118), the item is replaced by an Automatic XY scaling item that when 
selected scales both the axes simultaneously. 


Phase/Gain 


This item is available for the Frequency Response plots only (see Section, 
Frequency Responses on page 117). For plots displaying the Gain 
functions, the item label Phase is shown and selecting it replaces the Gain 
plots with plots of the Phase functions. The opposite holds for plots 
displaying the Phase functions. 


Show steady-state gains 


This item is available for step-response plots only. When the item is 
selected, the steady state gains of one of the models plotted in the window 
are either displayed or hidden in a toggle fashion. When the gains are 
displayed, a check-mark is displayed at the item. When a new step- 
response window is created, the gains are displayed by default. 
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Windows 


Contains one menu item for each open MATLAB figure window and one item 
for the MATLAB command window. By selecting an item, the corresponding 
window is made active. 


Simulated/Predicted Outputs 


Plots where the measured outputs contained in the evaluation data set is drawn 
together with either the simulated or predicted model outputs can be created using 
the commands in the Model output menu (see the Section, Plot on page 103). 


The simulations and predictions use all the inputs available for the models (the 
predictions use also the outputs), so the plot window will have only one column of 
plot areas (cf. Figure 3-30 on page 109). The check-boxes in the General Plot Setup 
window (Figure 3-27 on page 100) used to select which Manipulated Outputs and 
Feed-Forward variables to include in plots thus have no effect on these plots. The 
check-boxes used to select Process Variables control which signals that are plotted 
but do not influence how the plots are calculated (the output predictions always use 
all the outputs of the model). 


Another difference to Figure 3-30 is that the sub-menu labels in the Curves menu 
(see Section, The menu bar of plot windows on page 112) contains together with 
the file names also a numeric figure providing a measure of the agreement between 
the measured and model outputs. The measure is calculated as 


ny, N 
Quality measure = y? 4 y (9, (K) - yAk)) (3.1) 
fom 
i=1 'k=1 


where y, are the calculated (predicted or simulated) outputs and y, the outputs in 
the evaluation data set. 6; is the variance of process output #i as estimated from the 
evaluation data. 


The first item in the Curves menu does not have any sub-menus and only displays 
the color and line-style used for the plot of the measured outputs. 
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Prediction Error Variances 


Plots of the prediction error variances for a model can be created using the 
commands in the Prediction Error Variances menu (see the Section, The menu 
bar of plot windows on page 112). 


The plots show a measure of the prediction error variance as a function of the 
prediction length. For each prediction length, the measure is calculated as in 
Equation (3.1) on page 116. The predictions are calculated using all signals included 
in the model, irrespective of which that are checked in the General Plot Setup 
Interface window, but in the sum over Process Variables, only the checked ones are 
included. 


The created plot windows differ from the one in Figure 3-30 on page 109 in that 
they always consist of only one plot area. 


Step Responses 


Step Response plots of all transfer functions of a model can be created using the 
commands in the Step Response menu (see the Section, The menu bar of plot 
windows on page 112). 


Frequency Responses 


This function is available only if the Options > Full plot menu of the Model 
Evaluation control panel (see Figure 3-26 on page 98) is checked. 


Frequency Response plots of all transfer functions of a model can be created using 
the commands in the Frequency Response menu (see the Section, The menu bar of 
plot windows on page 112). The plots in the window displays either the Gain or the 
Phase functions, which one controlled through the Options > Gain/Phase menu of 
the plot window. By clicking with the mouse in one of the plot areas, the Bode plot 
(including both Gain and Phase) of the actual transfer function is displayed in a 
separate window. 
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Impulse Responses 


This function is available only if the Options > Full plot menu of the Model 
Evaluation control panel (see Figure 3-26 on page 98) is checked. 


Impulse Response plots of all transfer functions of a model can be created using the 
commands in the Impulse Response menu (see the Section, The menu bar of plot 
windows on page 112). 


Noise Spectra 


This function is available only if the Options > Full plot menu of the Model 
Evaluation control panel (see Figure 3-26 on page 98) is checked. 


Plots of the power spectra of the noise acting on the model outputs can be created 
using the commands in the Noise Spectrum menu (see the Section, The menu bar 
of plot windows on page 112). 


As one spectrum describes all noises acting on one output, the created plot windows 
will have only one column of plots (cf. Figure 3-30 on page 109). 


Poles and Zeros 


This function is available only if the Options > Full plot menu of the Model 
Evaluation control panel (see Figure 3-26 on page 98) is checked. 


The poles and zeros of all transfer functions of a model can be plotted using the 
commands in the Poles Zeros menu (see the Section, The menu bar of plot 
windows on page 112). 


The positions of the poles and zeros in the complex plane are displayed by putting 
marks at the corresponding positions of x-y plots. The marks are ‘x’ for poles and 
‘o’ for zeros. This means that only colors are used to distinguish the poles and zeros 
of different models (see the Models managing section of Section, The menu bar of 
the General Plot Setup window on page 105). As only ten different colors are used, 
some models will share the same color if more than ten models are plotted. 
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3.5.3.3 Model Merging 


Overview of the Model Merging Tool 


Sometimes it is desirable to merge model descriptions of parts of a process into a 
full process model. If all the partial descriptions are obtained through identification, 
the correct way to obtain the full process model is to use the possibility to merge 
data files prior to the identification (see Section 3.5.2.17 on page 73). Note that the 
merging of data files is possible also when the files are obtained through 
identification experiments performed on different parts of the process (that is, the 
files doesn’t need to contain exactly the same signals). 


If some of the sub-models on the other hand are obtained by other means (for 
instance as user defined models, see Section 3.5.3.4 on page 125), the merging must 
be done afterwards, thus the need for a separate Model Merging Tool. 


The Model Merging Tool handles the merging of two models only, so with more 
models, the construction of the final model is done in steps. 


Another limitation with the tool is that the result, in the form of a model file, is not 
complete enough to serve as an initial model for a prediction error identification 
(see Section, Prediction error identification on page 96) and it does not provide 
confidence intervals. The merging thus should be the last step in forming a model 
used in the design tool included in the AdvaControl MIMO Controller (see Section 
3.5.4 on page 130). 


The Model Merging Tool is started from the Tools > Modeling > Merge Models 
menu of the Engineering Tools control panel (see Section 3.3 on page 33). The 
appearance of the main window of the tool is shown in Figure 3-31 on page 120 
with descriptions of the different items in Table 3-15 on page 121. 
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Figure 3-31. Layout of the Model Merging Tool main window. 
The different items in the display are described in Table 3-15 on page 121. 
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Table 3-15. Descriptions of the items of the Model Merging Tool 


main window, Figure 3-31 


Item 


Description 


The File menu. It includes Save, Save as and Close commands 
used to save a created model and to quit the tool. 


The Windows menu. Contains one menu item for each open MATLAB 
figure window and one item for the MATLAB command window. By 
selecting an item, the corresponding window is made active. 


Radio buttons to activate “Signal mode” for a model. If “Signal mode” 
is activated, all transfer functions for that model associated with a 
specific signal can be selected/deselected at once (see item 10). The 
radio buttons have a slightly non-standard behavior in that one or no 
button may be selected) 


Indicators showing the colors used to indicate the state of the model 
selectors (see items 10 and 11). Each model is associated with one 
color given in two shades, one indicating that the model is selected 
(will be included in the created model) and one indicating that it is 
deselected (will be excluded from the created model). 


Button to initiate a model merge. 


Pop-up menus used to select the models to merge. The menus 
contain a list of all available models in the current application. 


An editable field used to enter a description of the final model. The 
description will be saved together with the model. 


Status information. When two models are merged, it is done in two 
steps where the first step consists of creating two models based on 
the two initial models, but with the deselected transfer functions 
removed. Those two models are then merged together to form the 
final model. The status lines here shows the number of states of the 
two models before and after the deselected transfer functions are 
removed. 


Status line showing the number of states of the combined model. 
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Table 3-15. Descriptions of the items of the Model Merging Tool 


main window, Figure 3-31 (Continued) 


Item 


Description 


10. 


Model selectors to select/deselect all transfer functions of a model 
associated with one signal. Which model the transfer functions 
belong to is determined by the radio buttons described in item 3 and 
is indicated by the colors of the indicators. 


The model selectors reflect the current state for the signal. If all 
transfer functions of the model for the specific signal are selected or 
deselected, the selector will have a single shade (see item 4), 
otherwise half of the selector will have the shade corresponding to 
selected transfer functions and the other half will have the shade 
corresponding to deselected functions. 


Single shaded selectors work in a toggle fashion. Clicking on a 
selector for a signal with no transfer functions selected, selects all the 
transfer functions for that signal and vice versa for a signal with all 
functions selected. 


Clicking on the half of a double shaded selector that have the shade 
corresponding to deselected functions will select all functions not 
already selected, and the selector will change into a single shaded 
one. If the other half of the selector is clicked the opposite happens. 


If “Signal mode” is not active (none of the radio buttons described in 
item 3 is selected) the signal related model selectors are not shown. 
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Table 3-15. Descriptions of the items of the Model Merging Tool 


main window, Figure 3-31 (Continued) 


Item 


Description 


11. 


Model selector for a specific transfer function. The position in the 
chart defines which transfer function (see item 12 and 13) 


If both models contain the transfer function corresponding to the 
selector, it will be double colored, with the shades of the colors 
showing if the transfer function will be included in the resulting model 
or not (see item 4). 


If the transfer function is included in only one of the models, the 
selector is single colored, with the shade showing if the transfer 
function is to be included in the combined model or not. If the transfer 
function is not included in any of the models, the selector is not 
shown. 


When a single colored selector is clicked, the corresponding transfer 
function is selected or deselected depending on its original state. The 
same is true when a half of a double colored selector is clicked. The 
difference is that the transfer functions for both models cannot be 
selected at the same time, so if a transfer function for one of the 
models is selected when the transfer function for the other model 
already was selected, the other transfer function will be deselected. 


The names of the process outputs associated with the transfer 
functions in the corresponding rows of the chart. 


13. 


The names of the process inputs associated with the transfer 
functions in the corresponding columns of the chart. 


(1) Usually, exactly one button in a group of radio buttons must be selected. 


How to use the Model Merging Tool 


Basic file handling in the Model Merging Tool 


The models to merge are loaded by selecting two files using the pop-up menus of 
the main display (item 6 of Figure 3-31 on page 120). The pop-up menus contain 
lists of all models available in the current application. The file lists are constructed 
when the tool starts and are updated when the tool saves a model. 
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A combined model is saved to disk by using the File > Save or File > Save as 
menus (item | of Figure 3-31 on page 120). When Save as is used, a dialog is 
brought up where the name of the file to store the model in, is entered. When the 
Save command is used, the model is saved in the same file used the last time a 
model was saved during the session. If no model has been saved yet, the dialog to 
enter a file name is brought up. 


When the tool is exited (using the File > Close menu) a dialog is brought up 
permitting a save (as done by the Save command described above) to be performed 
before exiting. The other alternatives available through the dialog are to exit without 
saving or to cancel the exit process, returning to the tool. 


Merging the models 


When the two models to merge has been loaded (as described in Section, Basic file 
handling in the Model Merging Tool on page 123), the transfer functions to include 
in the resulting model, and from which models they should be taken, are defined 
using the model selectors in the main window (items 10 and 11 of Figure 3-31 on 
page 120. 


When the resulting model is saved to disk as described in Section, Basic file 
handling in the Model Merging Tool on page 123, the actual calculations to merge 
the models are performed. A description to save together with the model can be 
entered in an edit field in the main window (item 7 of Figure 3-31). 


If it is desirable to check the number of states of the resulting model, that may differ 
depending on from which model the included transfer functions are taken, the 
combined model can be pre-calculated using the Combine push button (item 5 of 
Figure 3-31. The resulting number of states is then shown in the main window, item 
». 
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3.5.3.4 User Defined Models 


Overview of the User Defined Models Tool 


To be able to incorporate process knowledge into the modeling, a simple tool is 
provided that permits the creation and editing of models where each transfer 
function is a first order model with delay. 


The models created with the tool cannot be used as initial models for prediction 
error identification (see Section, Prediction error identification on page 96). 


The User Defined Models Tool is started from the Tools > Modeling > Edit First 
Order Models menu of the Engineering Tools control panel (see Section 3.3 on 
page 33). 


The main window of the tool is shown in Figure 3-32 on page 126 with descriptions 
of the items in Table 3-16 on page 127. The menus of the main window are 
described in Section, The menu bar of the User Defined Models Tool main window 
on page 128. 
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Edit first order model: NoName | |OF x] 


Figure 3-32. User Defined Models Tool main window. 
The items are described in Table 3-16 on page 127. 
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Table 3-16. Descriptions of the items of the User Defined Models Tool 


main window, Figure 3-32. 


Item 


Description 


The menu bar. Described in Section, The menu bar of the User 
Defined Models Tool main window on page 128. 


Button that brings up a small window (item 5) where the model 
parameters for one transfer function is entered. There is one button 
for each transfer function, with the buttons arranged so that each 
column represents a process input (Manipulated Output or Feed- 
Forward Variable) and each row represents a Process Variable. 


Editable field for the sampling interval. The transfer functions are 
defined in continuous time and the given sampling interval is used 
when the discrete time model is created. 


Editable field for measurement noise variance estimate. For each 
Process Variable, an estimate of the noise variance should be given. 
It is used for the observer design (see Section, 3.5.4.1 MPC Tuning 
Tool on page 130). 


Window where the definitions of one of the transfer functions are 
entered. The window is brought up when any of the transfer function 
buttons of the main window (item 5) is pushed. Each transfer function 
is defined as a continuous time first order lag with a pure delay?) and 
the window contains editable fields where the Gain, Time Constant 
and Delay are entered. The window also contains an OK button that 
closes the window after the transfer function parameters been saved 
in the main window’), an Apply button that just saves the 
parameters in the main window and a Cancel button that closes the 


window without saving the parameters. 
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The creation of the discrete time model and saving to disk is not performed until ordered via the 


File > Save or File > Save As... menus (see Section, The menu bar of the User Defined 
Models Tool main window on page 128). 
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The menu bar of the User Defined Models Tool main window 


File 


New 


Resets the definitions currently stored in the main window and changes the 
window title to NoName. If the definitions that are stored in the window are 
changed compared to when they were created by Open..., New or when the 
tool was started, a dialog is presented prompting the user to accept the 
discarding of the old definitions. If this is not accepted, no action is performed. 


Open... 


If the definitions that are stored in the window are changed compared to when 
they were created by Open..., New or when the tool was started, a dialog is 
presented prompting the user to accept the discarding of the old definitions. If 
this is not accepted, no action is performed. If the action is accepted, a file 
dialog window is presented where the user can select a model file from where 
to read model definitions. If the selected file is identified to contain a model not 
created by this tool, a warning message is issued and no action is performed, 
otherwise the model definition of the file is read into the main window. 


Save 


The model definitions currently stored in the window are saved to the currently 
selected file. If no model is selected (through Open, Save As... or Save, or if a 
New has been performed) a file dialog is presented where the user can enter the 
file name. The model is saved, and the window title is changed to the name of 
the selected file. 


Save As... 


A file dialog is presented where the user can enter the name of a file to save the 
model definition in. The model is saved, and the window title is changed to the 
name of the selected file. 
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Exit 


If the definitions that are stored in the window are changed compared to when 
they were created by Open..., New or when the tool was started, a dialog is 
presented prompting the user to decide if the changed data should be saved to 
disk or not. The options are: 


yes 
The data is saved to disk and the tool is exited. 
no 
The tool is exited without saving the changed data. 
cancel 
No action is performed. 
Tools 
Edit Non-linear Functions 


The tool to edit and create definitions of non-linear transformations is started, 
permitting the creation of such definitions for any of the involved signals. See 
Section 3.5.2.10 on page 59 and Section 3.5.2.11 on page 64) for a description 
of this tool. 


How to use the User Defined Models Tool 


The tool can be used both to create new models from scratch and to edit models 
created earlier with the tool. 


The number of transfer functions in the model is determined by the configuration of 
the current application (see Section 3.5.1 on page 34). Each transfer function is 
defined separately in continuous time as a first order lag plus a pure delay (see 
Figure 3-32 on page 126, items 2 and 5). 


When the model is saved to disk (see Section, The menu bar of the User Defined 
Models Tool main window on page 128), the continuous time model is discretized 
with a user supplied sampling interval (see Figure 3-32 on page 126, item 3). 
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3.5.4 Tuning Tools 


Tuning specifications parameters are stored on disk files. The tools for MPC tuning 
consists of a editor for specification data files, and a tool for control system analysis. 


3.5.4.1. MPC Tuning Tool 


The MPC Tuning Tool is used to create and modify specification data. The user 
interaction is carried out through a number of windows. The windows in this section 
are described using an example with seven process variables, three manipulated 
outputs, and two feed-forward variables. 


The tuning tool consists of five windows: 


e The main window for definition of the model used for design, and for entering 
a description of the design. 


e The process variables tuning specification definition window. 
e The manipulated outputs tuning specification definition window. 
¢ The feed-forward variable tuning specification definition window. 


¢ The observer design definition window. 
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Main Window 


Figure 3-33 shows the main window for the MPC Tuning Tool. Its different 
components are explained in Table 3-17. 


Figure 3-33. Main window for the MPC Tuning Tool 


Table 3-17. Description of the items in the main window for the MPC tuning tool. 


Item no. Description 


1 Push-button to open a file browser to select a model for the 
design. The name of the chosen model is displayed on the button 
and the sampling interval is displayed below the button. 


2 Editable field to enter a description of the actual design. 
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Table 3-18 shows the pull-down menu functions for this window. The View-menu 
and the Tools-menu are identical for all five windows in the Tuning Tool. 


Table 3-18. Pull-down menu functions for the main window of the 
MPC Tuning Tool 


Menu 


Menu item 


Function 


File 


New 


Set default values for a new 
specification file. 


Open 


Prompt for a file name of an 
existing specification data file, to 
open the file for editing. 


Save 


Save the actual specification data 
file. Prompt for a file name if it is 
undefined. 


Save As 


Prompt for a new file name. Then 
save the actual specification in 
this file. 


Exit 


Close all windows of the tuning 
tool and exit. 
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Table 3-18. Pull-down menu functions for the main window of the 
MPC Tuning Tool (Continued) 


Menu Menu item Function 
View Main Design Raise the main design window 

Process Variables Design Raise the process variables 
tuning specification definition 
window. 

Manipulated Outputs Design Raise the manipulated outputs 
tuning specification definition 
window. 

Feed-Forward Variables Design | Raise the feed-forward variable 
tuning specification definition 
window. 

Observer Design Raise the observer design 
window. 

Raise Sub-Windows Raise all windows. 

Hide Sub-Windows Hide all windows, with exception 
for the main design window. 

Tools MPC Analysis Invoke the MPC analysis tool. 

Suggest Tuning Specifications Use an algorithm to provide initial 
tuning specifications. 

Options 3 DOF design Check mark to select one degree 
or three degree of freedom 
design. 


The process variables tuning specification definition window 


Figure 3-34 shows the window with tuning specifications for process variables in 
the MPC Tuning Tool. Its different components are explained in Table 3-19 below. 
Editable fields are used unless anything else is mentioned. Some fields are shaded 
because these parameters are not available in this particular example 
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Figure 3-34. Window with tuning specifications for process variables. 


Table 3-19. Description of the components in the window 
for process variables 


Item no. Description 
1 Pull-down menu to hide the window. 
2 Pull-down menus with menu-items as described in Table 3-18 
3 List of process variables in the actual configuration. Process 


variables that are not included in the selected model are dimmed. 


4 Pop-up menus to select the specModePV variables that define the 
highest ranked mode that could be selected. 


Check-box to indicate that a process variable is critical. 


Lower indices for the prediction horizons. 


Upper indices for the prediction horizons. 


CO;N!] O|]o 


Equivalent feedback control errors. (Inverse of penalty.) 
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Table 3-19. Description of the components in the window 
for process variables (Continued) 


Item no. Description 


9 Equivalent set-point control errors. (Inverse of penalty.) In case of 
a one degree of freedom design these fields are dimmed and 
updated with the contents of the fields described in item 8. 


10 Equivalent feed-forward control errors. (Inverse of penalty.) In 
case of a one degree of freedom design these fields are dimmed 
and updated with the contents of the fields described in item 8. 


11 Dead-band for the feedback control 


12 Lower values for soft limits for process variables. 
13 Priorities for lower soft limits for process variables. 
14 Upper values for soft limits for process variables. 
15 Priorities for upper soft limits for process variables. 
16 Lower values for soft limits for deviations. 

17 Priorities for lower soft limits for deviations. 

18 Upper values for soft limits in for deviations. 

19 Priorities for upper soft limits for deviations. 

20 Lower limit for set-point speed (negative). 


21 Upper limit for set-point speed (positive). 


22 Pop-up menus to select the response type for the set-point filters. 
Available response type are Normal (C=0.7), Damped (¢=1.0), and 
Off (no filter used). 


23 Rise times for the set-point filters. A rise time less than the 
sampling interval means that the set-point filter should be omitted. 
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The manipulated outputs tuning specification definition window 


Figure 3-35 shows the window with tuning specifications for manipulated outputs 
in the MPC Tuning Tool. Its different components are explained in Table 3-20 
below. Editable fields are used unless anything else is mentioned. Some fields are 
shaded because these parameters are not available in this particular example. 
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Figure 3-35. Window with tuning specifications for manipulated outputs. 


Table 3-20. Description of the components in the window 
for manipulated outputs 


Item no. Description 

1 Pull-down menu to hide the window. 

2 Pull-down menus with menu-items as described in Table 3-18 

3 List of manipulated outputs in the actual configuration. 
Manipulated outputs that are not included in the selected model 
are dimmed. 

4 Pop-up menus to select the specModeOuT variables that define 
the highest ranked mode that could be selected. 

5 Check-box to indicate that a manipulated output is critical. 

6 Upper indices for the control horizons. 
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Table 3-20. Description of the components in the window 
for manipulated outputs (Continued) 


Item no. Description 
7 Equivalent feedback control move. (Inverse of penalty.) 
8 Equivalent set-point control move. (Inverse of penalty.) In case of 


a one degree of freedom design these fields are dimmed and 
updated with the contents of the fields described in item 7. 


AS) Equivalent feed-forward control move. (Inverse of penalty.) In 
case of a one degree of freedom design these fields are dimmed 
and updated with the contents of the fields described in item 7. 


10 Cost factors for manipulated outputs. 

11 Lower values for soft limits for manipulated outputs. 

12 Priorities for lower soft limits for manipulated outputs. 
13 Upper values for soft limits for manipulated outputs. 

14 Priorities for upper soft limits for manipulated outputs. 
15 Lower values for hard limits for manipulated outputs 

16 Upper values for hard limits in for manipulated outputs. 
17 Lower limit for speed of manipulated outputs (negative). 
18 Upper limit for speed of manipulated outputs (positive). 
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The feed-forward variables tuning specification definition window 


Figure 3-36 shows the window with tuning specifications for feed-forward variables 
in the MPC Tuning Tool. Its different components are explained in Table 3-21 
below. Editable fields are used unless anything else is mentioned. . 
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Figure 3-36. Window with tuning specifications for feed-forward variables. 


Table 3-21. Description of the components in the window 
for feed-forward variables. 


Item no. Description 

1 Pull-down menu to hide the window. 

2 Pull-down menus with menu-items as described in Table 3-18 

3 List of feed-forward variables in the actual configuration. Feed- 
forward variables that are not included in the selected model are 
dimmed. 

4 Pop-up menus to select the specModeFF variables that define the 
highest ranked mode that could be selected. 

5 Check-box to indicate that a feed-forward variable is critical. 
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Observer Design Window 


Figure 3-37 shows the window for observer design in the MPC tuning tool. Its 
different components are explained in Table 3-22 below. Editable fields are used 
unless anything else is mentioned. 
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Figure 3-37. Window for definition of observer design specifications. 


Table 3-22. Description of the components in the window 
for observer design 


Item no. Description 
1 Pull-down menu to hide the window. 
2 Pull-down menus with menu-items as described in Table 3-18 
3 Name of selected process model. 
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Table 3-22. Description of the components in the window 
for observer design (Continued) 


Item no. Description 


4 Selector for type of observer design. Available types are 

1: Restricted radius design with covariance R2 for model. 

2: Restricted radius design with user defined covariance R2. 

3: Kalman filter design with covariances R1, R12, and R2 for model. 
4: Kalman filter design with user defined covariances R1 and R2. 


5 Fields for restricted radius observer design parameters. Dimmed 
when observer type is 3 or 4. 

6 Magnitude of slowest open loop pole. 

7 Bandwidth in radians/s for slowest open loop pole. 


8 (*) Editable field for ratio between bandwidth for the slowest observer 
pole and for the slowest open loop pole. (Ratio between item no 9 
and item no 7.) 


9 (*) Editable field for bandwidth in radians/s for slowest observer pole. 


10(*) | Editable field for magnitude of slowest observer pole. 


11 Fields for user defined covariances. 


12 Editable field for user defined covariance R2. Dimmed when observer 
type is 1 and 3. 


13 Editable field for user defined covariance R1. Dimmed when observer 
type is 1,2, and 3. 


Note: (*) Changes in any of the fields identified by item no 8, 9, and 10 will cause 
appropriate update of the other two fields. 
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The MPC analysis tool is used to analyze the underlying linear control system 
design. The system is shown in Figure 3-38. Time and frequency responses can be 
generated for the transfer functions of the control system listed in Table 3-23. 


PV 


FF 
SP 
——| SP-control FF-control 
RS 
RC OUT 
FB-control |———™|__ Process m(5) 
-] ~ee 


Figure 3-38. Structure of the control system for linear analysis. 


Table 3-23. Transfer functions for linear analysis 


Section 3.5.4 Tuning Tools 


Transfer Function 


Outputs 


Inputs 


Open Loop Model (G) 


Process variables, PV 


Manipulated outputs, 
manOUT. Measurable 
disturbances, FF. 


Feedback controller (C) 


Manipulated outputs, 
OUT 


Process variables, PV. 


points. (Ggp) 


Closed loop response to set- 


Process variables, PV, 
Manipulated outputs, 
OUT * 


Set-points, SP. 


Closed loop response to 
measurable disturbances 


(Git) 


Process variables, PV, 
Manipulated outputs, 
OUT * 


Measurable 
disturbances, FF. 
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Table 3-23. Transfer functions for linear analysis (Continued) 


Transfer Function 


Outputs Inputs 


Complementary sensitivity 
function (T) 


Process variables, PV, | Process variable 
Manipulated outputs, disturbances, RC. 
OUT * 


Sensitivity function (S) 


Process variables, PV, | Process variable 
Manipulated outputs, disturbances, RS. 
OUT * 


Loop gain (L) 


Process variables, PV 
(after process) 


Process variables, PV 
(before controller) 


In the closed loop step responses, manipulated outputs are also considered as process 


outputs.This definition differs from the normal text-book definition. 


The following analysis functions are available in the MPC analysis tool. 


Table 3-24. Analysis functions. 


Analysis function 


Description 


Bode diagram 


Create a Bode diagram (amplitude plot) for each 
individual transfer function in the multivariable 
system. 


Nyquist diagram 


Create a Nyquist diagram for each individual 
transfer function in the multivariable system. 


Singular value frequency 
response 


Create a singular value frequency plot for the 
multivariable transfer function. 


Characteristic loci 


Plot the characteristic loci for the 
multivariable transfer function. 


Step response 


Create step responses for each individual transfer 
function in the multivariable system. 


Not all combinations of transfer functions and analysis function are meaningful. 
Only those that are useful are available. 
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Main Window 


Figure 3-39 on page 144 shows the main MPC analysis window. The window holds 
data for the analysis functions in terms of editable fields and pop-up menus. The 
user can modify these data to obtain the desired characteristics of the analysis. The 
different components in the window are explained in Table 3-25 below. 
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Figure 3-39. Main window for MPC analysis. 


144 3BSE 011 572R101 


3dMPC User's Guide 
Section 3.5.4 Tuning Tools 


Table 3-25. Description of the components in the main window 


for MPC analysis 
Item no. Description 
1 Check-boxes with process variables that are available in the 


configuration. The internal name of the signal is displayed. The 
signal is shown in the analysis plots if the box is checked. 


2 Pop-up menus to select the operator defined modes for the 
process variables. Available options are AutoSP, Predict, and Off. 


3 Editable field that defines the step amplitude for simulation of step 
responses in set-point, in the complementary sensitivity function, 
and in the sensitivity function. 


4 Editable field that defines the step length for simulation of step 
responses in set-point, in the complementary sensitivity function, 
and in the sensitivity function. 


5 Check-boxes with manipulated outputs that are available in the 
configuration. The internal name of the signal is displayed. The 
signal is shown in the analysis plots if the box is checked. 


6 Pop-up menus to select the operator defined modes for the 
manipulated outputs. Available options are Auto and Man. 


7 Editable field that defines the step amplitude for simulation of step 
responses for the open loop model. 


8 Editable field that defines the step length for simulation of step 
responses for the open loop model. 


9 Check-boxes with feed-forward variables that are available in the 
configuration. The internal name of the signal is displayed. The 
signal is shown in the analysis plots if the box is checked. 


10 Pop-up menus to select the operator defined modes for the feed- 
forward variables. Available options are Active, Predict, and Off. 


11 Editable field that defines the step amplitude for measurable 
disturbances. This is used both for simulation of open loop step 
responses and for simulation of closed loop step responses for 
steps in the measurable disturbances. 
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Table 3-25. Description of the components in the main window 
for MPC analysis (Continued) 


Item no. Description 


12 Editable field that defines the step length for measurable 
disturbances. This is used both for simulation of open loop step 
responses and for simulation of closed loop step responses for 
steps in the measurable disturbances. 


13 Pop-up menus for definition of specification files that are available 
in the pull-down menus for analysis. 


14 Pop-up menus for definition of line styles that correspond to 
particular specification files. 


15 Pop-up menus for definition of models that are available as 
processes for evaluation in the pull-down menus for analysis. 


16 Pop-up menus for definition of line colors that correspond to 
particular evaluation models. 


17 Axis definition for Nyquist plots. 


18 Axis definition for amplitude diagrams of Bode plots and singular 
value frequency responses. 

19 Frequency definition for frequency responses. 

20 Pull-down menu functions. See Table 3-26. 

21 Display of used time unit. 


Table 3-26 shows the first and the second level pull-down menu functions for the 
MPC analysis tool main window. Here also the File pull-down menu functions are 
described. 
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Apart from the File pull-down menu the first level menus in Table 3-26 correspond 
to the transfer functions defined in Table 3-23. For clarification, the symbolic 
definitions in Table 3-23 also appears in Table 3-26. 


Table 3-26. Pull-down menus for the main window of the 


MPC analysis tool 
First Level Second Level Function 
Menu Item Menu Item 
File Hide Hide the window. 
Exit Exit and close the window. 
Open Loop Bode Plot Bode-diagrams for the open loop model. 
mode! Nyquist Plot Nyquist-diagrams for the open loop 
(G) model. 
SV Frsp Plot singular value frequency responses for 


the open loop model. 


Step Response 


Plot step-responses for the open loop model. 


FB-Controller 
(C) 


Bode 


Plot Bode-diagrams for the controller. 


Nyquist Plot Nyquist-diagrams for the controller. 
SV Frsp Plot singular value frequency responses for 
the controller. 
SP-design Bode Plot Bode-diagrams for the closed loop 
(Gp) system from set-points to process variables. 
Nyquist Plot Nyquist-diagrams for the closed loop 
system from set-points to process variables. 
SV Frsp Plot singular value frequency responses for 


the closed loop system from set-points to 
process variables. 


Step Response 


Plot step-responses for the closed loop 
system from set-points to process variables. 
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Table 3-26. Pull-down menus for the main window of the 
MPC analysis tool (Continued) 


First Level Second Level Function 
Menu Item Menu Item 
FF-design Bode Plot Bode-diagrams for the closed loop 
(Ga) system from feed-forward variables to 
process variables. 

Nyquist Plot Nyquist-diagrams for the closed loop 
system from feed-forward variables to 
process variables. 

SV Frsp Plot singular value frequency responses for 


the closed loop system from feed-forward 
variables to process variables. 


Step Response 


Plot step-responses for the closed loop 
system from feed-forward variables to 
process variables. 


FB-compl-sens 
(T) 


Bode Plot Bode-diagrams for the complementary 
sensitivity function. 

Nyquist Plot Nyquist-diagrams for the complementary 
sensitivity function. 

SV Frsp Plot singular value frequency responses for 


the complementary sensitivity function. 


Step Response 


Plot step-responses for the complementary 
sensitivity function. 
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Table 3-26. Pull-down menus for the main window of the 
MPC analysis tool (Continued) 


First Level Second Level Function 
Menu Item Menu Item 
FB-sens Bode Plot Bode-diagrams for the sensitivity 
(S) function. 
Nyquist Plot Nyquist-diagram for the sensitivity 
function. 
SV Frsp Plot singular value frequency response 
diagram for the sensitivity function. 
Step Response | Plot step-response diagram for the sensitivity 
function. 
FB-loop-gain Bode Plot Bode-diagrams for the loop gain. 
(L) Nyquist Plot Nyquist-diagrams for the loop gain. 
SV Frsp Plot singular value frequency responses for 
the loop gain. 
Char Loci Plot characteristic loci for the loop gain. 


Further levels of the pull-down menus depend on selected designs and evaluation 
models in Figure 3-39. 
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Windows for Displaying Analysis Results 


The analysis results are displayed in graphical form. Different windows are used for 
the different analysis functions. All these windows have a similar menu bar for pull 
down menus, see Table 3-27. 


Table 3-27. First and second level pull-down menus for the analysis windows. 


First Level Second Level Function 
Menu Item Menu Item 
File Hide Hide the window. 
Close Close the window. 
Signals System Outputs | Information about the transfer function 
outputs defined by the check-boxes in 
Figure 3-39. 
System Inputs Information about the transfer function 
inputs defined by the check-boxes in 
Figure 3-39. 
Curves curve-1 Description of curves in the analysis 
window. 
curve-2 ; 
Under each curve-item there are menus 
for deleting and for modification of color 
and style of the curves for this transfer 
function. 


An example of a window with analysis results is shown in Figure 3-40. It shows the 
ten step responses for the extended complementary sensitivity function 
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Figure 3-40. Example window for step response analysis. 
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Chapter 4 Runtime Operation 


4.1 Operating Overview 


This chapter describes the interfaces concerning the on-line controller, that is the 
Start-up and Shut-down Interface and the Operator Interface. The Start-up and Shut- 
down Interface is used to Start and Stop a specific controller on the local machine. 
The Operator Interface is used to control a local or remote controller. For more 
details see below. 


4.2 Runtime Tutorial 
See Appendix B 


4.3 Start-Up and Shut-Down Interface 


The Start-Up and Shut-Down Interface is started from the windows start-menu. That 
is, press the Start button, then select Advant | 3dMPC, then select On-Line 
controller as shown in Figure 4-1. 


Documentation 
dh Engineering Tool 
& On-Line Controller 


& Operator Interface 


: “SS Enginecing 

a Compaq Information Center > (@, MIMO Development Project 
= rc cca 

[a Netscape Sinatllodits (3), Software Administration 


(s WinSDE 1.21 
| | New Office Document 


Figure 4-1. Start-Up and Shut-Down Interface in windows start menu 
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ig. 3dMPC Start-Up and Shut-Down 


Select an application and push start 
The Shell controller is running 

The simulator is running 
The Shell controller is OK! 


IV Simulate process 


Figure 4-2. Start-Up and Shut-Down Interface window 


When the Start-up and Shut-down window appears, see Figure 4-2, the user can 
select an application in the left field, see item 1, Table 4-1, of the window by 
clicking on the application name. If the application is running the Stop button will 
be enabled otherwise the Start button will be enabled. The current status of the 
controller is displayed in the middle field, see item 2, Table 4-1. When Start is 
pressed the controller will be started. If the process should be simulated internally in 
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the controller the Simulate Process check box should be checked, otherwise the 
user will be prompted for an OPC server to be used, see Figure 4-3. 


i, Select OPC Server |. {OE x} 


9 _______f¥ Read/write Asyne 


3 —_____=— OK | Cancel 


Figure 4-3. OPC server selection window 


All available OPC servers are displayed in the list box, see item 7, Table 4-1. When 
a Server is selected and the OK button is pushed the controller will be connected to 
the OPC server and then started. 


The last thing that is done at start-up of the controller is to set the password which is 
used when the operator wants to get Engineer Authority see Figure 4-4. 


ia, Enter password for the controller |. {OF x! 


Enter password Ok | 13 
11 —_____ Cancel | 14 


Confirm password 


— 


12 


Figure 4-4. Enter Password window 


If no password is entered the shift between Operator/Engineer Authority will not be 
protected by a password. 
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To shut-down a controller the user should select an application and then press Stop 
in the Start-up and Shut-down display, see Figure 4-2. The shut-down is then 
confirmed in a yes/no dialog which makes it harder to shut-down the controller by 


mistake. 


Table 4-1. Descriptions of the items of the Start-Up and Shut-Down 


and OPC Select interfaces, Figure 4-2, Figure 4-3, and Figure 4-4 


Item 


Description 


1. 


List of all the applications, controllers, that could be run on the 
machine 


Status field for the selected controller 


Stop button used to stop the controller 


Start button used to start the controller 


Close button used to close the window 


@| a} x] |p 


Simulation checkbox. When this is checked the controller will 
simulate the process internally in the controller. 


List of available OPC servers on the local machine that the controller 
could be connected to 


OK button to accept the OPC server 


Checkbox that determines if the connection between the 3d MPC and 
the OPC server should be asynchronous or synchronous 


Cancel button, exits the window without connecting to the OPC 
server 


Edit field to enter password in 


Edit field to confirm password in 


OK button, enters the password to the controller and closes the 
dialog 


Cancel button, closes the dialog without entering a password 
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4.4 Operator Interface 


The purpose of the Operator functions is to provide an interface permitting a process 
operator to supervise and interact with the 3dMPC and the controlled process in an 
easy and intuitive way. There are parts in the Operator Interface that only an 
operator with Engineer Authority can change. Such parameters are marked 
“Engineer only” in the description below. This means that the part is disabled or 
invisible for a regular operator, but visible when Authority is set to Engineer. 


The Operator Interface is started from the windows start-menu. That is, click the 
Start button, select Advant | 3dMPC and then select Operator Interface as shown 
in Figure 4-5. 


fS 3dMPC Documentation 
(3) Engineering dh Engineering Tool 
y rc} MIMO Development Project & On-Line Controller 


rc} Software Administration & Operator Interface 
(is) WinSDE 1.21 


Figure 4-5. Operator Interface in windows start menu 


At the start-up of the interface the user has to connect to a controller on a specific 
computer, see Figure 4-6. The existing computers on the network are displayed in 
item 1, Table 4-2. To connect the Operator Interface to a 3d MPC controller, select 
the computer where the controller is executing. Press Find Controllers, item 3, 
Table 4-2. The running applications on that computer will be displayed in the 
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listbox in the lower part of the window, item 2, Table 4-2. When an application is 
selected press Connect, item 4, Table 4-2, to start the Operator Interface. 


‘. 3dMPC Controller Connection 


Local Machine 3 
58- Microsoft Windows Network eee 


Figure 4-6. Controller Connection window 


Table 4-2. Descriptions of the items of the Controller Connection window, 


Figure 4-6 
Item Description 

1. List of all the existing computers on the network 

2. List of 83¢MPCs executing on the particular computer see item 1 

3. Find Controllers button, displays the executing controllers, on the 
specified computer, in item 2 

4. Connect button, connects the Operator Interface to the selected 
controller and closes the window 

5. Cancel button used to close the window without connecting to a 
controller 


When the Operator Interface is connected to a controller, the Operator main 
window will be displayed, see Figure 4-9. The Operator main window is the main 
window of the controller and it contains the most important information of the 
controller and all the other windows can be opened by clicking in the Controller 
Windows tree, see item 5 in Table 4-3. 
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The Operator Interface is used to change and view values in the 3d MPC On-line 
controller. It contains several windows that manage different kind of information. 
All data that is changeable is presented in a button which, when pushed, displays a 
dialog in which the new value could be entered, see Figure 4-7. In cases where the 
new value is an arbitrary numerical value display a) in Figure 4-7 is used. In cases 
where the new value is restricted to be one of a group of predefined values, such as 
modes and priorities the list box in display b) Figure 4-7 is used to display the 
alternatives for the new value. 


a) ig. Set AuotSP for PV1 | - Ox! 


Variable Old Value New Value 


AutoSP [220 oo [220.0 Set | List | 


b) is. Set Mode for P¥1 lel Es 


Variable Old Value New Value 


ModePY i Off Jott ba Set | [ist | 


Figure 4-7. Change value and mode dialogs 


The buttons has the same meaning for both dialogs. When the Set button is pushed, 
the change is performed immediately. When the List button is pushed!, the change 
command is saved into a list permitting several changes to be performed 
simultaneously. When the Cancel button is pushed, the change is inhibited. To the 
left in the dialogs, a string indicating what parameter to change is displayed together 
with the old value of the parameter. 


1. This function is not yet implemented. 
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Figure 4-8 shows an overview of the windows in the interface. The windows and 
their functions are described later in this chapter. 


Controller window 


Backup and restore window 


Controller signals window 


Close-up window 


Web plot window 


Operator main window 


Identification experiment window 


Tuning window 


|_| Extended Trim Curve window 


Simulator window 


Figure 4-8. Graphical hierarchy of the Operator Interface 
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The function of each block in Figure 4-8 is summarized below. 
Operator main window 


The Operator main window is the main display in which the most important 
information of the controller is displayed. From this window all the other 
windows can be reached. 


Controller window 


This window shows the main data of the controller and the parameter 
scheduling information. 


Backup and restore window 


This window shows the warm start option and handles the backup and restore 
of the controller. 


Controller signals window 


This window shows all the process signals (Process Variables, Manipulated 
Outputs and Feed-Forward Variables). The window also contains GUI controls 
to change parameters associated with the signals, such as set points and modes. 


Close-up window 


This window shows one signal at a time and all the major information for that 
signal. The window also contains GUI controls to change parameters 
associated with the signals, such as set points and modes. 


Webplot display window 


This window shows the so called webplots of the Process Variables (measured 
process outputs), the Manipulated Outputs (the controller outputs) and of the 
control deviations. 


Identification Experiment window 


This window shows the conditions under which the start-up and identification 
experiments take place. The window also contains GUI controls to change 
parameters associated with the identification experiment setup. 
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Tuning window 


This window shows for each signal, all avaliable tuning parameters, such as 

prediction or control horizons, equivalent error or move intervals, limits, and 
priorities for limits. The window also contains GUI controls to change these 

parameters. 


Extended trim curve window 


This window shows one out of four configurable plots of process signals. Each 
plot can contain up to five signals and covers a specified time span, from the 
current time and backwards. The plots also contain the available predictions of 
the signals. 


Simulator window 


This window shows the simulator model, simulation sampling interval, levels 
for the Process Variables, Manipulated Outputs and Feed Forward variables 
and step size for Feed Forward variables. The window also contains GUI 
controls to change these parameters. 
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4.4.1 Operator main window 


This window is the main window for the controller, see Figure 4-9. It is the first 


Backup and Rebtore 
©) Controller Signals 
Process Variabl 
Manipulated Ou 
Feed Forwards 
+ Close Up 
+ Process Variabl 


(- Manipulated Ov 
Feed Forwards 
Identification Experim: 
Tuning 
Process Vaviabl 
Manipulated Ov 
Feed Forwards 
Web Plot 
> Extended Trim Curve 
First 
Second 
Third 
Fouth 


‘Waim Start: Disabled Authority: Engineer Time: Fri Oct 01 09:02:15 y 


Figure 4-9. Operator main window 


window that is displayed when the Operator Interface connects to a controller. It is a 
Multiple Document Interface (MDD window which means that all the other 
windows will be displayed inside this one. All the information in this window is 
described in Tat 3. Different Controller Windows to be viewed are selected 
from the tree in the left part of the Operator Main window. The window will be 
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displayed inside the Operator main window. In the Option menu the Authority can 
be changed between Operator and Engineer. To get Engineer Authority, which 
enables the user to change more parameters in the controller, the user has to enter 
the password for the controller. The Simulator menu is only available when internal 
simulation is selected at start-up of the controller. From this menu the user can view 
the Simulator window, see page 182. 


Table 4-3. Descriptions of the items of the Operator main window, Figure 4-9 


Item Description 
1. Controller Data. Shows the most important information of the 
controller 
2. File menu, containing the Exit submenu 
3. Option menu, containing the Authority and Decimals submenus. 


The Authority menu is used to select Operator/Engineer Authority. 
The Decimals menu is used to select the number of decimals to 
display in the Operator Interface. 


4. Simulate menu, containing the SetUp submenu (only displayed 
when internal simulation is chosen at the start-up) 

5. Controller windows tree. Shows all the windows available for the 
controller 

6. Status bar. Shows additional information of the controller 


7. Windows display area 
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4.4.2 Controller window 


The Controller window handles the main data of the controller and the parameter 
scheduling data, see Figure 4-10. The parameters in the window are described in 


1 2 3 A 3 i & 8 9 10 11 


is, Controller 


; Description ~ Parameter Scheduling ~ Active Controller 
Default configuration, crated a 
98-11-09 15:53:12 PS None om Off |pestart | File: default. mt ‘i 
Model: idmodt_P¥/6|p.mxm 
= Save 
Eee Tsamp: 10.00 
Interval Active Specification files Gave as | 
1 ane 
> Description 
~ Basic Tsamp —} 5 05 | = Select| none: Prediction horizon and weights by 
Specitied 10.00 s fa Select | hone 
Sos} 
Actual 10.01 s | eS) es 
~ Alarms —25_| ee 
Activate | i | 
| 


fcked 
| 


I> 16 


Figure 4-10. Controller window 


Table 4-4. When the parameter scheduling is on, the value of the parameter 
scheduling variable decides which controller that is executed, but when the 
parameter scheduling is off the operator is allowed to manually select one of the 
four controllers by clicking on the Select interval buttons, item 16, Table 4-4. 
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Table 4-4. Descriptions of the items of the Controller window, Figure 4-10 


Item 


Description 


1. 


Description of the application 


Basic sampling interval, specified and actual 


Parameter scheduling variable name 


Parameter scheduling variable value 


Parameter scheduling variable unit 


Parameter scheduling variable status, On/Off 


Parameter scheduling Start/Stop (Engineer only) 


Active controller file 


CO} CO) N) Oy) oa; |) wo) rh 


Active controller model 


st 
o 


Active controller specified sampling interval 


Open, Save and Save as for active controller (Engineer only) 


Indication alarm Blocked/Activate'") 


Ordering alarm Blocked/Activate (Engineer only)? 


Limits for parameter scheduling interval (Engineer only) 


Indication for active parameter scheduling interval 


Ordering of active controller 


Controller name used in the interval 


18. 


Description of active controller 


(1) Function not yet implemented. 
(2) Function not yet implemented. 
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4.4.3 Backup and Restore window 


This window handles the data that is used for the backup, restore and warm start of 
the controller, see Figure 4-11. All the parameters in the window are described in 


ij, Backup and Restore BEE 


> Backup File 
Close | 4 
1 —__ 5 [Backup_1 9990603_085214.mxb i 
Restore File 


2 — a5 [Backup_19990520_1 64822. mxb i Restore From File | 


Warm Start 


L [Disabled bd 


Figure 4-11. Back-up and Restore window 


Table 4-5. In order to save current controller to use at warm start the warm start 
option, item 3 in Table 4-5 has to be set to Pre-Defined. 


Table 4-5. Descriptions of the items of the Back-up and Restore window, 


Figure 4-11 
Item Description 
1. Backup the active configuration to the file in the edit field 
2 Restore the configuration from the file in the edit field (Engineer only) 
3. Warm start option (Engineer only) 
4 Close button to close window 
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4.4.4 Controller Signals window 


This window displays all the signals in the controller, see Figure 4-12, Figure 4-13 
and Figure 4-14. All major data for the first 16 Process Variables, Manipulated 
Outputs or Feed Forward variables can be viewed at the same time. If there are more 
than 16 Process Variables there will be a second tab for the 17th and the rest of the 
PVs. All data in the window is described in Table 4-6. 


= 
N 
WwW 
KR 
Nn 
fon 
~ 
oo 
Ne} 
— 
S 
— 
= 
ja 
N 


~/Olx 
Manipulated Outputs d Forward ariables 
Stats Cyitical Npme Upit Mode Signgls 
Act Opr PY WSP PredPY AutoSP ExtRef 
J ok Nw yl [Aus Auto | [050 [oso [ 050 050 |[ None 
(21) ok Nor x2 ny2 [Auto auto | [000 [ooo [ 000 [None 
ia] ok No ¥3 nya [Predict Auto os = [ 043 [ 043 [None 
4] ok No YA ny [Predict auto | [ ov? [077 [077 000 |[ None 
is] 0k No ¥5 ny [Predict Auto os [oar [ o37 [None 
is] ok No Ys ny6 [Auto auto | [005 [000 [005 ooo |[ None 
7] ok No oY ny? [Predict Auto 004 [008 [004 000 |[ Nore 


~ Modes for All PV 


1 3 Auto E1 Predict Of 


Figure 4-12. Controller Signal window, Process Variables 
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Figure 4-13. Controller Signal window, Manipulated Outputs 
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123 4 a 6 67 8 
a -|o]x 
Process \Vpriables | Manipulated Outputs 
Status Ciitical ~Npme + Mody Signpls 
Act Opr FF 
eal Ok No D1 nd Active Active 0.00 
[2] 0k No 02 nd2 Active _ Active 0.00 


Modes All FF 


i 3 Active Predict Ott 


Figure 4-14. Controller Signal window, Feed Forward Variables 


Table 4-6. Descriptions of the items of the Controller Signal window, 
Figure 4-12, Figure 4-13 and Figure 4-14 


Item Description 


1. Signal number. Each group (PV, OUT and FF) numbered 
independently. The numbers are used in the web plots to identify the 
different rays. If a button is pushed the corresponding signal’s Close- 
up display (see Section Close-up window on page 172) will be 


displayed. 

2. Signal status. A strings “OK” or “Fail”. 

3. Indication for critical variables (see Appendix B.6.1, Configuration 
Tutorial) 

4. Signal name 
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Table 4-6. Descriptions of the items of the Controller Signal window, 
Figure 4-12, Figure 4-13 and Figure 4-14 (Continued) 


Item Description 
5. Signal value unit 
6 Indication of the current mode for the signal 
7 Indication of the mode chosen by the operator for a signal. 
8. Numeric value for the signal 
9 Indication of the Working Setpoint (the setpoint currently used by the 


controller) for a Process Variable 


10. Indication of the predicted value at the end of the prediction horizon 
for a Process Variable 

11. Indication of the Setpoint chosen by the operator for a Process 
Variable. 

12. Indication of the external setpoint for a Process Variable. 

13. Buttons to set all Signals of the selected type to a specific mode. 

14. Indication of the operator chosen OUT value in Man mode. 

15. Indication of the Manipulated Output value in E3 mode. 

16. Indication of the current preferred OUT value for a Manipulated 
Output. 

17. Indication of the operator chosen preferred OUT value for a 


Manipulated Output. 


18. Indication of the externally provided preferred OUT value for a 
Manipulated Output. 
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4.4.5 Close-up window 


This window handles all the non specification data for the 
Figure 4-16 and Figure 4-17. All the data in the window is 


signals, see Figure 4-15, 
described in Table 4-7. 


w=. Close-up [-[ol<J 
if Output Variables if Feed Forward Variables } 
Mode Variables 
2 Name Description Units RangeLo RangeHi iitical Mode Vi 
a Actual Operator _Extemal Spec 
%H) Vee ie wl “1.00 1.00 Be [Auto Auto [[-e1 [Auto 
Signal | Status and Supervision 
3 Status Supervision Bad Data 10 
pv WSP _PredPV_—_autoSP extfefl _dbandFB Lo Hi Ditt Actual Limit 
0.04 [000 [ 0.00 000 | [ None [ 000 Ok 0.95 085 7.00 i) 10 
Soft Limits for PV Alarm Limits for PY PvE 
4 Mode Actual Limits Operator Limits External Limits 
Actual Operator Lo Hi Lo Hi Lo Hi Hi 
Of ott__| [None _ [None None None None None None | __None_| 
;- Soft Limits for DEV ;-Alarm Limits for DEV 
cs) Mode Actual Limits Operator Limits External Limits 
Actual Operator Lo i Lo Hi Lo Hi Lo Hi 
Of oft__| [None [7 None None None None None None _|__None_| 
;- Speed Limits for Set Points, Units/s 
6 Mode Actual Limits Operator Limits Extemal Limits 
Actual Operator Lo Hi Lo Hi Lo Hi 
Off of | [None [” None None None None None 


7.68 


Figure 4-15. Close-up window, Process Variables 
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12 


‘=. Close-up 


L-[oy x] 
Process Variables if [ Feed Forward Variables 
OurT1.U1 ~ 
General ~ Mode Variables: 
Name Description Units RangeLo RangeHi  Ciitical Mode ManOutFF 1 1 
Actual Operator External Spec__—_ Enabled 
he Ci 
2 ae ba 00 [100 N° | | Awa | Auto |[ Auto [Auto| __ Yes 
Signal Set Points for Manipulated Output 
Refut 13 
OUT ManOUT ExtRef3 BalRef AddOut Mode Actual Operator External 
0.00 o.00 None None 0.00 Off 0.00 o.00 None 
Sofft Limits for OUT 
Mode Actual Limits Operator Limits Extemal Limits 
Actual Operator Lo Hi Lo Hi Lo Hi 
Off ot | [None [” None None None None None 
Hard Limits for OUT 
Mode Actual Limits Operator Limits External Limits 
Actual Operator Lo Hi Lo Hi Lo Hi 
Off ot | [None [~ None None None None None 
Speed Limits for OUT. Units/s- 
Mode Actual Limits Operator Limits External Limits Blocking 
Actual Operator Lo Hi Lo Hi Lo Hi Dec Inc 
Off on | [None [None None None None None 


Figure 4-16. Close-up window, Output Variables 
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1 
w=. Close-up [-[ol<J 
Process Variables Output Variables i iFeed Forward Variables 
ys | oe Description Units RangeLo RangeHi tical Meg See 11 
D1 Inter Reflux Duty nd [700 ["7007- = No [or of |[-Actve ~ [Active 
3 ~~ T es ‘oe Sem Hi Ditt iis Limit 10 
[oso ok [aos [os [10 [0 10 
| 
0 
Figure 4-17. Close-up window, Feed Forward Variables 
Table 4-7. Descriptions of the items of the Close-up window, 
Figure 4-15, Figure 4-16 and Figure 4-17 
Item Description 
Te Displayed variable 
2 General information for the variable (configuration data) 
3 Signal data 
4. Soft limit for the variables (Operator Limits, engineer only) 
5 Soft limit for the deviation (Operator Limits, engineer only) 
6 Speed limit for the variable (Operator Limits, engineer only) 
7 Alarm limit for Process Variable value (Engineer only) 
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Table 4-7. Descriptions of the items of the Close-up window, 
Figure 4-15, Figure 4-16 and Figure 4-17 (Continued) 


Item Description 

8. Alarm limit for Process Variable deviation (Engineer only) 

9. Trim curve for the variable 

10. Status and supervision information for the variable (Limit, engineer 

only) 

Ty, Mode information for the variable 

12. Hard limit for Output Variable (Operator Limits, engineer only) 

13. Setpoint information for the Output Variable 


4.4.6 Web Plot window 


This window is used to get an overview of the process that the controller is 
controlling. All the values and deviations of the Process Variables and all the values 
of the Manipulated Outputs are displayed in the same window, Figure 4-18. All the 


1 ae 1 5 4 6 67 8 


ia, Web Plot 


Figure 4-18. Web Plot window 
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data in the window is described in Table 4-8. If the mouse is moved to the text at the 
tip of an axes (item | and 6 in Table 4-8), the name of the variable will be displayed. 
If the text is clicked, the corresponding close-up window will be displayed. 


Table 4-8. Descriptions of the items of the Web Plot window, Figure 4-18 


Item Description 


1. Process Variable represented as a ray in the web plot 


Setpoint for the Process Variable represented with dashed line 


Value for the Process Variable represented with solid line 


Zero deviation for the Process Variable represented with dashed line 


Deviation for the Process Variable represented with solid line 


Manipulated Output Variable represented as a ray in the web plot 


NN] @) or) B) oO) hm 


Reference value for the Manipulated Output represented with a 
dashed line 


8. Value for the Manipulated Output represented with solid line 
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4.4.7 \dentification Experiment window 


This window is used to setup the Identification experiment, see Figure 4-19. All the 


Logging Data PRBS Data 
Status Running Status 
Loaging Sampling Interval (s} 10 PRBS Sampling Interval (s) W251 20 
Logging File Name | Stari EAES Stop PRBS | 1 9 
Logging Description - 4 
None 

Basic Data 

Basic Sampling Interval (s) 


2 16.1 


Figure 4-19. Setup Identification Experiment window 
data in the window is described in Table 4-9. In order to use the PRBS (17 and 19 in 


Table 4-9) to excite the process the user has to have the Engineer Authority. With 
Operator Authority this interface can only be used to log a process. 


3BSE 011 572R101 177 


3dMPC User’s Guide 
Chapter 4 Runtime Operation 


Table 4-9. Descriptions of the items of the Setup Identification 
Experiment window, Figure 4-19 


Item Description 
Ty Manipulated Output Variable name 
2. Manipulated Output Variable unit 
3. Manipulated Output Variable mode 
4. Manipulated Output Variable value 
5. Excitation signal value 
6. Amplitude for the PRBS signal (Engineer only) 
7. Manipulated Output Variable hard limit low 
8. Manipulated Output Variable hard limit high 
9. Manipulated Output Variable speed limit low 
10. Manipulated Output Variable speed limit high 
11. Logging status 
12. Description of log 


13. Start log (Engineer only if PRBS is on) 


14. Stop log (Engineer only if PRBS is on) 


15. Name of log file 


16. Logging sampling interval 


17. Start PRBS signal (Engineer only) 


18. Basic sampling interval 


19. Stop PRBS signal (Engineer only) 


20. PRBS sampling interval (Engineer only) 
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4.4.8 Tuning window 


This window is used to change specification data on-line, see Figure 4-20, 
Figure 4-21 and Figure 4-22. All the data in the window is described in Table 4-10. 


Figure 4-20. Tuning window, Process Variables 
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In order to change any of the data in this window the user has to have Engineer 
Authority. 


Figure 4-21. Tuning window, Manipulated Outputs 
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Process Variables Manipulated Outputs: 


FFI. D1 bd 
Feed Forward 
e) Index Name FF Unit Critical 
1 D1 0.00 nd1 No 
~ Mode 
4 Actual Spec 
Of Active 


Figure 4-22. Tuning window, Feed Forward Variables 


Table 4-10. Descriptions of the items of the Tuning windows, 
Figure 4-20, Figure 4-21 and Figure 4-22 


Item Description 


is Selected signal. This list box is used to select a signal of the type that 


window is used. 


is displayed. To select a different signal type, the tabs in the top of the 


2. General information of the signal. All values are boolean except the 


item 5. (Engineer only) 


Control Horizon High Limit that must be an integer greater than zero 
and less or equal to the largest Prediction Horizon for the PVs, see 


3. Set point filter. The filter Type and the rise time (Tr) can be set and the 


rise time must be greater than zero. (Engineer only) 
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Table 4-10. Descriptions of the items of the Tuning windows, 


Figure 4-20, Figure 4-21 and Figure 4-22 (Continued) 


Item 


Description 


Specification mode for the signal. (Engineer only) 


Prediction Horizon. These values are integers and the low value must 
be greater than zero and lower than or equal to the high value which 
when entered must be greater than or equal to the low value. 
(Engineer only) 


Equivalent Error Intervals. All values must be greater or equal to zero. 
(Engineer only) 


Dead band for the Process Variable value. (Engineer only) 


Speed limits for the signal value. The Lo limits must be less than zero 
and the High Limit must be greater than zero. (Engineer only) 


Soft limits for the deviation. The limits must be between zero and 
+(Range Hi - Range Lo) for the signal. The Priorities must be a value 
from 1 to 15. (Engineer only) 


Soft limits for the signal value. The limits must be between the signal 
ranges. The Priorities must be a value from 1 to 15. (Engineer only) 


Equivalent Moves. All values must be greater or equal to zero. 
(Engineer only) 


Cost, gives a penalty to the square of the difference 
(Manipulated Output value - RefOut). (Engineer only) 


13. 


Hard limits for the signal value. The actual limit is shown in red if the 
signal value is outside it. The limits must be between the signal 
ranges and inside the soft limits. (Engineer only) 


4.4.9 Simulator window 


This window handles all the data concerning the internal simulator, see Figure 4-23, 
Figure 4-24, Figure 4-25 and Figure 4-26. This window can be reached through the 
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Simulator menu in the Operator main window. All the data in the window is 
described in Table 4-11. 


, Manipulated eed Forward 
Process Variaples J ___.Qutputs._. Variables 
Name Unit Level 
1u1 nul 0.00 | 
2 U2 nu2 0.00 | 


0.00 | 


Figure 4-24. Simulator window, Manipulated Outputs 
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. Simulator Setup 


Level 


Step Size 


i nd 0.00 | 0.00 | 
0.00 | 0.00 | 


Figure 4-25. Simulator window, Feed Forward Variables 


Feed Forward 
Variables 


Simulator sample Interval 
[seconds] 


— =] 


Figure 4-26. Simulator window, Model 
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Table 4-11. Descriptions of the items of the Simulator window, 
Figure 4-23, Figure 4-24, Figure 4-25 and Figure 4-26 


Item Description 
1; Variable name 
2. Unit of the variable value 
3. Simulated signal level 
4. Simulated step level (added to level) 
5. Simulator model 
6. Simulator execution interval 


4.5 Runtime Operation Menus 


See Operator main window. 
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Chapter 5 Maintenance 


5.1 Error Messages 


5.1.1 On-line messages 


There are different types of messages coming from the on-line controller. The 
following types of messages are recorded in files in the application directory: 


¢ Messages are recorded in the file messages .txt. They describe the internal 
status of the algorithms in the on-line controller. 


e Alarms are recorded in the file alarm. 1ist. Occurrence of alarms can be 
noticed in Status field in the operator interface. 


¢ An OPC connection report is found in the file oPCConnection.txt after 
start up of the on-line controller. Here it is reported whether the on-line 
controller succeed or failed to connect to the configured OCS/DCS tags. 


Erroneous inputs to the operator interface are reported using standard Windows NT 
message-boxes. 


The internal simulator in the on-line controller produces messages and alarms 
similarly to an on-line controller connected to an OPC server. 


5.1.2 Off-line messages 


Error messages from the engineering tools are reported using standard Windows NT 
message-boxes. 


Some useful information for the experienced user may be available in the MATLAB 
console. 
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5.2 Disk-file Organization 


5.2.1 Product Directory Structure 
The product 3dMPC is located in the directory 


SAdvantRoot%\SystemModules\3dMPC\1.1\ 


where %AdvantRoot% was defined at the installation. This directory contains five 
sub-directories. They are described in Table 5-1. 


Table 5-1. Directory structure for 3dMPC product. 


Directory name Description 


Applications Application directories. Two examples are shipped 
with the product: PDistiller and Shell. 


Each application directory contains files and sub- 
directories. These are described in Table 5-2. 


bin Executable files. 

doc User’s Guide in PDF format. 
etc Miscellaneous files 

matlab MATLAB files. 


After installation, the environment variable sMPCApplicationsRoot% points out 
the Applications directory. This environment variable can be modified to point out 
another location for the applications. Note that the engineering license file must 
reside in the directory pointed out by the sMPCApplicationsRoot% variable. 
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5.2.2 Application Directory Structure 


Data for a 3dMPC application is stored in several files in a directory structure. The 
files and the directories are described in Table 5-2. The dot “.” refers to the 
application root directory, that has the directory defined by 
SMPCApplicationsRoot% as parent directory. . 


Table 5-2. Directory structure for a 3dMPC application. 


Directory Description of contents 


Configuration data file config.mxc. 
Application license file License.a 
Message files, described in Section 5.1.1, On-line 


messages. 
\backup Backup files * .mxb. 

Warm-start file warmstart .mxb. 
\data Data files * .mxd used for identification. 
\engineer Various files for engineering interface functions. 
A\log Log-files * . log from on-line logging. 
.\models Model files * .mxm. 
.\operator Various files for operator interface functions. 
\specifications Tuning specification files * .mxt. 
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A backup function allows that data, defined by the operator, can be stored on disk 
for later retrieval by the restore or the warm-start function. 


The on-line controller will have all manipulated outputs operating in manual after a 
restore or a warm-start. 


Table 5-3, Table 5-4, Table 5-5, and Table 5-6 list variables saved on disk in a 
backup. 


Table 5-3. General variables saved in a backup. 


Variable Description 
warmStartOpt Options how to restart controller. 
blockAlarms"") Block alarms for whole controller. 
parSchedOn Boolean that enables parameter scheduling. 
specFile1 Name of specifications file for scheduling interval 1. 
specFile2 Name of specifications file for scheduling interval 2. 
specFile3 Name of specifications file for scheduling interval 3. 
specFile4 Name of specifications file for scheduling interval 4 
parSchedLim12 Limit between scheduling interval 1 and 2. 
parSchedLim23 Limit between scheduling interval 2 and 3. 
parSchedLim34 Limit between scheduling interval 3 and 4. 
loggingTsamp Sampling interval for logging function. 
loggingDesc Text describing the logging. 
logFileName Name of file where log data is stored. 
prosTsamp Sampling interval for PRBS function. 


(1) Function not yet implemented. 
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Table 5-4. Variables related to Manipulated Outputs that are backed up. 


Variable Description 

oprModeOUT Operator mode choice. 

enbIManOUTff Enables use in the observer of the manipulated 
variable 

manOUT Operator defined value of manipulated variable. 

oprRefOUT Operator defined value of refOUT, used for static 
optimization. 

modeRefOUT Mode that determines how refOUT is defined 

oprModeLimOUT1 Mode that determines how soft limits are defined 

oprModeLimHardOUT | Mode that determines how hard limits are defined 

oprModeLimSpeedOUT | Mode that determines how speed limits are defined. 
See oprModeLimSpeedOUT. 

prosAmpl Excitation amplitude for PRBS excitation. 
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Table 5-5. Variables related to Process Variables that are backed up. 


Variable Description 

badDataLimitPV Allowed number of consecutive samples with read 
failure. 

oprModePV Operator mode choice. 

autoSP Operator defined set-point. 

oprModeLimPV Mode that determines how soft limits are defined: 

oprModeLimDEV Mode that determines how soft limits for the deviation 
are defined: 

oprModeLimSpeedSP | Mode that determines how set-point speed limits are 
defined: 

alarmLoPV Alarm limit for PV. 

alarmHiPV Alarm limit for PV. 

alarmLoDEV Alarm limit for DEV. 

alarmHiDEV Alarm limit for DEV. 


Table 5-6. Variables related to Feed-Forward Variables that are backed up. 


Variable Description 
badDataLimitFF Allowed number of consecutive samples with read 
failure. 
oprModeFF Operator mode choice 
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Appendix A Theory 


A.1 Parametric identification in 3dMPC 


The identification algorithms in 3d MPC finds the parameters of the discrete time 
state-space innovations model 


x(t+1) = Ax(t) + Bu(t) + Ke(t) 
y(t) = Cx(t) + e(t) 


where f is the time index, counted as number of samples, u the process inputs 
(Manipulated Outputs and Feed-Forward signals), y the process outputs (Process 
Variables) and e the innovations (the part of y that cannot be predicted from past 
data). The algorithms find the matrices A, B, C, K and an estimate of the covariance 
of e. 


(A.1) 


A.1.1 Subspace based identification 


This chapter will give a very brief introduction to the ideas behind the subspace 
based identification algorithm used in the Modeling Tools of 3dMPC. A detailed 
description of subspace identification is given in [VODM 96]. 


The basic idea behind the algorithm is to, from the input-output data, create state 
sequences x(t) and x(t+ 1) obeying a relation 


x(t-+ 1) = Ax(t) + Bu(t) + €,(2) 


i (A.2) 
y(t) = Cx(t) + €,(4) 


with the residuals €,(t) and €,(t) independent of x(1) and u(t) . The relation is then 


an ordinary linear regression from where unbiased estimates of the matrices A, B 
and C can be retrieved. 


The first step in the algorithm is to find the best linear predictor of the form 


yt) = Ly, (0) + Lu, (t) + Lyuft) (A.3) 
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where y(t) is a vector of predicted future process outputs, yO) , u(t) vectors with 
old process outputs and inputs respectively and u(t) a vector with future process 
inputs! The number of used old and future values are given by two user supplied 
parameters B and y and the vectors in (A.3) are formed as below 


ao y(t - B) 
yt) = | E+) y,(t) = y¢—-B +1) 
yety-V) at re 
u(t — B) u(t) 
ufo = [Me-B+D ufo = | Herd 


u(t — 1) u(t+y—1) 


The matrices L,, L, and L, determines how the old process inputs and outputs and 
the future inputs are weighed together to form the predictions. One of the 
predictions would be? 


B 
yt+i- 1) = YP (Ly y(t-(B-k+ 1D) +L,,ut-(B-k+ 1) + 
k=1 
Y 


>. L, wt+k—1) 
k=1 
To be of any use when trying to determine the matrices of (A.2), the predictor must 


be related to the state-space formulation. To find this relation we assume that x(t) is 


1.The predictor assumes future process inputs known. This is usually not the case, but the 
predictor is just an intermediate tool not used in practice, and as the identification is 
performed off-line, the data can always be organized to make the assumption fulfilled. 

2.With vector valued signals y and u, each Ly, (and correspondingly for L, and L, ) is a matrix. 
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the best estimate of the state vector based on y,(t) and u pt) . By iterating (A.2), 
setting future (unknown) residuals to zero, the predictor can be written as 


yt) = Vat) + Hu fo) 


C 0 0... 0 0 
a CB 0... 0 0 (A.5) 
ry = A, = CAB CBO 0 
cay! a _ 
(ers: rn CBO 


By comparing (A.3) and (A.5) we see that 
Pyr@ = Lyy, + Lou, (A.6) 


In (A.6), the right hand side is known (when the predictor is determined) but none of 
the terms of the left hand side is, not even the number of states (equalling the 
number of columns of ry ). 


Furthermore, if the state dimension is fixed, their is still an ambiguity in the 
representation, as any given state-space representation, with respect to its input- 
output relation, can be changed to an equivalent one by a simple change of state 
variables, which amounts to multiplying the state-vector from the left by an 
invertible transformation matrix. 


A multiplication of x(t) from the left could just as well be seen as a multiplication 
of I’, from the right. This freedom means that, relating to one valid I,, all other 


matrices spanning the same column space are valid alternatives, as you can go from 
one to another by multiplying it with the appropriate transformation matrix. To find 
this column space, (A.6) is repeated as many time instants as the data allows, with 
each equality put as one column in a large matrix as below 


Pyle xe) ...) = Lily, y,¢+D ..] +L2[u,0 u,¢+D.) AD 


I’, can now be found by finding a set of columns that spans the column space of the 
(known) right hand side of (A.7). 


A numerically sound way to find the column space of a matrix is through the 
singular value decomposition of the matrix. The number of linearly independent 
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columns of the matrix equals the number of non-zero singular values, and this could 
in theory be used to determine the number of states. The problem is that as we never 
will have found the exact predictor (A.3), no singular values will be exactly zero, 
and in practice, because the process from where the data is collected never obeys the 
assumption of being linear, finite dimensional and time invariant, their is usually not 
even a distinguishable threshold in the singular values. The only reliable way to 
decide an appropriate state dimension is to compare models of different dimensions 
against new data. 


In the 3dMPC algorithm, the state dimension must be supplied as an input. The base 
for the column space of I’, is then taken as the left singular vectors corresponding 
to the largest singular values, with the number of columns equalling the required 
number of states. 

With I’, known, the state sequence x(2) can be found by (pseudo)-inverting (A.7). 


As the states by construction are the best estimates based on available data (yO 


and u pit) ), the second relation of (A.2) is automatically satisfied (if €,(f) was not 


independent of x(t) , the estimate could have been improved). 


To proceed we now need a state sequence x(t + 1) that satisfies the first relation. To 
use the same sequence as above, with just having it shifted one step, will not do as 


each x(t + 1) will lack dependence of y(t—) and u(t— ) existing in x(t). Setting 


up the relation as in (A.2) will then not give €,(t) independent of x(t) . Instead we 


construct a new predictor that predicts y(t + 1) to y(t + y- 1) based on y(t) ; u(t) ; 
y(t) and u(t). By repeating exactly the same steps as before we get a state sequence 


x(t +1) that together with the earlier obtained x(f) sequence obey (A.2). 


By setting up (A.2) for all time instants possible, A, B and C are obtained by solving 
the equations in least-squares sense. A Kalman filter is then calculated based on the 
system matrices and using the statistical estimates of the covariances of the 
residuals €,(¢) and €,(t) as noise covariances. This gives us the K matrix of (A./), 
and as side result of the calculations also the covariance of the innovations. 


3BSE 011 572R101 


3dMPC User’s Guide 
Section A. 1.2 Prediction error identification 


A.1.2 Prediction error identification 


In the prediction error method, each candidate model is used to calculate the optimal 
predictor (given the candidate is the true model) and by using this predictor, the 
prediction errors are calculated by looping through the data. The model giving the 
smallest value of a norm of the prediction errors is chosen. 


With the model structure of (A.1), our use of the method can be summarized as: 


Find A, B, C, K minimizing 
N 


det( e(t)e (1) 
f=-1 
with 
e(t) = y(t) — Cx(t) 
x(t +1) = Ax(t) + Bu(t) + Ke(t) 


Large residuals e are damped to decrease the sensitivity to outliers and specific 
methods are used to handle the fact that the minimized function does not have a 
unique minimum!, 


A comprehensive treatment of prediction error identification is given in [Lju 87]. 


1.due to the non-uniqueness of the state-space representations 
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Appendix B Configuration Tutorial 


B.1 Introduction 


This tutorial is a step by step guide through the different phases of a controller 
design. The main steps can be summarized as: 


¢ Collect data revealing the dynamics of the process. 

e — Process the data to increase the relative level of usable information. 
¢ Identify a model from the data. 

e Design a controller based on the model. 


The material is organized with the purpose that the reader performs the different 
steps of the tutorial using the 3d MPC Engineering Tools. Most of the files that will 
be produced as intermediate results during such a session are included with the 
3dMPC distribution in the Applications\Shell catalog and can be used to verify that 
the exercises have been carried through correctly and to make it possible to skip 
some of the exercises. 
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The process model used in the tutorial is the so called Shell Heavy Oil Fractionator 
problem presented in [PrMor 87]. A schematic picture of the process is given in 
Figure B-1 on page 200. 


PV3 (A)PV1 Top draw 


DC > 
OUT1 


Upper reflux 
FF2 


PVv4 (T)— 


(A)PV2 Side draw 
> 


Pvé (T)}— 


Intermediate reflux 


Bottoms reflux 
OUT3 


Bottoms draw 


Feed 


Figure B-1. Schematic picture of the Shell process. 
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The main controlled variables are the product concentrations in the top and side 
draws. The main control variables are the top and side draw flow rates together with 
the heat duty of the bottoms reflux heat exchanger. Also available as measurements 
are the top temperature, the upper reflux temperature, the side draw temperature, the 
intermediate reflux temperature and the bottoms reflux temperature. Measurable 
and used as feed-forward variables are the heat duties of the intermediate and top 
reflux heat exchangers. 


All units used in the tutorial are normalized except the time unit that is in seconds. 
The time scale of seconds corresponds to minutes in the real process. The reason of 
this change is that it simplifies real-time simulations. 


Table B-1 lists the abbreviations used for the different signals throughout the 


tutorial. 
Table B-1. Abbreviated signal names used in the tutorial 
Signal Abbreviation 
Product concentration in top draw PV1 
Product concentration in side draw PV2 
Top temperature PV3 
Upper reflux temperature PV4 
Side draw temperature PV5 
Intermediate reflux temperature PV6 
Bottoms reflux temperature PV7 
Top draw flow rate OUT1 
Side draw flow rate OUT2 
Heat duty of bottoms reflux heat exchanger OUT3 
Heat duty of intermediate reflux heat exchanger FF1 
Heat duty of top reflux heat exchanger FF2 


Throughout the tutorial it is assumed that the Engineering Tools are running with 
Shell as the application (see Section 3.3 on page 33 of how to start the interface). 
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B.2 Control objective 


The objectives of the control system to be designed are: 


¢ Keep the controlled variables PV1 and PV2 as close as possible to their set- 
points while minimizing OUT3. 


This objective should be achieved subject to the constraints: 
e All the Manipulated Outputs must be in the interval [-0.5 , 0.5]. 


e The rate of change of all the Manipulated Outputs is limited to an absolute 
value of 0.05 units/sec. 


° PV1 must be in the interval [-0.5 , 0.5] 
° PV7 must be above -0.5. 


B.3 Identification experiment 
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The purpose of the identification experiment is to collect data that can be used to 
identify the process. To make certain that the data contains enough information, 
excitation signals are added to the ordinary controller outputs. The experiment is set 
up from the Identification Experiment Setup window in the 3dMPC Operator 
Interface (see Figure 4-19 on page 177). The additions to each of the controller 
outputs are Pseudo Random Binary (PRBS) sequences where the user specifies the 
interval with which new values are calculated and for each signal the amplitude of 
the sequence. The sampling interval used to log the experiment result is also defined 
in the window, independent of the PRBS interval. Both intervals must be multiples 
of the basic sampling interval used by the controller. 


The sampling interval for the log should preferably be chosen as the sampling 
interval that eventually will be used for the controller. 


The most important parameters to consider when choosing sampling interval are the 
desired bandwidth of the closed loop system and the bandwidth of the open loop 
system. The bandwidth of the closed loop system determines which disturbance 
frequencies that can be attenuated efficiently. It is normally a bad idea however to 
try to make the closed loop system significantly faster than the open system as this 
tends to make the solution very sensitive to modeling errors. 
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When designing the controller, a crucial parameter is the prediction horizon (see 
Section 3.5.4.1, MPC Tuning Tool). The controller outputs are calculated so that a 
criteria containing terms penalizing control deviations and controller output 
increments is minimized over the prediction horizon. For the minimization to be 
meaningful, the horizon must cover the important part of the impulse response of 
the system (so that the true effects of the calculated increments are known to the 
algorithm). The length of the horizon, in terms of sampling instants, increases the 
dimension of the minimization problem making it desirable to have it as short as 
possible. Reasonable values are in the interval 10 to 20 samples. The desire to keep 
down the prediction horizon in terms of samples while still covering enough time, 
puts a lower limit on the length of the sampling interval while the bandwidth 
demand puts an upper limit on it. 


The interval for the PRBS generation influences the energy content of the input. 
Each output sequence will be a signal switching between two levels, with a 
seemingly random time between the switches. The shortest length between two 
switches is the interval chosen for the generation, with considerably longer times 
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being rare. Figure B-2 shows a typical example. The interval should be chosen so 


i i 1 i i i 1 i i 
0 100 200 300 400 500 600 700 800 900 1000 
Time [s] 


Figure B-2. Typical PRBS sequence. 
The interval for the generation was 20 seconds. 


that the shortest steps give visible responses in the process outputs. 


As the discussion above makes clear, some knowledge of the process is needed to 
set up a successful identification experiment. In this tutorial, it will be assumed that 
the step response registrations shown in Figure B-3 on page 205 are available. 


The response time (time constant plus delay) ranges from a couple of seconds to 
more than a minute, so one sensibly designed controller probably couldn’t handle 
both the fast and slow dynamics (if possible, a process containing dynamics of very 
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different time scales, and where it is desirable to attenuate disturbances in several 
time scales, a cascade solution should be used). In this tutorial we will aim at a 
design concentrating on the slower dynamics. 


The file Shell\log\LOG_00.log contains simulated data of the same format as 
created by the 3d MPC On-line controller during an identification experiment. The 
data is sampled with an interval of 10 seconds and the PRBS sequences are 
generated with an interval of 20 seconds. To be able to identify the transfer 
functions from the feed-forward signals, square-wave signals were applied to those 
signals during the simulation. 
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Figure B-3. Open loop registrations of the step responses of the Shell process. 
All input stepsizes equal 0.2 . 
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B.4 Data pre-processing 


A successful modeling usually requires some amount of pre-processing of the data 
used in the identification. The engineering tools of the 3d MPC contains a Pre- 
processing tool for this purpose. The tool is started from the Tools menu of the 
Engineering Tools control panel. The main window of the Pre-processing tool is 
shown in Figure B-4 on page 206. 


Spectrum (rads) 

[aaa —] [A] 

= | yi 
i-| [i Ey 


Figure B-4. Main window of the Pre-processing tool. 
Before the tool can be used to process the data, the log file must be converted to the 


format used by the Engineering Tools. The conversion tool is launched from the 
Tools->Log-file import menu of the main Pre-processing tool window. When 
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launched, the conversion tool brings up the window shown in Figure B-5 on page 
207. 


Log-file Import 


log 
1999 4 14 13 15 2 


i] PRBS amplitudes=0.1. All OUT FI 


Figure B-5. Log-file import window. 


The file to convert is loaded by pushing the button in the upper part of the window. 
We select the file LOG_00.log. When loaded, the name of the loaded file is 
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displayed on the push-button and, if existing in the log file, a description of the log 
is shown in the text field in the lower part of the window!. 


Also shown in the window are the number of samples in the log and some statistics 
regarding the sampling interval. The statistics are the median value obtained for the 
sampling interval, and the values bounding the low and high quartiles. These values 
are based on time-stamps saved together with the logged samples in the file. Above 
these figures there is an editable field where the sampling interval to associate with 
the data set is entered. Usually, the median value is a reasonable measure of the 
interval with which the samples actually are taken, and should most often be used. 
But as the data used here is created in a simulation, the figure is just a measure of 
how fast the simulation program was able to run. Instead we should enter the 
simulated sampling interval (10 seconds). 


When the sampling interval is entered, the file conversion is performed by selecting 
the File->Save as... menu and selecting a file name in the brought up file dialog. 
The default file name (the same name as the log file but with the extension changed 
from .log to .mxd) is OK. When the data file is saved, the Log-file import window 
can be closed. 


The data file to pre-process is loaded into the tool through the File->Open menu in 
the tools main window by selecting the file in the file dialog brought up. When the 
file is loaded, the name of the file will be displayed in the title bar of the main 
window. 


To display the data, the signals to display must be selected. Up to five process 
outputs (Process Variables) and five process inputs (Manipulated Output and Feed- 
Forward Variables) can be simultaneously displayed in the small plot areas in the 
right part of the main window (see Figure B-4 on page 206). The signals are 
selected via the pop-up menus above the plots and are displayed in the form of their 
power spectra. One process input and one process output can then, by clicking with 
the mouse in the corresponding plots, also be displayed in the larger plot areas in the 
left part of the window. There are four larger plots, two that just are magnified 
versions of the smaller spectra plots and two that display the signal values 
themselves (time plots). 


Apart from just providing a more detailed view of the signal’, mean values or mean 
values and trends can be removed individually from the selected signals. This is 


1. The description is entered in the Identification Experiment Setup window of the Operator Interface when a 
log is started. 
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done through the pop-up menus above the plots. Usually, all the signals are 
processed uniformly as defined by the pop-up menus in the bottom left corner of the 
window (linear filtering must be performed uniformly as filters used on separate 
signals will turn up as parts of the identified process, or their inverses if applied on 
process inputs). The reason to allow means and trends to be removed from 
individual signals is that they can be the result of non-measurable disturbances that 
affects only some of the signals. 


A standard measure to take, if no trends need to be removed, is to globally remove 
the means from all the signals. 


If the data set is long enough, it is a good practice to split it into two parts, one to be 
used for the identification and the other used for evaluating the identified models. 


2. Even more detailed views can be obtained by clicking in a plot. This displays the plot in a separate window 
where it is possible to zoom in on specific parts by dragging with the mouse over the area of interrest. 
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To split the data, the menu command Tools->Split is used. The command brings up 
the Split File window showed in Figure B-6 on page 210. 


MIMO Pre-processing Tool: Split File LOG_00.mxd into 2 New Data Sets. 


Fill in the Sample time, and press OK. Current mouse mode: Zoom 


O 02 O4 O06 08 1 2 14 16 18 2 


Spit at Sample time: 


Figure B-6. Split File window of the Pre-processing tool. 


If any signals are selected in the main window of the Pre-processing tool they are 
displayed in the two plot areas of the Split File window. The signals to display can 
also be selected directly in the window using the pop-up menus above each of the 
plot areas. 


The point where to split the file is entered in the editable field below the plot areas 
as the time relative the first sample in the file. If the Return key is hit, a vertical line 
is shown in each of the two plot areas indicating the split point. If accepted, the OK 
button is pushed. This brings up a file dialog where the name of the file where to 
save the first part of the file is entered. When that file is saved, a new dialog to save 
the second part of the file is brought up. The files LOG_00_1.mxd and 
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LOG_00_2.mxd created by splitting LOG_00.mxd are available in the Shell\data 
catalog. 


Generally, the main purpose of the pre-processing is to perform linear filtering to 
accentuate certain frequency intervals of the data. The reason for this is that the real 
processes of interest usually are much more complicated than what can be described 
by the (reasonably) low-dimensional linear models that we want to use for the 
controller design. The design does not, however, rely on models being accurate in 
the whole frequency range. Modeling errors at frequencies above the bandwidth of 
the closed-loop system are not so important and the same can be said about 
modeling errors at very low frequencies, as they are accounted for by the integral 
action of the controller. The filtering makes it possible to fit a simpler model with 
the required accuracy in the limited mid-range frequency interval. 


Assume we like to efficiently attenuate disturbances up to frequencies of 0.05 
rad/sec (the bandwidth of a first order system with a time constant of 20 seconds). 
The model must at least cover a couple of octaves above that, so an upper limit of 
the filter bandwidth of 0.2 rad/seconds could be appropriate. This is quite fast 
compared with some of the dominating dynamics in the process (see Figure B-3 on 
page 205), so in order to be certain that also these parts of the dynamics are 
included, we choose at this stage! to use a low-pass filter. 


The type of filter is defined by the pop-up menu in the bottom left corner of the 
main Pre-processing window, and two editable fields to the right are used to define 
the breakpoints of the filter. For low-pass or high-pass filters, only one of the fields 
is used. The filter order can be set by clicking the filter order push-button in the 
bottom right corner of the window. Here we accept the default order. 


1. When evaluating a controller design, the analysis might show that better model quality is needed in certain 
frequency intervals, making it necessary to go back and create new filtered data files. 
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The plots shown in the main window includes both the filtered and the original data 
(see Figure B-7 on page 212). When a filter definition has been changed, the tool 


MIMO Pre-processing Tool - LOG_00_1f.mxd 


File Edit Tools Windows 
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Figure B-7. Main window of the Pre-processing tool displaying filtered data. 


directly recalculates the displayed plots. 


It is often a good idea to save the filtered data in new files, so that both the filtered 
and unfiltered data are available for model evaluation. The filtered versions of 
LOG_00_1.mxd and LOG_00_2.mxd are available in the Shell\data catalog as the 
files LOG_00_1f.mxd and LOG_00_2f.mxd. 
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B.5 Identification 


B.5.1 Introduction 


The identification tool included in the Engineering Tools contains two different 
identification algorithms, both producing state-space models (see Section A.1 on 
page 193). One is a prediction error method that iteratively searches the parameter 
space to find a model that minimizes a criteria based on the one step-ahead 
prediction errors obtained by the state-space model. 


The other is a so called subspace based method that without any search directly 
finds a state space model describing the input output relations of the data used for 
the identification. 


The reason to have both these methods available is that the subspace method gives 
biased estimates when applied to data collected in closed loop and that the 
prediction error search is rather sensitive to the initial parameter guess and is much 
slower than the subspace method. 


The methodology is to use the subspace method (even in closed loop) to decide on a 
suitable model structure in terms of system dimension and number of input and 
output signals and to find a preliminary initial model. This model is then used as a 
start for a prediction error search. 


B.5.2 Initial model scan 


The identification tool is started from the Tools->Modeling->Identification menu 
of the Engineering Tools control panel. 
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When the tool is started, the window showed in Figure B-8 on page 214 is brought 
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Figure B-8. Initial appearance of the MIMO Identification window. 


up. Initially most of the menus and controls of the window are disabled. To get them 
enabled, a data set to be used for the identification must be loaded. This is done via 
the 

File->Load identification data set... menu. Another data set to be used for a 
preliminary model evaluation may be loaded through the File->Load evaluation 
data set... menu. We load our pre-processed data LOG_00_1f.mxd and 
LOG_00_2f.mxd as identification data and evaluation data respectively. 


With an identification data set loaded, the disabled menus and controls of the 
window will be enabled. The three columns of check-boxes are used to select which 
signals to include in the identification. For each signal, any known delays can be 
entered (as number of samples), and a multiplicative scaling factor can be applied to 
the signal before the identification. The reason to use scalings is to get better 
numerical conditioning when the magnitude of the different signals differ very 
much. As an help in deciding appropriate scalings, the standard deviations of the 
signals, as estimated from the identification data, and the sizes of the signal ranges, 
as defined in the application configuration (see Section 3.5.1 on page 34) are 
displayed for each signal. We check all the signals to be included in the 
identification, but do not add any delays or apply any scalings. 


When starting from a new data set, we must use the subspace based method to get 
an initial model. The method takes three input parameters: n, beta and gamma. The 
parameter n is the number of states in the identified state-space model and beta and 
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gamma are used internally in the algorithm to create a preliminary predictor model. 
The predictor uses beta old process input and output vectors and gamma future 
process input vectors to predict gamma future process output vectors. Gamma also 
puts a limit on the number of states n in that n must be less or equal to gamma - 1. 


To aid in the selection of the parameters n, beta and gamma, an automatic scan over 
intervals of the parameters can be performed. The intervals are entered in the 
editable fields below the columns of check-boxes used to select the signals. The 
gamma values used in the scan can be chosen to be given as relative to the beta 
values, by checking the Add beta box, or independent to beta. The reason for this 
possibility is that gamma values close or equal to beta often have proven to be good 
choices. The number of combinations tried in the scan can in that way be drastically 
reduced. The tool can, after the signals to include in the identification have been 
selected, provide standard intervals.The standard intervals restrict the parameter 
combination to those having gamma = beta. With many signals included, the upper 
limit of the intervals may seem unreasonably large. A good advice is then to 
decrease it to something found more reasonable. If the result of the scan indicates 
that the best models are achieved for parameters at the upper boundary of the 
intervals, a new scan can be performed by increasing the upper limits again, but at 
the same time also increasing the lower limits to decrease the number of tested 
combinations. With the current example, the upper limits for the intervals suggested 
by the tool are 70 for n and 71 for beta (and gamma). We decrease both to 20. 


A scan is started from the Subspace identification->Scan n, beta and gamma 
menu of the MIMO Identification window. For each model created during a scan, a 
measure is calculated consisting of a weighted sum of the prediction error variances 
produced by the model using the evaluation data. The sum is taken both over the 
different Process Variables and over different prediction lengths. No models are 
saved during the scan, but the result is presented in the form of the window shown 
in Figure B-9 on page 216. 


The window consists of a plot of the prediction error variance measure over all of 
the scanned models and a list of the 20 models having the lowest values for the 
measure. Each row of the list contains a check-box used to select the model to be re- 
created and saved to disk. The rows also contain the values for the triplet n, beta, 
gamma, the value of the prediction error variance measure, an editable field where 
the name to use when saving the file to disk is entered and a push-button that when 
pushed brings up an edit window where a description to save together with the 
model is entered. 
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Variance 
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Figure B-9. Scan result window. 


Brought up as a result of a scan over n, beta and gamma in the MIMO Identification tool. 


216 


The file name fields are preset to names containing the date and time the scan was 
finished and the values for n, beta and gamma. The prefix containing date and time 
can be changed at once for all the fields by entering a new value in the editable field 
at the bottom of the window. 


When selecting models to create, it is a good practice to take models with different 
values of n, choosing the one having the lowest variance value in each group. 


The plot commands of the Model Evaluation Tool (see Section 3.5.3.2, Model 
Evaluation and Section B.5.3, Model evaluation) are directly available from the 
Evaluate selected models menu of the window. Selecting any of the commands 
creates a plot of the specified type containing all the selected models. 


The models are created from the Identification->Start subspace identification 
with saving menu in the window. 


3BSE 011 572R101 


3dMPC User’s Guide 
Section B.5.3 Model evaluation 


The Shell\models catalog contains a number of saved models named 
idmod_n<i>_b<j>_g<k>.mxm with i, j and k the values used for n, beta and gamma 
respectively. 


B.5.3 Model evaluation 


To examine the preliminary models created by the scan procedure described in 
Section B.5.2, Initial model scan, the Evaluation tool is used. The tool is started 
from the Tools->Modeling->Evaluation menu of the Engineering Tools control 
panel (or directly from the Evaluate selected models menu of the scan result 
window). It brings up two windows, the Model Evaluation control panel (see Figure 
B-10 on page 217) and the General Plot Setup Interface window (see Figure B-11 


Evaluation of models | — (OO x} 
File Plot Information Options Windows 


Evaluation data set : LOG 00 2t Update the popupmenu | 


Figure B-10. Control panel of the Model Evaluation tool. 


on page 218). 


In the General Plot Setup window, the default behavior of the plots created by the 
tool is defined. That includes selecting what signals to plot and some specifications 
regarding specific types of plots. When the window is first brought up, all signals 
are checked to be included in the plots. The window also administrates the colors 
and line-styles used to distinguish to which models the different plots belong. 


After the settings in the window has been completed, it can for convenience be 
hidden (minimized/iconified) but to keep consistent matchings of models and 
colors/line-styles it must not be closed. 


The creation of the plots are controlled from the Model Evaluation control panel. 
Some of the plot algorithms use in addition to a model also a data set. Which data 
set to use is defined through the pop-up menu in the panel. The menu contains all 
the data files available when the panel was first brought up. If data sets are added or 
removed during an evaluation session, the menu can be updated by pushing the 
button to the right of it. 
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General Plot Setup Interface 
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Figure B-11. General Plot Setup Interface window. 


The purpose of the model is to serve as a predictor in the controller, so a key feature 
is its ability to produce accurate predictions. To compare different models in this 
respect, the Prediction Error Variance plot is used. The plots are created per model, 
either in a new window or in an existing one. To plot in a new window, the 
Plot->Prediction Error Variance->New plot-><model name> menu of the Model 
Evaluation control panel is used. Figure B-12 on page 219 shows such a menu 
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selection. To direct a plot to an existing window, the 
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Figure B-12. Menu selection in the Model Evaluation control panel. 
Creates a Prediction Error Variance plot in a new window. 


Plot->Prediction Error Variance-><window name>-><model name> menu, as 
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shown in Figure B-13 on page 220, is used. The other types of plots are created in a 
perfectly analogous manner. 


Evaluation of models 
File Information Options ‘Windows 
Eval Model Output > 


—— Step Response Prediction Error Variances 1 > 


Figure B-13. Menu selection in the Model Evaluation control panel. 
Creates a Prediction Error Variance plot in an existing window. 


Figure B-14 on page 221 shows the Prediction Error Variances for the models 
created in Section B.5.2 on page 213 with the file Shell\data\LOG_00_2f.mxd used 
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as evaluation data. The figure indicates that the model idmodel1_n7b8g8.mxd has 


PREDICTION ERROR VARIANCES _ Prediction length = 250 s_ Sampling Time= 10 s Data set: LOG_O0_2f.mad 
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Figure B-14. Prediction Error Variance plots. 
Models created in Section B.5.2 on page 213. The dash-dotted line is model 
idmodel1_n7b8g8 


the best prediction capabilities. Before finally choosing a model though, it is 
advisable to plot also the step responses to rule out models with very strange 
behaviors, and to make a Model Output plot to get a feel of how well the data is 
explained by the model. 


B.5.4 Selection of feed-forward and support signals 


The 3dMPC has the ability to use measurable disturbances as feed-forward signals. 
In addition to this, it is also possible to connect Process Variables that are not 
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controlled, but only used to improve the estimation of the internal state of the 
process. This might drastically improve the ability to counteract non-measurable 
disturbances. 


Unfortunately it is often in practice hard to know both which measurable 
disturbances that are worth considering, and, even more so, which additional 
Process Variables (here called support signals) that could be of any use. 


In the Shell process, PV1 and PV2 are the controlled variables while PV3-PV7 are 
potential support signals. In addition, the measurable disturbances FF1 and FF2 are 
available as possible feed-forward signals. To determine the usefulness of these 
signals we identify new models where each of the models lack one of the signals. 
These models are then compared with the model chosen in Section B.5.3 on page 
217. If any of the models lacking a signal seems to perform just as well, or better, 
that is an indication that the full model has not picked up any useful information 
regarding that signal! 


To create the new models, a batch identification feature of the Identification tool can 
be used. Default, the MIMO Identification window (as seen in Figure B-8 on page 
214) has this feature hidden and to activate it, the Batch identification toggle- 
button in the lower left corner of the window is pushed. 


With the toggle-button pressed, a lower part of the window is folded out to form the 
window of Figure B-15 on page 223. The main part of the added window section 
contains five rows of controls. The rows constitute slots where input parameters to 
the subspace identification can be loaded. The window actually contains 15 slots, 
with five shown at a time. A pop-up menu in the lower right part of the window is 
used to select which five slots to show. 


To load a set of input parameters into a slot, it is first selected by clicking one of the 
radio-buttons. The values for n, beta and gamma and the name of the file where to 
save the model are entered in the edit fields to the right of the radio-button. The 
push-button to the right of the model name field brings up a window where a 
description to save together with the model can be entered. The signals to include in 
the identification, together with any known delays and wanted signal scalings are 
defined in the upper part of the window in the same way as was done in Section 
B.5.2 on page 213. The parameters are finally loaded by pushing the load a 
structure button in the lower right corner of the window. The slots containing a set 


1. The reason might be either that the signal is not important, or that the identification procedure failed in finding 
the (existing) relation. 
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MIMO Identification 


Figure B-15. MIMO Identification window with the 
batch identification feature activated. 


of parameters are indicated by the string “Structure Number <i> loaded” being 
displayed directly to the right of the corresponding radio-button. A slot can be 
cleared by selecting it with a radio-button and pushing the delete a structure 
button, or several slots can be cleared at once by selecting them with the check- 
boxes to the right of the slot rows before pushing the Do button directly above the 
column of check-boxes. 


Models corresponding to all the loaded slots are created by selecting the 
Subspace identification->Batch identification->All structures menu of the 
MIMO Identification window. If the 

Subspace identification->Batch identification->Current structure menu is 
selected, only the model corresponding to the parameters of the slot currently 
selected is created. 


The Shell\models catalog contains models named idmodel1_<xx>.mxm with <xx> 
replaced by a string indicating the missing signal. 


The models are compared with the same type of Prediction Error Variance plot used 
in Section B.5.3 on page 217. However, in Figure B-14 on page 221, the prediction 
error variances are calculated as a weighted sum including all the seven process 

output variables. The variables included in the sum are those checked in the General 
Plot Setup Interface window (see Figure B-11 on page 218). To compare the models 
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created here, we want only the controlled variables PV1 and PV2 to be included in 
the sum as it is the prediction of those variables that we want to improve by the use 


of feed-forward and support signals. Figure B-16 on page 224 contains the resulting 
plot. 


PREDICTION ERROR VARIANCES _ Prediction length = 250 s_ Sampling Time= 10 s Data set: LOG_O0_2f.mad 
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Figure B-16. Prediction Error Variance plots. 

Only the variances in PVI and PV2 are considered. The lowest solid line is for the 
model idmodel1_n7b8g8 containing all signals. All the other miss one signal with 
the lowest dotted line missing PV6. The upper solid line is for the model missing 
PV4 and the upper dash-dotted line is for the model missing PV3. 


The plot indicates that PV6 does not improve the predictions and that PV3 and PV4 
are the only signals that significantly improves it. 
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B.5.5 Prediction error identification 


When the model structure, in terms of number of states and included signals, is 
settled, and a preliminary model exists, this model is used as a start for a prediction 
error identification. The identification is initialized by selecting the 

Prediction error identification->Start with model-><model name> menu of the 
MIMO Identification window. This brings up the window shown in Figure B-17 on 
page 225. 


MIMO Identification Tool: Prediction error identification from : idmod1_P¥6. 


Name for the model : idmod1_PV6_p1 Max iterations : | 10 


Description 


This PE¥ model was based on: idmod1_PV6 . 


Press OK to start the prediction error identification 


OK | Cancel | 


Figure B-17. Window from where a Prediction error identification is started. 


The window contains a field where the name of the file where to save the model is 
entered, a field where a description to save together with the model can be entered, 
and a field where a limit on the number of iterations taken can be given. The name 
field initially contains a name based on the start model appended with the string 

““ p<i>” with i a number making the file name unique. The identification is started 
by pushing the OK button. 
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The file Shell\models\idmod1_PV6_p.mxm contains a model created through a 
prediction error search starting from the model idmod1_PV6. 


Usually it is a good idea to construct a number of models based on different data 
sets in order to be able to test the sensitivity of the controller to modeling errors. 


B.6 Controller design 


The controller design can be divided in two parts. In the first part, prediction and 
control horizons together with the weights on control deviations and Manipulated 
Output moves are specified. The tuning is done with the help of the MPC analysis 
tool where the linear system resulting if no limits or dead-bands are present, is 
analyzed. 


The second part consists of tuning/testing the effects of limits and dead-bands. This 
testing is done by running an instance of the on-line controller as a simulator (see 
Section Simulator window4.4.9 on page 182). This kind of simulation is also the 
only way (apart from tuning on-line) to tune costs on the Manipulated Output levels. 


In this tutorial, only one degree of freedom tuning will be considered. 


B.6.1 Linear tuning 


When a model is selected, the controller design is carried out in the MPC Design 
Tool. The tool is started from the Tools->Tuning->MPC design menu of the 
Engineering Tools control panel. 
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When the tool starts, the main window shown in Figure B-18 on page 227 is brought 
up. 


MPC Design (Shell): NoName 


<a 


Figure B-18. Main window of the MPC Design tool. 


A prerequisite to start a design is that a model is available. The model to use for the 
design is defined by pushing the push-button in the upper part of the window. Here 
we choose the model idmod1_PV6_p.mxm created in Section B.5.5 on page 225. 
When a model is chosen, four new windows are displayed. The windows are used to 
design the state observer! and to define the design parameters concerning the 
Process Variables, the Manipulated Outputs and the Feed-Forward variables 


1. The predictor used in the controller uses an estimate of the internal process state. The observer provides this 
estimate. 
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respectively. The windows are shown in Figure B-19 on page 228 to Figure B-22 


4 Observer Design 


Restricted-radius observer design using estimated measurement noise covariance R2 


Figure B-19. Observer design window. 


on page 231. 


There are two different observer design methods available, one is an ordinary 
Kalman design and the other is a method in which all filter poles are restricted to lie 
inside a circle with a specified radius. The poles of the system already lying inside 
the circle are left at their positions and the ones lying outside are being mirrored into 
the circle. The ordinary Kalman design uses two noise covariance matrices (process 
noise and measurement noise) while the pole restriction method uses only the 
process noise covariance. For both methods it can be chosen if the noise covariances 
should be taken from the model file (they are produced by the identification 
routines) or specified by the user. Which method to use and the source of the 
covariance matrices are selected via the pop-up menu at the top of the Observer 
Design window. The default is to use the pole restriction method with the noise 
covariance taken from the model file. This is recommended in most cases. 
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With the pole restriction method, the size of the circle can be defined in three 
different ways: as the radius of the circle, as the frequency corresponding to the 
breakpoint of a first order lag with a pole at the restricting circle or as the ratio 
between this frequency and the corresponding frequency for the slowest of the open 
loop poles. The specification is entered in any of the fields labelled “radius”, “speed 
ratio” or “bandw”’. When the contents of one field is changed, the other two are 
updated automatically. A small restricting circle gives a faster filter (the effect of 
process noise is detected quicker) but also a filter more sensitive to measurement 
noise and modeling errors. Usually, a value of the speed ratio between one and two 
is recommended. The default value is a radius of 0.8, and for our model that 
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Figure B-20. Tuning window for Process Variables. 


corresponds to a ratio of 1.77. 


The windows where design parameters are defined contains fields and pop-up 
menus where signal modes, prediction horizons, control horizons, weights to put on 
control deviations and moves in Manipulated Outputs, costs on Manipulated 
Outputs levels, control deviation dead-bands, limits and set-point filters are defined. 
Variables considered as critical for the function of the controller can also be marked 
in the windows (a Process Variable marked as critical must be in mode Auto or E1, 
a Manipulated Output in mode Auto and a Feed-Forward in mode Act for any 
Manipulated Output to be allowed to be in mode Auto) 


The signal modes defined at design time will restrict what modes that can be 
selected at run-time. Note that modes for signals not included in the model (PV6 in 
our case) are fixed to Off. We set PV1 and PV2 in Auto, the rest of the Process 
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MPC Design -- OUT 


File 


View Tools 


Variables (except PV6) in Predict (they are used only for state estimation, what we 
called support signals in Section B.5.4), all three Manipulated Outputs in Auto and 
both Feed-Forward signals in Act (used as Feed-Forward signals). 


The prediction horizons are given as the first and last samples (counted into the 
future) to consider in the optimization. The first sample to consider for a variable is 
chosen ideally as the beginning of the part of the response moving towards the 
steady state level!, and the last sample is chosen so that the important part of the 
impulse response (from all inputs) is covered by the horizon. 
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Figure B-21. Tuning window for Manipulated Outputs. 


The weights to put on control deviations are given as so called Equivalent Errors. 
The interpretation is that an error equalling the Equivalent Error gives a unity 
contribution (each sample of the prediction horizon) to the minimized loss function. 
Variables to be tightly controlled should therefore have small Equivalent Error 
values. The weights put on moves in Manipulated Outputs are analogously given as 
Equivalent Move Intervals (a move equal to an Equivalent Move Interval gives a 
unity contribution to the loss function). If the three degree of freedom feature is 
activated (by checking the Options->3 Degree of Freedom Design menu of the 
MPC Design main window) different weights can be assigned for the feedback, the 


1. that is, past any pure delay and any initial response in the opposite direction relative the steady state 
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set-point and the feed-forward designs. For a detailed description of the Design 
Tool, see Section 3.5.4 on page 130. 
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of Filmi 


Figure B-22. Tuning window for Feed-Forward variables. 


Reasonable initial values for the prediction horizons can be obtained through the 
Tools->Suggest Tuning Specifications->Prediction Horizon menu in any of the 
MPC Design windows. All the prediction horizons will be set to the same values 
chosen to include the important parts of the impulse responses from the 
Manipulated Outputs (with the impulse sizes scaled according to the ranges for the 
Manipulated Outputs). 
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Before the new horizons are applied, the small plot of the impulse responses shown 
in Figure B-23 on page 232 is presented. 


Suggested Prediction Horizon Acceptance | _ |Oy x] 


Accept | Cancel | 


Suggested Prediction horizon phi=21 


Normalized PV 


Samples 


Figure B-23. Window suggesting initial values for the prediction horizons. 


The window with the plot contains push-buttons permitting the acceptance or 
rejection of the suggested horizon. 


No corresponding tool to suggest the control horizon is available, but the control 


horizon is not dependent on the process to the same extent as the prediction horizon, 
and the default value (four) is often acceptable. 


Suggestions for the Equivalent Errors and Moves are obtained through the menu 
Tools->Suggest Tuning Specifications-> Weights. 


The file Shell\specifications\default1.mxt contains the controller specifications 
obtained by accepting the suggested prediction horizon and weights. 
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To analyze the design we start the analysis tool from the Tools-> MPC analyse 
menu in any of the MPC Design windows. This brings up the window of Figure 
B-24 on page 233. 


a979009 4 


Figure B-24. Main window of the MPC Analysis tool. 


The analysis tool has a collection of menus creating different types of plots (Bode, 
Nyquist, step-responses and so on) describing the design. For each plot, both a 
design specification and a model must be specified (except the plots concerning 
only the open-loop process or only the controller). This flexibility is needed to be 
able to test the sensitivity against modeling errors by connecting the controller to a 
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model other than the one it was designed for. What designs and models that should 
be available in the menus are controlled by a set of pop-up menus in the MPC 
Analysis main window where also colors and line-styles used to identify models and 
designs in the plots are specified. The model used for a design is always available 
(called Nominal Model in the menus) and if the MPC Design tool is open, the 
current design in the tool is available (called Actual Design in the menus). 


As in the model evaluation tool, Section B.5.3 on page 217, each plot can either be 
shown in a new figure or directed to an existing one. Figure B-25 on page 234 
shows an example of a menu selection producing a step-response plot for a 
Complementary Sensitivity function. 


fesse FB-sens FB-loop-gain 


— Bode » 
Nyquist > Hts Feed-Forwe 
SV Frey » 
PSSSERAS OSS SSE Step Step 


Nominal Model >| New Figure 


Step Response >} ActualDesign > 


i default! > = modi > Figure 8 
vu detuned1_1 E 
detuned1_2 


Figure B-25. Menu selection in the analysis tool. 
Creates a step-response plot for the Complementary Sensitivity function in a new 
window using the current design and its nominal model. 


With one degree of freedom, specifications are given for the feed-back design only, 
with the specifications automatically copied to the set-point and feed-forward 
designs. 


Only one degree of freedom design is treated in this tutorial, but even when carrying 
out a three degree of freedom design, it is usually advisable to start with the 
feedback design as that is the one determining the system stability. 


The purpose of the feedback is to compensate for non-measurable disturbances and 
model errors (with a perfect model and all disturbances measurable, pure feed- 
forward would do fine). To check the responses to disturbances, the Sensitivity and 
Complementary Sensitivity functions are used (see Section 3.5.4.1 on page 130). 
The Sensitivity function describes the response at the process outputs to 
disturbances also at the process outputs. The Complementary Sensitivity function 
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describes the responses to measurement disturbances, but by convention with 
reversed sign! 


Figure B-26 on page 235 shows the step responses of two transfer functions of the 
Complementary Sensitivity function with the default design (prediction horizon and 
weights as suggested by the tool) and the nominal model. 


Magnify FB-compl-sens-Step Response BEE bob Ld FB- Bscon sens- Shoe isso 


300400 BOD. 


Figure B-26. Complementary Sensitivity function step responses. 
Shows the responses for the default design. Left plot the response in PVI toa 
disturbance in PV], right plot the response in PV2 to a disturbance in PV2. 


One rule of thumb when designing a multivariable control system is to have it 
balanced in the respect that all Manipulated Outputs should be used to 
approximately the same extent, and the responses in the controlled variables should 
be approximately equally fast. Another rule is that it is usually not possible to 
design the feed-back to give a closed-loop system significantly faster than the open- 
loop system”. Here we see that the response in PV2 is considerably faster both than 
that in PV1 and the open-loop response (see Figure B-3 on page 205). 


The file Shell\specifications\detuned1_1.mxt contains a design where the 
Equivalent Error for PV2 has been approximately doubled (increased from .527 to 


1. The controller will compensate in the opposite direction of a measurement disturbance so in order to get a 
transfer function with positive gain, the minus sign is included. 
2. The set-point and feed-forward designs might very well be faster. 
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1.0). Figure B-27 on page 236 contains the corresponding step-responses for the 


| Magnify FB-compl-sens-Step Response | (OF x} Magnify FB-compl-sens-Step Response 


Figure B-27. Complementary Sensitivity function step responses. 
shows the responses for a design with increased Equivalent Error for PV2. The first 
row contains responses in PV1, the second row responses in PV2, the first column 
contains responses to a step in PV1, the second column responses to a step in PV2. 


Complementary Sensitivity function of this design and Figure B-28 on page 237 
contains the responses in the Manipulated Outputs. 
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mpl ponse -jo)xi| a ly F esponse 


Magnify FB-compl-sens-Step Response Pe) ES |) ¢ Magnify FB-compl-sens-Step Response |. | OF x! 


Figure B-28. Complementary Sensitivity function step responses. 
shows the responses for a design with increased Equivalent Error for PV2. The rows 
contain responses in OUTI, OUT2 and OUT3) respectively, the columns contain 
responses to a steps in PVI and PV2 respectively. 


(1) The responses in controller outputs are not included in the standard definition of the 
Complementary Sensitivity function. 


To test the sensitivity of the design to modeling errors, the controller is connected to 
a model other than the one used for the design. It is therefore preferable to create 
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models using different data-sets covering working conditions that should be 
handled. Here we have not created such models, instead we use the model 
Shell\models\mod1.mxm! to this purpose. 


Figure B-29 on page 238 shows responses for the same design as in Figure B-27 and 


Magnify FB-compl-sens-Step Response |. | OF x! Magnify FB-compl-sens-Step Response 


300—:—“‘x C8 


Figure B-29. Complementary Sensitivity function step responses. 
shows the responses for the design showed in Figure B-27 and Figure B-28. The left 
plot shows responses in PV1 to a step in PV] and the right plot shows responses in 
PV2 to a step in PV2. The solid line is the result with the controller connected to the 
nominal model and the dashed lines the result when connected to mod1. 


Figure B-28, both with the controller connected to the nominal model and 
connected to the model mod1.mxm. 


The responses indicate that the design probably is too sensitive. One standard way 
to detune an MPC is to decrease the control horizon, often all the way to 1. The file 
Shell\specifications\detuned1_2.mxt contains the design resulting if the weights of 


1. This is actually the model that created the log file used in the identification. 


3BSE 011 572R101 


3dMPC User's Guide 
Section B.6.2 Non-linear tuning 


detuned1_1 is used, but all control horizons are set to one. Figure B-30 on page 239 


Magnify FB-compl-sens-Step Response BEE Magnify FB-compl-sens-Step Response 


02 

; 0.15 

0.1 

0.05 

Q 
0 10 20 300 400 500 300 «400500 


Figure B-30. Complementary Sensitivity function step responses. 
The same plot as in Figure B-29 on page 238, but with all control horizons set to 
one. 


contains the responses corresponding to Figure B-29 obtained with this design. 


B.6.2 Non-linear tuning 


The linear analysis of Section B.6.1 does not consider the effects of limits put on 
Manipulated Outputs and Process Values or the effects of the static optimization 
performed if costs are put on Manipulated Output levels (see Section 3.5.4.1 on 
page 130). These effects must be checked by simulation. The simulation is 
performed by connecting the Operator Interface to an 3d MPC On-line controller in 
simulation mode (see Section 4.3 on page 153). 


The simulation should verify that the objectives of Section B.2 on page 202 are 
achieved in an acceptable way. The tuning basically amounts to set priorities to the 
soft limits and choosing costs for Manipulated Output levels. In our case we only 
have cost on one of the outputs, so the value of the cost actually doesn’t matter. 


When satisfactory behavior is found, the modified controller settings are saved (this 
can be done directly in the Operator Interface see Section 4.4 on page 157). 
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B.6.3 On-line tuning 


No matter how well tuned and analyzed the controller is, some surprises may 
always come when it is actually connected to the real process. Thus some fine- 
tuning on-line during commissioning is often needed. This tuning is performed 
equivalently to the non-line tuning of Section B.6.2. 
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Appendix C Technical Details 


This appendix summarizes some technical details for the 3d MPC product. 


C.1 Connections 


This section lists all variables that are accessed though the OPC server interface. 
The variables in the tables below are vectors unless not anything else is mentioned. 
Each component (vector element) can be connected to a specific tag in the 
OCS/DCS. The variables marked with * must be connected. 


Table C-1 lists variables related to manipulated outputs, that can be configured as 
inputs to the on-line controller. 


Table C-1. Inputs to 3d MPC related to manipulated outputs. 


Name Value Description 
extRef3 (range) External values for manipulated outputs. Used 
in E3-mode. 
balRef (range) External values for manipulated outputs. Used 
in Bal-mode. 
refOUT (range) External reference/set-points for manipulated 


outputs. Used in Ext-mode. 


loOUT1 (range) External values for low soft limits for 
manipulated outputs. Used in Ext-mode. 


hiOUT1 (range) External values for high soft limits for 
manipulated outputs. Used in Ext-mode. 


speedLoOUT (range) External values for low speed limits for 
manipulated outputs. Used in Ext-mode. 


speedHiOUT (range) External values for high speed limits for 
manipulated outputs. Used in Ext-mode. 


hardLoOUT (range) External values for low hard limits for 
manipulated outputs. Used in Ext-mode. 
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Table C-1. Inputs to 3dMPC related to manipulated outputs. (Continued) 


Name 


Value 


Description 


hardHiOUT 


(range) 


External values for high hard limits for 
manipulated outputs. Used in Ext-mode. 


blockDec 


0,1 


Boolean inputs. Prevent decrease in the 
corresponding manipulated output for 
components that are one. 


blockInc 


0,1 


Boolean inputs. Prevent increase in the 
corresponding manipulated output for 
components that are one. 


extModeOUT 


1,2,3,4 


External values for definition of mode for 
manipulated outputs. The values have the 
following interpretation: 


1: Auto-mode 
2: Man-mode 
3: E3-mode 
4: Bal-mode 


Table C-2 lists variables related to process variables, that can be configured as 
inputs to the on-line controller. 


Table C-2. Inputs to 3dMPC related to process variables. 


Name Value Description 

PV * (range) Process variables used for feedback. 

extRef1 (range) External set-points for process variables. Used 
in E1-mode. 

loPV (range) External values for low soft limits for process 
variables. Used in Ext-mode. 

hiPV (range) External values for high soft limits for process 
variables. Used in Ext-mode. 
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Table C-2. Inputs to 3dMPC related to process variables. (Continued) 


Name Value Description 


loDEV (range) External values for low soft limits for process 
variable deviations. Used in Ext-mode. 


hiDEV (range) External values for high soft limits for process 
variables deviations. Used in Ext-mode. 


extModePV 2,3,4,5 External values for definition of mode for 
process variables. The values have the 
following interpretation: 

2: E1-mode 

3: Auto-mode 

4: Predict-mode 

5:Off-mode 


Table C-3 lists variables related to feed-forward variables, that can be configured as 
inputs to the on-line controller. 


Table C-3. Inputs to 3dMPC related to feed-forward variables. 


Name Value Description 


FF (range) Feed-forward variables used for feed-forward 
from measurable disturbances. 


extModeFF 1,2,3 External values for definition of mode for feed- 
forward variables. The values have the 
following interpretation: 


1: Active-mode 
2: Predict-mode 
3:Off-mode 
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Table C-4 lists other variables, that can be configured as inputs to the on-line 


controller. 
Table C-4. Other inputs to 3d MPC. 
Name Value Description 
parSched (range) Scalar variable for parameter scheduling 
watchdogIn real Scalar variable for monitoring of on-line 
controller from the OCS/DCS. See variable 
watchDogOut in Table C-5. 


Table C-5 lists variables, that can be configured as outputs from the on-line 


controller. 


Table C-5. Outputs from 3d MPC. 


Name 


Value 


Description 


OUT * 


(range) 


Manipulated outputs 


watchdogOut 


real 


Scalar variable for monitoring of on-line 
controller from the OCS/DCS. 


The scalar input variable watchDoglIn is 
returned to the OCS/DCS. 


modeOUT 


1,2,3,4 


Vector with actual modes of manipulated 
outputs. 


modePV 


1,2,3,4,5 


Vector with actual modes of process variables. 


modeFF 


1,2,3 


Vector with actual modes of feed-forward 
variables. 


Status 


Controller status. An eight digit integer value, 
ABCDEFGH, with interpretation of the different 
digits in Table C-6. To extract the information 
in the Status variable in the OCS/DCS, the 
corresponding OCS/DCS object must be able 
to represent integer numbers between 0 and 
30000000 with full accuracy. 
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Table C-6. Interpretation of digits in 3dMPC status output. 


Digit 


Code 
{ 


Description 


No OUT in Auto or no PV in Auto. 


Some OUT in Auto and some PV in Auto. 


All OUT in Auto and some PV in Auto. 


Execution with delay < 0.25 tsamp. 


Execution with 0.25 tsamp < delay < 0.50 tsamp. 


Execution with 0.50 tsamp < delay < 1.00 tsamp. 


Execution with delay = 1.00 tsamp. 


OPC server OK. 


OPC server failure. 


Optimal solution achieved for quadratic optimization. 


Feasible solution not found for quadratic optimization. 


No signal failure. 


Signal failure for one or more PV. 


Signal failure for one or more FF. 


Signal failure for one or more PV and FF. 


No alarms 


Alarm for one or more PV. 


Alarm for one or more DEV. 


Alarm for one or more PV and DEV. 
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Table C-6. Interpretation of digits in 3dMPC status output. (Continued) 


Digit | Code Description 
G 0 No inconsistency for extLloPV and extHiPV. 
1 Inconsistency for extLoPV or extHiPV. 
2 Inconsistency for extLoDEV or extHiDEV. 
3 Both of above inconsistencies are satisfied. 
H 0 No inconsistency for ext*OUT. 
1 Inconsistency for extLoHardOUT or extHiHardOUT. 
2 Inconsistency for extLoSoftOUT or extHiSoftOUT. 
4 Inconsistency for extLoSpeedOUT or extHiSpeedOUT. 
3,5, Combination of inconsistencies for ext*OUT formed from 
6,7 the sum of the three above values for inconsistencies. 


C.1.1 OCS Connection Example 


246 


Normally the manipulated outputs from 3d MPC are external set points for 


controllers in the underlying OCS. Figure C-1 shows an example of how 3dMPC is 


connected to a cascaded PID-controller in the ABB Advant OCS with Master 
Software. Other OCS/DCS are assumed to be connected analogously. There is one 
such connection structure for each manipulated output. The signals are transferred 
using the OPC interface. 


The connection structure in Figure C-1 provides the following interaction between 
3dMPC and the PID-controller PIDCONA: 


If the PID controller PIDCONA is operating in El-mode it gets it set-point 
from 3d MPC. The extModeOUT(i) variable is then equal to 1 and does not 
force the i:th manipulated output to operate in any mode. 


The blockInc(i) and the blockDec(i) inputs become active when the PID 
controller output saturates. 
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3dMPC 
; PIDCONA 
blockDec(i) 
'_ pe} blockInc(i) OUT=LL |— 
OUT=HL]} —— 
OUT(i)} | EXTREF1 SW-C 
El ae | 1 
; —p 11 
|) extModeOUT(i) 1p» REVACT 4—e D> 
r balRef(i) 1 —y| BUMPLESS ee 


Figure C-1. Example Configuration with ABB Advant OCS with Master 
Software for the i:th manipulated output 


e — Ifthe PID controller PIDCONA is not operating in El-mode, the 
extModeOUT(i) variable is then equal to 4 and forces the i:th manipulated 
output to operate in bal-mode. The OUT(i) would then follow the WSP for the 
PID controller by using the balRef() input. 


° If REVACT=0, then the PID controller has direct action, and the connections to 
blockInc(i) and blockDec(i) have to be interchanged. 


¢ The BUMPLESS input in the PID controller assures that the WSP output from 
the PID controller is correctly updated when the PID controller is not running 
in automatic modes (set-point tracking is enforced during non-automatic 
modes). 
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C.2 Format of log-files 


Logged data is stored in ascii files. These contain a header before the logged data 
appears. The header contains the following information: 


Number of manipulated outputs. 
Number of process variables. 

Number of feed-forward variables. 
Names of manipulated outputs. 

Names of process variables. 

Names of feed-forward variables. 
Descriptive text that is given by the user. 


Each row then in the log-file consists of data described in Table C-7. The size values 
refer to the number of columns each vector occupies in the log file. 


Table C-7. Description of data in log-files. 


Data Size Description 
date 1 Actual date, for example 19970418 
time 1 Seconds elapsed since mid-night. 
Info 1 Definition of actual excitation. No excitation: 0. PRBS 
excitation: 10 
addOUT nOUT | Vector with actual values of the excitation signal. 
OUT nOUT | Vector with actual values of manipulated outputs. 
PV nPV Vector with actual values of process variables. 
WSP nPV Vector with actual values of working set-points. 
FF nFF Vector with actual values of feed forward variables. 


modeOUT |nOUT | Vector with actual modes of manipulated outputs. 


modePV nPV Vector with actual modes of process variables. 


modeFF nFF Vector with actual modes of feed forward variables. 
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P 

page selector toggle buttons, model evaluation 
tool 111 
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specModeOUT 136 

specModePV 134 

speed limits, manipulated outputs 137 

speed limits, set-point 135 

speed-ratio, observer design 140, 229 

split data 57, 209 

standard functions, non-linear transformation 
editor 65 

start controller 154 

start-menu 33 

start-up and shut-down 153 

state-space model 193 

static optimization 11 

status codes 244 

steady-state gain indication 111, 115 

step amplitude, MPC analysis tool 145 

step length, MPC analysis tool 145 

step response plots 117, 142 

step size button, model evaluation tool 111 

step sizes, default 101 

step-response lengths, default 102 

stop controller 154, 156 

subspace identification 193, 213 

support signals 221, 230 
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three degree of freedom 11, 133, 230 
time constant, user defined models 127 
time-stamps 208 

trends, remove 50, 209 

tuning window 162, 179 

tuning, MPC 130, 226 


U 


user defined models 125 


W 


warm start option 167 
web plot window 175 
webplot display window 161 


X 


x-interval button, model evaluation tool 110 


Y 


y-interval button, model evaluation tool 110 


Z 


zoom, cut data tool 57 

zoom, model evaluation tool 110 

zoom, remove outliers tool 68 

zoom, signal close-up in pre-processing tool 53 
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